Ignore:
Timestamp:
2018-11-25T15:24:21+01:00 (22 months ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 5b: by default, suppress global communication in stpctl, see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icewri.F90

    r10314 r10358  
    227227 
    228228  
    229    SUBROUTINE ice_wri_state( kt, kid, kh_i ) 
     229   SUBROUTINE ice_wri_state( kid ) 
    230230      !!--------------------------------------------------------------------- 
    231231      !!                 ***  ROUTINE ice_wri_state  *** 
     
    238238      !! History :   4.0  !  2013-06  (C. Rousset) 
    239239      !!---------------------------------------------------------------------- 
    240       INTEGER, INTENT( in ) ::   kt               ! ocean time-step index 
    241       INTEGER, INTENT( in ) ::   kid , kh_i 
    242       INTEGER               ::   nz_i, jl 
    243       REAL(wp), DIMENSION(jpl) ::   jcat 
     240      INTEGER, INTENT( in ) ::   kid  
    244241      !!---------------------------------------------------------------------- 
    245242      ! 
    246       DO jl = 1, jpl 
    247          jcat(jl) = REAL(jl) 
    248       END DO 
    249        
    250       CALL histvert( kid, "ncatice", "Ice Categories","", jpl, jcat, nz_i, "up") 
    251  
    252       CALL histdef( kid, "sithic", "Ice thickness"          , "m"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    253       CALL histdef( kid, "siconc", "Ice concentration"      , "%"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    254       CALL histdef( kid, "sitemp", "Ice temperature"        , "C"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    255       CALL histdef( kid, "sivelu", "i-Ice speed "           , "m/s"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    256       CALL histdef( kid, "sivelv", "j-Ice speed "           , "m/s"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    257       CALL histdef( kid, "sistru", "i-Wind stress over ice" , "Pa"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    258       CALL histdef( kid, "sistrv", "j-Wind stress over ice" , "Pa"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    259       CALL histdef( kid, "sisflx", "Solar flx over ocean"   , "W/m2"   , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    260       CALL histdef( kid, "sinflx", "NonSolar flx over ocean", "W/m2"   , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    261       CALL histdef( kid, "snwpre", "Snow precipitation"     , "kg/m2/s", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    262       CALL histdef( kid, "sisali", "Ice salinity"           , "PSU"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    263       CALL histdef( kid, "sivolu", "Ice volume"             , "m"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    264       CALL histdef( kid, "sidive", "Ice divergence"         , "10-8s-1", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    265       CALL histdef( kid, "si_amp", "Melt pond fraction"     , "%"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    266       CALL histdef( kid, "si_vmp", "Melt pond volume"       ,  "m"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    267       ! 
    268       CALL histdef( kid, "sithicat", "Ice thickness"        , "m"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    269       CALL histdef( kid, "siconcat", "Ice concentration"    , "%"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    270       CALL histdef( kid, "sisalcat", "Ice salinity"         , ""       , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    271       CALL histdef( kid, "snthicat", "Snw thickness"        , "m"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    272  
    273       CALL histend( kid, snc4set )   ! end of the file definition 
    274  
    275       CALL histwrite( kid, "sithic", kt, hm_i          , jpi*jpj, (/1/) )     
    276       CALL histwrite( kid, "siconc", kt, at_i          , jpi*jpj, (/1/) ) 
    277       CALL histwrite( kid, "sitemp", kt, tm_i - rt0    , jpi*jpj, (/1/) ) 
    278       CALL histwrite( kid, "sivelu", kt, u_ice         , jpi*jpj, (/1/) ) 
    279       CALL histwrite( kid, "sivelv", kt, v_ice         , jpi*jpj, (/1/) ) 
    280       CALL histwrite( kid, "sistru", kt, utau_ice      , jpi*jpj, (/1/) ) 
    281       CALL histwrite( kid, "sistrv", kt, vtau_ice      , jpi*jpj, (/1/) ) 
    282       CALL histwrite( kid, "sisflx", kt, qsr           , jpi*jpj, (/1/) ) 
    283       CALL histwrite( kid, "sinflx", kt, qns           , jpi*jpj, (/1/) ) 
    284       CALL histwrite( kid, "snwpre", kt, sprecip       , jpi*jpj, (/1/) ) 
    285       CALL histwrite( kid, "sisali", kt, sm_i          , jpi*jpj, (/1/) ) 
    286       CALL histwrite( kid, "sivolu", kt, vt_i          , jpi*jpj, (/1/) ) 
    287       CALL histwrite( kid, "sidive", kt, divu_i*1.0e8  , jpi*jpj, (/1/) ) 
    288       CALL histwrite( kid, "si_amp", kt, at_ip         , jpi*jpj, (/1/) ) 
    289       CALL histwrite( kid, "si_vmp", kt, vt_ip         , jpi*jpj, (/1/) ) 
    290       ! 
    291       CALL histwrite( kid, "sithicat", kt, h_i         , jpi*jpj*jpl, (/1/) )     
    292       CALL histwrite( kid, "siconcat", kt, a_i         , jpi*jpj*jpl, (/1/) )     
    293       CALL histwrite( kid, "sisalcat", kt, s_i         , jpi*jpj*jpl, (/1/) )     
    294       CALL histwrite( kid, "snthicat", kt, h_s         , jpi*jpj*jpl, (/1/) )     
    295  
    296       !! The file is closed in dia_wri_state (ocean routine) 
    297       !! CALL histclo( kid ) 
    298       ! 
     243      !! The file is open in dia_wri_state (ocean routine) 
     244 
     245      CALL iom_rstput( 0, 0, kid, 'sithic', hm_i         )   ! Ice thickness 
     246      CALL iom_rstput( 0, 0, kid, 'siconc', at_i         )   ! Ice concentration 
     247      CALL iom_rstput( 0, 0, kid, 'sitemp', tm_i - rt0   )   ! Ice temperature 
     248      CALL iom_rstput( 0, 0, kid, 'sivelu', u_ice        )   ! i-Ice speed 
     249      CALL iom_rstput( 0, 0, kid, 'sivelv', v_ice        )   ! j-Ice speed 
     250      CALL iom_rstput( 0, 0, kid, 'sistru', utau_ice     )   ! i-Wind stress over ice 
     251      CALL iom_rstput( 0, 0, kid, 'sistrv', vtau_ice     )   ! i-Wind stress over ice 
     252      CALL iom_rstput( 0, 0, kid, 'sisflx', qsr          )   ! Solar flx over ocean 
     253      CALL iom_rstput( 0, 0, kid, 'sinflx', qns          )   ! NonSolar flx over ocean 
     254      CALL iom_rstput( 0, 0, kid, 'snwpre', sprecip      )   ! Snow precipitation 
     255      CALL iom_rstput( 0, 0, kid, 'sisali', sm_i         )   ! Ice salinity 
     256      CALL iom_rstput( 0, 0, kid, 'sivolu', vt_i         )   ! Ice volume 
     257      CALL iom_rstput( 0, 0, kid, 'sidive', divu_i*1.0e8 )   ! Ice divergence 
     258      CALL iom_rstput( 0, 0, kid, 'si_amp', at_ip        )   ! Melt pond fraction 
     259      CALL iom_rstput( 0, 0, kid, 'si_vmp', vt_ip        )   ! Melt pond volume 
     260      CALL iom_rstput( 0, 0, kid, 'sithicat', h_i        )   ! Ice thickness 
     261      CALL iom_rstput( 0, 0, kid, 'siconcat', a_i        )   ! Ice concentration 
     262      CALL iom_rstput( 0, 0, kid, 'sisalcat', s_i        )   ! Ice salinity 
     263      CALL iom_rstput( 0, 0, kid, 'snthicat', h_s        )   ! Snw thickness 
     264 
    299265    END SUBROUTINE ice_wri_state 
    300266 
Note: See TracChangeset for help on using the changeset viewer.