Changeset 12401 for NEMO/branches/UKMO/NEMO_4.0.1_add_pond_lids/src/ICE/ice.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/ice.F90 (modified) (3 diffs)

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

@@ -70 +70 @@
    !! a_ip        |      -      |    Ice pond concentration       |       |
    !! v_ip        |      -      |    Ice pond volume per unit area| m     |
+   !! lh_ip       !    lh_ip_1d !    Ice pond lid thickness       ! m     !
    !!                                                                     |
    !!-------------|-------------|---------------------------------|-------|
@@ -332 +333 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   a_ip_frac  !: melt pond fraction (a_ip/a_i)
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   h_ip       !: melt pond depth                          [m]
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   lh_ip      !: melt pond lid thickness                  [m]
 
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   at_ip      !: total melt pond concentration
@@ -448 +450 @@
 
       ii = ii + 1
-      ALLOCATE( a_ip(jpi,jpj,jpl) , v_ip(jpi,jpj,jpl) , a_ip_frac(jpi,jpj,jpl) , h_ip(jpi,jpj,jpl) , STAT = ierr(ii) )
+      ALLOCATE( a_ip(jpi,jpj,jpl) , v_ip(jpi,jpj,jpl) , a_ip_frac(jpi,jpj,jpl) , h_ip(jpi,jpj,jpl), lh_ip(jpi,jpj,jpl) , STAT = ierr(ii) )
 
       ii = ii + 1