Changeset 1359 for trunk/LMDZ.GENERIC

Timestamp:

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

Author:

aslmd

Message:

LMDZ.GENERIC : added an option to allow to change preff and use topo. corrected a bug with isotherm and hydrostatic integrations.

File:

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

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

@@ -578 +578 @@
       write(*,*)
       write(*,*) 'flat : no topography ("aquaplanet")'
+      write(*,*) 'set_ps_to_preff : used if changing preff with topo'
       write(*,*) 'nuketharsis : no Tharsis bulge'
       write(*,*) 'bilball : uniform albedo and thermal inertia'
@@ -663 +664 @@
              enddo
          end if
+
+c       'set_ps_to_preff' : used if changing preff with topo  
+c       ----------------------------------------------------
+        else if (trim(modif) .eq. 'set_ps_to_preff') then
+          do j=1,jjp1
+           do i=1,iip1
+             ps(i,j)=preff
+           enddo
+          enddo
 
 c       'nuketharsis : no tharsis bulge for Early Mars'
@@ -1311 +1321 @@
           Tset(1:iip1,1:jjp1,1:llm)=Tiso
           flagtset=.true.
-          
+
+          t(1:iip1,1:jjp1,1:llm)=Tiso
+          !! otherwise hydrost. integrations below
+          !! use the wrong temperature
+          !! -- NB: Tset might be useless
+        
           ucov(1:iip1,1:jjp1,1:llm)=0
           vcov(1:iip1,1:jjm,1:llm)=0