C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_cost_driver.F,v 1.16 2014/10/20 03:20:57 gforget Exp $ C $Name: $ #include "SEAICE_OPTIONS.h" #ifdef ALLOW_ECCO # include "ECCO_OPTIONS.h" #endif #ifdef ALLOW_CTRL # include "CTRL_OPTIONS.h" #endif #ifdef ALLOW_AUTODIFF # include "AUTODIFF_OPTIONS.h" #endif subroutine SEAICE_COST_DRIVER( myiter, mytime, mythid ) c ================================================================== c SUBROUTINE seaice_cost_driver c ================================================================== c c o Evaluate cost function contributions of c seaice observations c c ================================================================== c SUBROUTINE seaice_cost_driver c ================================================================== implicit none c == global variables == #include "EEPARAMS.h" #if (defined (ALLOW_ECCO) defined (ALLOW_SEAICE_COST_SMR_AREA)) #include "SIZE.h" #include "PARAMS.h" #include "GRID.h" #ifdef ALLOW_CAL # include "cal.h" #endif #ifdef ALLOW_CTRL # include "CTRL_SIZE.h" # include "ctrl.h" # include "ctrl_dummy.h" #endif #ifdef ALLOW_ECCO # include "ecco_cost.h" #endif #include "SEAICE_COST.h" #endif /* ALLOW_ECCO and ALLOW_SEAICE_COST_SMR_AREA */ c == routine arguments == integer myiter _RL mytime integer mythid #if (defined (ALLOW_ECCO) defined (ALLOW_SEAICE_COST_SMR_AREA)) c == local variables == character*(max_len_mbuf) msgbuf integer nnzbar integer nnzobs integer nrecloc _RL spminloc _RL spmaxloc _RL spzeroloc _RL localperiod c == end of interface == if (useSEAICE) then write(msgbuf,'(a)') 'ph-cost call cost_smrarea' call PRINT_MESSAGE( msgbuf, standardmessageunit, & SQUEEZE_RIGHT , mythid) nnzbar = 1 nnzobs = 1 nrecloc = ndaysrec spminloc = 0.0 spmaxloc = 1.0 spzeroloc = -1. localperiod = 86400. call SEAICE_COST_CONCENTRATION( & nnzbar, smrareabarfile, smrareabar, xx_smrareabar_mean_dummy, & nnzobs, smrareadatfile, smrareadat, mult_smrarea, & nrecloc, modelstartdate, localperiod, & maskC, wsmrarea, & spminloc, spmaxloc, spzeroloc, & objf_smrarea, num_smrarea, & myiter, mytime, mythid ) c read the area dat file and compare against the averaged sst file call SEAICE_COST_SST( & nnzbar, smrsstbarfile, smrsstbar, xx_smrsstbar_mean_dummy, & smrareabarfile, smrareabar, xx_smrareabar_mean_dummy, & nnzobs, smrareadatfile, smrareadat, mult_smrsst, & nrecloc, modelstartdate, localperiod, & maskC, wsmrarea, & spminloc, spmaxloc, spzeroloc, & objf_smrsst, num_smrsst, & myiter, mytime, mythid ) c read the area dat file and compare against the averaged salinity file call SEAICE_COST_SSS( & nnzbar, smrsssbarfile, smrsssbar, xx_smrsssbar_mean_dummy, & smrareabarfile, smrareabar, xx_smrareabar_mean_dummy, & nnzobs, smrareadatfile, smrareadat, mult_smrsss, & nrecloc, modelstartdate, localperiod, & maskC, wsmrarea, & spminloc, spmaxloc, spzeroloc, & objf_smrsss, num_smrsss, & myiter, mytime, mythid ) print *,' SICD cost smrarea/sst/sss ',objf_smrarea, & objf_smrsst,objf_smrsss print *,' SICD num points smrarea/sst/sss ',num_smrarea, & num_smrsst,num_smrsss # ifndef ALLOW_SEAICE_COST_AREASST c-- write(msgbuf,'(a)') 'ph-cost call cost_smrarea' c-- call print_message( msgbuf, standardmessageunit, c-- & SQUEEZE_RIGHT , mythid) c-- nnzbar = 1 c-- nnzobs = 1 c--cph( c-- nrecloc = ndaysrec c--cph hard-coding for ECCO c--cph nrecloc = 4383 c--cph) c-- spminloc = 0. c-- spmaxloc = 1. c-- spzeroloc = -9999. c-- localperiod = 86400. c-- call cost_generic( c-- & nnzbar, smrareabarfile, smrareabar, xx_smrareabar_mean_dummy, c-- & nnzobs, smrareadatfile, mult_smrarea, c-- & nrecloc, nrecloc, modelstartdate, localperiod, c-- & 'C', wsmrarea, spminloc, spmaxloc, spzeroloc, c-- & no_preproc, no_posproc, no_scalefile, c-- & 0, smrareadatfile, objf_smrarea, num_smrarea, c-- & myiter, mytime, mythid ) c-- write(msgbuf,'(a,1P2E18.10)') c-- & 'ph-cost smrarea ', objf_smrarea(1,1), num_smrarea(1,1) c-- call print_message( msgbuf, standardmessageunit, c-- & SQUEEZE_RIGHT , mythid) # endif /* ALLOW_SEAICE_COST_AREASST */ endif #endif /* ALLOW_ECCO and ALLOW_SEAICE_COST_SMR_AREA */ end