Changeset 996 for trunk/LMDZ.GENERIC/libf/phystd/newstart.F

Timestamp:

Jun 21, 2013, 1:54:50 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.GIANT R.I.P. (but moved def files in LMDZ.UNIVERSAL for the record). LMDZ.GENERIC added a mode in newstart to initialize with a 1D profile. added Ls in outputs. LMDZ.UNIVERSAL added def for saturn_tropo_strato

File:

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

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

@@ -540 +540 @@
       write(*,*) 'coldspole : cold subsurface and high albedo at S.pole'
       write(*,*) 'qname : change tracer name'
+      write(*,*) 't=profile  : read temperature profile in profile.in'
       write(*,*) 'q=0 : ALL tracer =zero'
       write(*,*) 'q=x : give a specific uniform value to one tracer'
@@ -892 +893 @@
                  qsurf(ig,iq)=val
              ENDDO
+
+c       t=profile : initialize temperature with a given profile
+c       -------------------------------------------------------
+        else if (trim(modif) .eq. 't=profile') then
+             write(*,*) 'Temperature profile from ASCII file'
+             write(*,*) "'profile.in' e.g. 1D output"
+             write(*,*) "(one value per line in file; starting with"
+             write(*,*) "surface value, the 1st atmospheric layer"
+             write(*,*) "followed by 2nd, etc. up to top of atmosphere)"
+             txt="profile.in"
+             open(33,file=trim(txt),status='old',form='formatted',
+     &            iostat=ierr)
+             if (ierr.eq.0) then
+               ! OK, found file 'profile_...', load the profile
+               do l=1,llm+1
+                 read(33,*,iostat=ierr) profile(l)
+                 write(*,*) profile(l)
+                 if (ierr.ne.0) then ! something went wrong
+                   exit ! quit loop
+                 endif
+               enddo
+               if (ierr.eq.0) then
+                 tsurf(1:ngridmx)=profile(1)
+                 tsoil(1:ngridmx,1:nsoilmx)=profile(1)
+                 do l=1,llm
+                   Tset(1:iip1,1:jjp1,l)=profile(l+1)
+                   flagtset=.true.
+                 enddo
+                 ucov(1:iip1,1:jjp1,1:llm)=0.
+                 vcov(1:iip1,1:jjm,1:llm)=0.
+                 q2(1:ngridmx,1:nlayermx+1)=0.
+               else
+                 write(*,*)'problem reading file ',trim(txt),' !'
+                 write(*,*)'No modifications to temperature'
+               endif
+             else
+               write(*,*)'Could not find file ',trim(txt),' !'
+               write(*,*)'No modifications to temperature'
+             endif
 
 c       q=profile : initialize tracer with a given profile