Changeset 6351


Ignore:
Timestamp:
2016-02-24T19:50:11+01:00 (5 years ago)
Author:
cetlod
Message:

trunk: move the output of scale factor before time swapping and output some variables needs for offline, see ticket #1682

Location:
trunk/NEMOGCM
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/SHARED/field_def.xml

    r6140 r6351  
    2323      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    2424         <field id="e3t"          long_name="T-cell thickness"   standard_name="cell_thickness"   unit="m"   grid_ref="grid_T_3D"/> 
     25         <field id="e3t_0"        long_name="Initial T-cell thickness"   standard_name="ref_cell_thickness"   unit="m"   grid_ref="grid_T_3D"/> 
    2526 
    2627         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
     
    181182      <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc --> 
    182183         <field id="empmr"        long_name="Net Upward Water Flux"                standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   /> 
     184         <field id="empbmr"       long_name="Net Upward Water Flux at pre. tstep"  standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   /> 
    183185         <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="1e-3/m2/s" /> 
    184186         <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   /> 
     
    466468        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" /> 
    467469 
    468         <!-- woce_eiv: available with key_trabbl_adv --> 
     470 
    469471        <field id="avt"          long_name="vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     472        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    470473        <field id="avm"          long_name="vertical eddy viscosity"     standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
    471474 
    472475        <!-- avs: available with key_zdfddm --> 
    473476        <field id="avs"          long_name="salt vertical eddy diffusivity"   standard_name="ocean_vertical_salt_diffusivity"   unit="m2/s" /> 
     477        <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    474478 
    475479        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r6140 r6351  
    140140      ENDIF 
    141141 
    142       IF( ln_linssh ) THEN 
    143          CALL iom_put( "e3t" , e3t_n(:,:,:) ) 
    144          CALL iom_put( "e3u" , e3u_n(:,:,:) ) 
    145          CALL iom_put( "e3v" , e3v_n(:,:,:) ) 
    146          CALL iom_put( "e3w" , e3w_n(:,:,:) ) 
    147       ENDIF 
     142      ! Output of initial vertical scale factor 
     143      CALL iom_put("e3t_0", e3t_0(:,:,:) ) 
     144      CALL iom_put("e3u_0", e3t_0(:,:,:) ) 
     145      CALL iom_put("e3v_0", e3t_0(:,:,:) ) 
     146      ! 
     147      CALL iom_put( "e3t" , fse3t_n(:,:,:) ) 
     148      CALL iom_put( "e3u" , fse3u_n(:,:,:) ) 
     149      CALL iom_put( "e3v" , fse3v_n(:,:,:) ) 
     150      CALL iom_put( "e3w" , fse3w_n(:,:,:) ) 
     151      IF( iom_use("e3tdef") )   & 
     152         CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 ) 
    148153 
    149154      CALL iom_put( "ssh" , sshn )                 ! sea surface height 
    150       if( iom_use('ssh2') )   CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) )   ! square of sea surface height 
    151155       
    152156      CALL iom_put( "toce", tsn(:,:,:,jp_tem) )    ! 3D temperature 
     
    228232      CALL iom_put( "avm" , avmu                       )    ! T vert. eddy visc. coef. 
    229233      CALL iom_put( "avs" , fsavs(:,:,:)               )    ! S vert. eddy diff. coef. (useful only with key_zdfddm) 
     234 
     235      IF( iom_use('logavt') )   CALL iom_put( "logavt", LOG( MAX( 1.e-20_wp, avt  (:,:,:) ) ) ) 
     236      IF( iom_use('logavs') )   CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, fsavs(:,:,:) ) ) ) 
    230237 
    231238      IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN 
     
    292299         CALL iom_put( "eken", rke )            
    293300      ENDIF 
    294           
     301      ! 
     302      CALL iom_put( "hdiv", hdivn )                  ! Horizontal divergence 
     303      ! 
    295304      IF( iom_use("u_masstr") .OR. iom_use("u_heattr") .OR. iom_use("u_salttr") ) THEN 
    296305         z3d(:,:,jpk) = 0.e0 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r6152 r6351  
    653653         ht_n(:,:) = ht_n(:,:) + e3t_n(:,:,jk) * tmask(:,:,jk) 
    654654      END DO 
    655  
    656       ! Write outputs 
    657       ! ============= 
    658       CALL iom_put(     "e3t",   e3t_n(:,:,:) ) 
    659       CALL iom_put(     "e3u",   e3u_n(:,:,:) ) 
    660       CALL iom_put(     "e3v",   e3v_n(:,:,:) ) 
    661       CALL iom_put(     "e3w",   e3w_n(:,:,:) ) 
    662       CALL iom_put( "tpt_dep", gde3w_n(:,:,:) ) 
    663       IF( iom_use("e3tdef") )   & 
    664          CALL iom_put( "e3tdef", ( ( e3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100. * tmask(:,:,:) ) ** 2 ) 
    665655 
    666656      ! write restart file 
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r6140 r6351  
    430430      !                                                ! ---------------------------------------- ! 
    431431      IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN 
    432          CALL iom_put( "empmr" , emp  - rnf )                   ! upward water flux 
     432         CALL iom_put( "empmr"  , emp    - rnf )                ! upward water flux 
     433         CALL iom_put( "empbmr" , emp_b  - rnf )                ! before upward water flux ( needed to recalculate the time evolution of ssh in offline ) 
    433434         CALL iom_put( "saltflx", sfx  )                        ! downward salt flux   
    434435                                                                ! (includes virtual salt flux beneath ice  
Note: See TracChangeset for help on using the changeset viewer.