Changeset 7293 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

Timestamp:

2016-11-21T16:27:53+01:00 (7 years ago)

Author:

vancop

Message:

Commit infrastructure for melt pond impletation

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90 (modified) (4 diffs)

branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

@@ -207 +207 @@
       CALL iom_put ('hfxspr'     , hfx_spr(:,:)         )   ! Heat content of snow precip 
 
+      ! MV MP 2016
+      IF ( ln_limMP ) THEN
+         CALL iom_put( "iceamp"  , at_ip  * zswi        )   ! melt pond total fraction
+         CALL iom_put( "icevmp"  , vt_ip  * zswi        )   ! melt pond total volume per unit area
+      ENDIF
+      ! END MV MP 2016
+
       
       !--------------------------------
@@ -223 +230 @@
       ! brine volume
       IF ( iom_use( "brinevol_cat" ) )  CALL iom_put( "brinevol_cat", bv_i * 100. * zswi2 )
+
+      ! MV MP 2016
+      IF ( ln_limMP ) THEN
+         IF ( iom_use( "iceamp_cat"  ) )  CALL iom_put( "iceamp_cat"     , a_ip   * zswi2   )       ! melt pond frac for categories
+         IF ( iom_use( "icevmp_cat"  ) )  CALL iom_put( "icevmp_cat"     , v_ip   * zswi2   )       ! melt pond frac for categories
+      ENDIF
+      ! END MV MP 2016
 
       !     !  Create an output files (output.lim.abort.nc) if S < 0 or u > 20 m/s
@@ -286 +300 @@
       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
 
+      ! MV MP 2016
+      IF ( ln_limMP ) THEN
+         CALL histdef( kid, "si_amp", "Melt pond fraction"      , "%"      ,   &
+      &         jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+         CALL histdef( kid, "si_vmp", "Melt pond volume"        ,  "m"     ,   &
+      &         jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      ENDIF
+      ! END MV MP 2016
+
       CALL histdef( kid, "vfxbog", "Ice bottom production"   , "m/s"    ,   &
       &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
@@ -330 +353 @@
       CALL histwrite( kid, "sidive", kt, divu_i*1.0e8   , jpi*jpj, (/1/) )
 
+      ! MV MP 2016
+      IF ( ln_limMP ) THEN
+         CALL histwrite( kid, "si_amp", kt, at_ip         , jpi*jpj, (/1/) )
+         CALL histwrite( kid, "si_vmp", kt, vt_ip         , jpi*jpj, (/1/) )
+      ENDIF
+      ! END MV MP 2016
+
       CALL histwrite( kid, "vfxbog", kt, wfx_bog        , jpi*jpj, (/1/) )
       CALL histwrite( kid, "vfxdyn", kt, wfx_dyn        , jpi*jpj, (/1/) )