Changeset 367 for trunk/LMDZ.GENERIC/libf/dyn3d

Timestamp:

Nov 9, 2011, 3:51:31 PM (13 years ago)

Author:

rwordsworth

Message:

ice initialization modified to make 33 kg/m2 everywhere at beginning.

File:

• trunk/LMDZ.GENERIC/libf/dyn3d/newstart.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/dyn3d/newstart.F

@@ -987 +987 @@
 
                  if (rlatu(j)*180./pi.gt.80.) then
-                    do isoil=1,nsoilmx
-                       ith(i,j,isoil) = 18000. ! thermal inertia
-                    enddo
+                    qsurf(ig,igcm_h2o_ice)=3.4e3
+                    !do isoil=1,nsoilmx
+                    !   ith(i,j,isoil) = 18000. ! thermal inertia
+                    !enddo
                    write(*,*)'     ==> Ice mesh North boundary (deg)= ',
      &                   rlatv(j-1)*180./pi
@@ -1028 +1029 @@
 
                  if (rlatu(j)*180./pi.lt.-80.) then
-                    do isoil=1,nsoilmx
-                       ith(i,j,isoil) = 18000. ! thermal inertia
-                    enddo
+                    qsurf(ig,igcm_h2o_ice)=3.4e3
+                    !do isoil=1,nsoilmx
+                    !   ith(i,j,isoil) = 18000. ! thermal inertia
+                    !enddo
                    write(*,*)'     ==> Ice mesh South boundary (deg)= ',
      &                   rlatv(j-1)*180./pi
@@ -1128 +1130 @@
 
                 if(phis(i,j).lt.-4000.*g)then
-                   qsurf(ig,igcm_h2o_ice)=1.e3
+                   !qsurf(ig,igcm_h2o_ice)=1.e3
+                   qsurf(ig,igcm_h2o_ice)=241.4 ! to make total 33 kg m^-2
                 endif
              ENDDO