Changeset 1369 for trunk/LMDZ.GENERIC

Timestamp:

Nov 19, 2014, 5:16:57 PM (10 years ago)

Author:

jbmadeleine

Message:

Added supercontinent option in newstart for generic model

File:

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

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

@@ -153 +153 @@
       character*1 yes
       logical :: flagtset=.false. ,  flagps0=.false.
-      real val, val2, val3 ! to store temporary variables
+      real val, val2, val3, val4 ! to store temporary variables
       real :: iceith=2000 ! thermal inertia of subterranean ice
       integer iref,jref
@@ -599 +599 @@
       write(*,*) 'oborealis : H2O ice across Vastitas Borealis'
       write(*,*) 'iceball   : Thick ice layer all over surface'
+      write(*,*) 'supercontinent: Create a continent of given Ab and TI'
       write(*,*) 'wetstart  : start with a wet atmosphere'
       write(*,*) 'isotherm  : Isothermal Temperatures, wind set to zero'
@@ -1305 +1306 @@
           CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi)
 
+c       supercontinent : H2O ice everywhere
+c       ----------------------------
+        else if (trim(modif) .eq. 'supercontinent') then
+          write(*,*) 'Minimum longitude (-180,180)'
+          read(*,*) val
+          write(*,*) 'Maximum longitude (-180,180)'
+          read(*,*) val2
+          write(*,*) 'Minimum latitude (-90,90)'
+          read(*,*) val3
+          write(*,*) 'Maximum latitude (-90,90)'
+          read(*,*) val4
+
+          do j=1,jjp1
+            do i=1,iip1
+              ig=1+(j-2)*iim +i
+              if(j.eq.1) ig=1
+              if(j.eq.jjp1) ig=ngridmx
+
+c             Supercontinent:
+              if (((rlatu(j)*180./pi.lt.val4).and.
+     &            (rlatu(j)*180./pi.gt.val3).and.
+     &            (rlonv(i)*180./pi.lt.val2).and.
+     &            (rlonv(i)*180./pi.gt.val))) then
+
+                rnat(ig)=1.
+                alb(i,j) = 0.3
+                do isoil=1,nsoilmx
+                  ith(i,j,isoil) = 2000.
+                enddo
+c             Ocean:
+              else
+                rnat(ig)=0.
+                alb(i,j) = 0.07
+                do isoil=1,nsoilmx
+                  ith(i,j,isoil) = 18000.
+                enddo
+              end if
+
+            enddo
+          enddo
+          CALL gr_dyn_fi(nsoilmx,iip1,jjp1,ngridmx,ith,ithfi)
+          CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi)
+
 c       isotherm : Isothermal temperatures and no winds
 c       -----------------------------------------------