Ignore:
Timestamp:
2017-05-29T16:58:42+02:00 (3 years ago)
Author:
vancop
Message:

More on melt ponds in LIM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

    r8063 r8085  
    171171      CALL iom_put( "vfxice"     , wfx_ice * ztmp       )        ! total ice growth/melt  
    172172 
    173       IF ( nn_pnd_scheme >= 0 ) & 
     173      IF ( ln_pnd .AND. ( nn_pnd_scheme >= 0 ) ) & 
    174174         CALL iom_put( "vfxpnd"  , wfx_pnd * ztmp       )        ! melt pond water flux 
    175175 
     
    211211 
    212212      ! MV MP 2016 
    213       IF ( nn_pnd_scheme >= 0 ) THEN ! change once we are sure 
     213      IF ( ln_pnd ) THEN 
    214214         CALL iom_put( "iceamp"  , at_ip  * zswi        )   ! melt pond total fraction 
    215215         CALL iom_put( "icevmp"  , vt_ip  * zswi        )   ! melt pond total volume per unit area 
     
    235235 
    236236      ! MV MP 2016 
    237       IF ( nn_pnd_scheme >= 0 ) THEN ! change once we are sure that case = 0 works 
     237      IF ( ln_pnd .AND. ( nn_pnd_scheme >= 0 ) ) THEN 
    238238         IF ( iom_use( "iceamp_cat"  ) )  CALL iom_put( "iceamp_cat"     , a_ip   * zswi2   )       ! melt pond frac for categories 
    239239         IF ( iom_use( "icevmp_cat"  ) )  CALL iom_put( "icevmp_cat"     , v_ip   * zswi2   )       ! melt pond frac for categories 
     
    304304 
    305305      ! MV MP 2016 
    306       IF ( nn_pnd_scheme >= 0 ) THEN 
     306      IF ( ln_pnd .AND. ( nn_pnd_scheme >= 0 ) ) THEN 
    307307         CALL histdef( kid, "si_amp", "Melt pond fraction"      , "%"      ,   & 
    308308      &         jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     
    357357 
    358358      ! MV MP 2016 
    359       IF ( nn_pnd_scheme >= 0 ) THEN 
     359      IF ( ln_pnd .AND. ( nn_pnd_scheme >= 0 ) ) THEN 
    360360         CALL histwrite( kid, "si_amp", kt, at_ip         , jpi*jpj, (/1/) ) 
    361361         CALL histwrite( kid, "si_vmp", kt, vt_ip         , jpi*jpj, (/1/) ) 
     
    370370      CALL histwrite( kid, "vfxbom", kt, wfx_bom        , jpi*jpj, (/1/) ) 
    371371      CALL histwrite( kid, "vfxsum", kt, wfx_sum        , jpi*jpj, (/1/) ) 
    372       IF ( nn_pnd_scheme >= 0 ) & 
     372      IF ( ln_pnd .AND. ( nn_pnd_scheme >= 0 ) ) & 
    373373         CALL histwrite( kid, "vfxpnd", kt, wfx_pnd     , jpi*jpj, (/1/) ) 
    374374 
Note: See TracChangeset for help on using the changeset viewer.