Changeset 12401 for NEMO/branches/UKMO/NEMO_4.0.1_add_pond_lids/src/ICE/icerst.F90

Timestamp:

2020-02-18T16:21:31+01:00 (4 years ago)

Author:

dancopsey

Message:

Add melt pond lid code.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_add_pond_lids/src/ICE/icerst.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_add_pond_lids/src/ICE/icerst.F90

@@ -132 +132 @@
       CALL iom_rstput( iter, nitrst, numriw, 'a_ip' , a_ip  )
       CALL iom_rstput( iter, nitrst, numriw, 'v_ip' , v_ip  )
+      CALL iom_rstput( iter, nitrst, numriw, 'lh_ip', lh_ip )
       ! Snow enthalpy
       DO jk = 1, nlay_s 
@@ -245 +246 @@
             CALL iom_get( numrir, jpdom_autoglo, 'a_ip' , a_ip )
             CALL iom_get( numrir, jpdom_autoglo, 'v_ip' , v_ip )
+            CALL iom_get( numrir, jpdom_autoglo, 'lh_ip', lh_ip)
          ELSE                                     ! start from rest
             IF(lwp) WRITE(numout,*) '   ==>>   previous run without melt ponds output then set it to zero'
             a_ip(:,:,:) = 0._wp
             v_ip(:,:,:) = 0._wp
+            lh_ip(:,:,:) = 0._wp
          ENDIF
          ! fields needed for Met Office (Jules) coupling