Changeset 9118


Ignore:
Timestamp:
2017-12-18T16:17:44+01:00 (3 years ago)
Author:
clem
Message:

debug ice test cases with thermo or dynamics only

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r9076 r9118  
    235235      CALL ice_itd_init                ! ice thickness distribution initialization 
    236236      ! 
    237       IF( ln_icethd ) THEN 
    238          CALL ice_thd_init             ! set ice thermodynics parameters (clem: important to call it first for melt ponds) 
    239       ENDIF    
     237      CALL ice_thd_init                ! set ice thermodynics parameters (clem: important to call it first for melt ponds) 
     238      ! 
    240239      !                                ! Initial sea-ice state 
    241240      IF( .NOT. ln_rstart ) THEN              ! start from rest: sea-ice deduced from sst 
     
    250249      CALL ice_forcing_init            ! set ice-ocean and ice-atm. coupling parameters 
    251250      ! 
    252       IF( ln_icedyn ) THEN 
    253          CALL ice_dyn_init             ! set ice dynamics parameters 
    254       ENDIF 
     251      CALL ice_dyn_init                ! set ice dynamics parameters 
    255252      ! 
    256253      CALL ice_update_init             ! ice surface boundary condition 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

    r9019 r9118  
    124124         tm_si(:,:) = SUM( t_si(:,:,:) * a_i(:,:,:) , dim=3 ) * z1_at_i(:,:) 
    125125         om_i (:,:) = SUM( oa_i(:,:,:)              , dim=3 ) * z1_at_i(:,:) 
     126         sm_i (:,:) = SUM( sv_i(:,:,:)              , dim=3 ) * z1_vt_i(:,:) 
    126127         ! 
    127128         tm_i(:,:) = 0._wp 
    128          sm_i(:,:) = 0._wp 
    129129         DO jl = 1, jpl 
    130130            DO jk = 1, nlay_i 
    131131               tm_i(:,:) = tm_i(:,:) + r1_nlay_i * t_i (:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:) 
    132                sm_i(:,:) = sm_i(:,:) + r1_nlay_i * sz_i(:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:) 
    133132            END DO 
    134133         END DO 
     
    298297         !            !---------------------------------------! 
    299298         sz_i(:,:,:,:) = rn_icesal 
    300          s_i(:,:,:)   = rn_icesal 
     299         s_i (:,:,:)   = rn_icesal 
    301300         ! 
    302301         !            !---------------------------------------------! 
     
    313312         !                                      ! Slope of the linear profile  
    314313         WHERE( h_i(:,:,:) > epsi20 )   ;   z_slope_s(:,:,:) = 2._wp * s_i(:,:,:) / h_i(:,:,:) 
    315          ELSEWHERE                       ;   z_slope_s(:,:,:) = 0._wp 
     314         ELSEWHERE                      ;   z_slope_s(:,:,:) = 0._wp 
    316315         END WHERE 
    317316         ! 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90

    r9019 r9118  
    248248      CALL histvert( kid, "ncatice", "Ice Categories","", jpl, jcat, nz_i, "up") 
    249249 
    250       CALL histdef( kid, "sithic", "Ice thickness"           , "m"      ,   & 
    251       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    252       CALL histdef( kid, "siconc", "Ice concentration"       , "%"      ,   & 
    253       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    254       CALL histdef( kid, "sitemp", "Ice temperature"         , "C"      ,   & 
    255       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    256       CALL histdef( kid, "sivelu", "i-Ice speed "            , "m/s"    ,   & 
    257       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    258       CALL histdef( kid, "sivelv", "j-Ice speed "            , "m/s"    ,   & 
    259       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    260       CALL histdef( kid, "sistru", "i-Wind stress over ice " , "Pa"     ,   & 
    261       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    262       CALL histdef( kid, "sistrv", "j-Wind stress over ice " , "Pa"     ,   & 
    263       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    264       CALL histdef( kid, "sisflx", "Solar flux over ocean"     , "w/m2" ,   & 
    265       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    266       CALL histdef( kid, "sinflx", "Non-solar flux over ocean" , "w/m2" ,   & 
    267       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    268       CALL histdef( kid, "isnowpre", "Snow precipitation"      , "kg/m2/s",   & 
    269       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    270       CALL histdef( kid, "sisali", "Ice salinity"            , "PSU"    ,   & 
    271       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    272       CALL histdef( kid, "sivolu", "Ice volume"              , "m"      ,   & 
    273       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    274       CALL histdef( kid, "sidive", "Ice divergence"          , "10-8s-1",   & 
    275       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
    276       CALL histdef( kid, "si_amp", "Melt pond fraction"      , "%"      ,   & 
    277       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    278       CALL histdef( kid, "si_vmp", "Melt pond volume"        ,  "m"     ,   & 
    279       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    280       CALL histdef( kid, "vfxbog", "Ice bottom production"   , "m/s"    ,   & 
    281       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    282       CALL histdef( kid, "vfxdyn", "Ice dynamic production"  , "m/s"    ,   & 
    283       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    284       CALL histdef( kid, "vfxopw", "Ice open water prod"     , "m/s"    ,   & 
    285       &       jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    286       CALL histdef( kid, "vfxsni", "Snow ice production "    , "m/s"    ,   & 
    287       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    288       CALL histdef( kid, "vfxres", "Ice prod from corrections" , "m/s"  ,   & 
    289       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    290       CALL histdef( kid, "vfxbom", "Ice bottom melt"         , "m/s"    ,   & 
    291       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    292       CALL histdef( kid, "vfxsum", "Ice surface melt"        , "m/s"    ,   & 
    293       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    294       CALL histdef( kid, "vfxpnd", "Ice melt ponds flux"     , "m/s"    ,   & 
    295       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    296  
    297       CALL histdef( kid, "sithicat", "Ice thickness"         , "m"      ,   & 
    298       &      jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    299       CALL histdef( kid, "siconcat", "Ice concentration"     , "%"      ,   & 
    300       &      jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    301       CALL histdef( kid, "sisalcat", "Ice salinity"           , ""      ,   & 
    302       &      jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    303       CALL histdef( kid, "snthicat", "Snw thickness"         , "m"      ,   & 
    304       &      jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
     250      CALL histdef( kid, "sithic", "Ice thickness"          , "m"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     251      CALL histdef( kid, "siconc", "Ice concentration"      , "%"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     252      CALL histdef( kid, "sitemp", "Ice temperature"        , "C"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     253      CALL histdef( kid, "sivelu", "i-Ice speed "           , "m/s"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     254      CALL histdef( kid, "sivelv", "j-Ice speed "           , "m/s"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     255      CALL histdef( kid, "sistru", "i-Wind stress over ice" , "Pa"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     256      CALL histdef( kid, "sistrv", "j-Wind stress over ice" , "Pa"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     257      CALL histdef( kid, "sisflx", "Solar flx over ocean"   , "W/m2"   , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     258      CALL histdef( kid, "sinflx", "NonSolar flx over ocean", "W/m2"   , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     259      CALL histdef( kid, "snwpre", "Snow precipitation"     , "kg/m2/s", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     260      CALL histdef( kid, "sisali", "Ice salinity"           , "PSU"    , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     261      CALL histdef( kid, "sivolu", "Ice volume"             , "m"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     262      CALL histdef( kid, "sidive", "Ice divergence"         , "10-8s-1", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )  
     263      CALL histdef( kid, "si_amp", "Melt pond fraction"     , "%"      , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     264      CALL histdef( kid, "si_vmp", "Melt pond volume"       ,  "m"     , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     265      ! 
     266      CALL histdef( kid, "sithicat", "Ice thickness"        , "m"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
     267      CALL histdef( kid, "siconcat", "Ice concentration"    , "%"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
     268      CALL histdef( kid, "sisalcat", "Ice salinity"         , ""       , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
     269      CALL histdef( kid, "snthicat", "Snw thickness"        , "m"      , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 
    305270 
    306271      CALL histend( kid, snc4set )   ! end of the file definition 
    307272 
    308       CALL histwrite( kid, "sithic", kt, hm_i           , jpi*jpj, (/1/) )     
    309       CALL histwrite( kid, "siconc", kt, at_i           , jpi*jpj, (/1/) ) 
    310       CALL histwrite( kid, "sitemp", kt, tm_i - rt0     , jpi*jpj, (/1/) ) 
    311       CALL histwrite( kid, "sivelu", kt, u_ice          , jpi*jpj, (/1/) ) 
    312       CALL histwrite( kid, "sivelv", kt, v_ice          , jpi*jpj, (/1/) ) 
    313       CALL histwrite( kid, "sistru", kt, utau_ice       , jpi*jpj, (/1/) ) 
    314       CALL histwrite( kid, "sistrv", kt, vtau_ice       , jpi*jpj, (/1/) ) 
    315       CALL histwrite( kid, "sisflx", kt, qsr , jpi*jpj, (/1/) ) 
    316       CALL histwrite( kid, "sinflx", kt, qns , jpi*jpj, (/1/) ) 
    317       CALL histwrite( kid, "isnowpre", kt, sprecip        , jpi*jpj, (/1/) ) 
    318       CALL histwrite( kid, "sisali", kt, sm_i           , jpi*jpj, (/1/) ) 
    319       CALL histwrite( kid, "sivolu", kt, vt_i           , jpi*jpj, (/1/) ) 
    320       CALL histwrite( kid, "sidive", kt, divu_i*1.0e8   , jpi*jpj, (/1/) ) 
    321  
     273      CALL histwrite( kid, "sithic", kt, hm_i          , jpi*jpj, (/1/) )     
     274      CALL histwrite( kid, "siconc", kt, at_i          , jpi*jpj, (/1/) ) 
     275      CALL histwrite( kid, "sitemp", kt, tm_i - rt0    , jpi*jpj, (/1/) ) 
     276      CALL histwrite( kid, "sivelu", kt, u_ice         , jpi*jpj, (/1/) ) 
     277      CALL histwrite( kid, "sivelv", kt, v_ice         , jpi*jpj, (/1/) ) 
     278      CALL histwrite( kid, "sistru", kt, utau_ice      , jpi*jpj, (/1/) ) 
     279      CALL histwrite( kid, "sistrv", kt, vtau_ice      , jpi*jpj, (/1/) ) 
     280      CALL histwrite( kid, "sisflx", kt, qsr           , jpi*jpj, (/1/) ) 
     281      CALL histwrite( kid, "sinflx", kt, qns           , jpi*jpj, (/1/) ) 
     282      CALL histwrite( kid, "snwpre", kt, sprecip       , jpi*jpj, (/1/) ) 
     283      CALL histwrite( kid, "sisali", kt, sm_i          , jpi*jpj, (/1/) ) 
     284      CALL histwrite( kid, "sivolu", kt, vt_i          , jpi*jpj, (/1/) ) 
     285      CALL histwrite( kid, "sidive", kt, divu_i*1.0e8  , jpi*jpj, (/1/) ) 
    322286      CALL histwrite( kid, "si_amp", kt, at_ip         , jpi*jpj, (/1/) ) 
    323287      CALL histwrite( kid, "si_vmp", kt, vt_ip         , jpi*jpj, (/1/) ) 
    324  
    325       CALL histwrite( kid, "vfxbog", kt, wfx_bog        , jpi*jpj, (/1/) ) 
    326       CALL histwrite( kid, "vfxdyn", kt, wfx_dyn        , jpi*jpj, (/1/) ) 
    327       CALL histwrite( kid, "vfxopw", kt, wfx_opw        , jpi*jpj, (/1/) ) 
    328       CALL histwrite( kid, "vfxsni", kt, wfx_sni        , jpi*jpj, (/1/) ) 
    329       CALL histwrite( kid, "vfxres", kt, wfx_res        , jpi*jpj, (/1/) ) 
    330       CALL histwrite( kid, "vfxbom", kt, wfx_bom        , jpi*jpj, (/1/) ) 
    331       CALL histwrite( kid, "vfxsum", kt, wfx_sum        , jpi*jpj, (/1/) ) 
    332       CALL histwrite( kid, "vfxpnd", kt, wfx_pnd        , jpi*jpj, (/1/) ) 
    333  
     288      ! 
    334289      CALL histwrite( kid, "sithicat", kt, h_i         , jpi*jpj*jpl, (/1/) )     
    335290      CALL histwrite( kid, "siconcat", kt, a_i         , jpi*jpj*jpl, (/1/) )     
Note: See TracChangeset for help on using the changeset viewer.