Changeset 996 for trunk/LMDZ.GENERIC

Timestamp:

Jun 21, 2013, 1:54:50 AM (12 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

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/newstart.F (modified) (2 diffs)
• libf/phystd/physiq.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/README

@@ -975 +975 @@
  - removed obsolete addit.F (and change corresponding lines in gcm)
  - remove unused "description.h" (and many places where it was "included")
+
+== 21/06/2013 == AS
+- Added possibility to initialize temperature from a rcm1d profile (t=profile)
+- Added output for Ls in diagfi

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

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -411 +411 @@
       integer num_run
       logical,save :: ice_update
-      
+
 !=======================================================================
 
@@ -1548 +1548 @@
       endif
 
+
          if (cursor == 1) then
            print*,'--> Ls =',zls*180./pi
@@ -1677 +1678 @@
 !        but its preferable to keep all the calls in one place...
 
+        call writediagfi(ngrid,"Ls","solar longitude","deg",0,zls*180./pi)
         call writediagfi(ngrid,"tsurf","Surface temperature","K",2,tsurf)
         call writediagfi(ngrid,"ps","Surface pressure","Pa",2,ps)