C $Header: /u/gcmpack/MITgcm/pkg/streamice/streamice_dump_ad.F,v 1.2 2014/05/25 08:16:14 dgoldberg Exp $ C $Name: $ #include "AD_CONFIG.h" #include "PACKAGES_CONFIG.h" #include "STREAMICE_OPTIONS.h" CBOP C !ROUTINE: adstreamice_dump C !INTERFACE: subroutine ADSTREAMICE_DUMP( mytime, myiter, myThid ) C !DESCRIPTION: \bv C *==========================================================* C | SUBROUTINE adstreamice_dump | C *==========================================================* C Extract adjoint variable from TAMC/TAF-generated C adjoint common blocks, contained in adcommon.h C and write fields to file; C Make sure common blocks in adcommon.h are up-to-date C w.r.t. current adjoint code. C *==========================================================* C \ev C !USES: IMPLICIT NONE C == Global variables === #include "SIZE.h" #include "EEPARAMS.h" #include "PARAMS.h" #include "STREAMICE.h" #include "GRID.h" #ifdef ALLOW_AUTODIFF_MONITOR # include "adcommon.h" #endif LOGICAL DIFFERENT_MULTIPLE EXTERNAL INTEGER IO_ERRCOUNT EXTERNAL C !INPUT/OUTPUT PARAMETERS: C == Routine arguments == C myIter - iteration counter for this thread C myTime - time counter for this thread C myThid - Thread number for this instance of the routine. integer myThid integer myiter _RL mytime #if (defined (ALLOW_ADJOINT_RUN) defined (ALLOW_ADMTLM)) C !LOCAL VARIABLES: c == local variables == C suff - Hold suffix part of a filename C beginIOErrCount - Begin and end IO error counts C endIOErrCount C msgBuf - Error message buffer CHARACTER*(MAX_LEN_FNAM) suff INTEGER beginIOErrCount INTEGER endIOErrCount CHARACTER*(MAX_LEN_MBUF) msgBuf double precision :: area_shelf_streamice_ad(1-olx:snx+olx,1-oly: $sny+oly,nsx,nsy) double precision :: b_glen_ad(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) double precision :: bdot_streamice_ad(1-olx:snx+olx,1-oly:sny+oly, $nsx,nsy) double precision :: c_basal_friction_ad(1-olx:snx+olx,1-oly:sny+ $oly,nsx,nsy) double precision :: float_frac_streamice_ad(1-olx:snx+olx,1-oly: $sny+oly,nsx,nsy) double precision :: h_streamice_ad(1-olx:snx+olx,1-oly:sny+oly, $nsx,nsy) double precision :: surf_el_streamice_ad(1-olx:snx+olx,1-oly:sny+ $oly,nsx,nsy) double precision :: tau_beta_eff_streamice_ad(1-olx:snx+olx,1-oly: $sny+oly,nsx,nsy) double precision :: u_streamice_ad(1-olx:snx+olx,1-oly:sny+oly, $nsx,nsy) double precision :: v_streamice_ad(1-olx:snx+olx,1-oly:sny+oly, $nsx,nsy) double precision :: visc_streamice_ad(1-olx:snx+olx,1-oly:sny+oly, $nsx,nsy) common /streamice_fields_rl_ad/ h_streamice_ad, u_streamice_ad, $v_streamice_ad, visc_streamice_ad, tau_beta_eff_streamice_ad, $float_frac_streamice_ad, surf_el_streamice_ad, $area_shelf_streamice_ad, c_basal_friction_ad, b_glen_ad, $bdot_streamice_ad #ifdef USE_ALT_RLOW double precision :: r_low_si_ad(1-olx:snx+olx,1-oly:sny+oly,nsx, $nsy) common /streamice_rlow_ad/ r_low_si_ad #endif c == end of interface == CEOP IF ( & DIFFERENT_MULTIPLE(streamice_adjDump,mytime,deltaTClock) & ) THEN C-- Set suffix for this set of data files. WRITE(suff,'(I10.10)') myIter writeBinaryPrec = writeStatePrec C-- Read IO error counter beginIOErrCount = IO_ERRCOUNT(myThid) CALL WRITE_REC_3D_RL( & 'ADJc_basal_friction.'//suff, writeBinaryPrec, & 1, c_basal_friction_ad, 1, myIter, myThid ) CALL WRITE_REC_3D_RL( & 'ADJh_streamice.'//suff, writeBinaryPrec, & 1, h_streamice_ad, 1, myIter, myThid ) #ifdef USE_ALT_RLOW CALL WRITE_REC_3D_RL( & 'ADJr_low.'//suff, writeBinaryPrec, & 1, r_low_si_ad, 1, myIter, myThid ) #endif CALL WRITE_REC_3D_RL( & 'ADJb_glen.'//suff, writeBinaryPrec, & 1, b_glen_ad, 1, myIter, myThid ) CALL WRITE_REC_3D_RL( & 'ADJbdot.'//suff, writeBinaryPrec, & 1, bdot_streamice_ad, 1, myIter, myThid ) ENDIF #endif /* ALLOW_ADJOINT_RUN */ RETURN END