Changeset 996

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

Files:

• LMDZ.GENERIC/README (modified) (1 diff)
• LMDZ.GENERIC/libf/phystd/newstart.F (modified) (2 diffs)
• LMDZ.GENERIC/libf/phystd/physiq.F90 (modified) (3 diffs)
• LMDZ.GIANT (deleted)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/callphys.def (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/compile (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/diagfi.def (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/gases.def (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/run.def (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/saturn_start (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/traceur.def (added)
• LMDZ.UNIVERSAL/deftank/PROTO_saturn_tropostrato/z2sig.def (added)
• LMDZ.UNIVERSAL/deftank/old_tests_giants (copied) (copied from trunk/LMDZ.GIANT/deftank)
• LMDZ.UNIVERSAL/deftank/old_tests_giants/launch (modified) (2 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)

trunk/LMDZ.UNIVERSAL/deftank/old_tests_giants/launch

@@ -6 +6 @@
 #PBS -N  JUPITER 
 #PBS -l "nodes=1:ppn=8"
+#PBS -l "mem=16gb"
 ##PBS -l "nodes=2:ppn=8"
 ##PBS -l "nodes=2:ppn=12"
@@ -33 +34 @@
 #mpirun wrf.exe
 
-mpirun gcm_80x120x15_phylmd_para.e 
-#mpirun gcm_160x240x15_phylmd_para.e
+ulimit -s unlimited
+
+##mpirun gcm_80x120x15_phylmd_para.e 
+#mpirun gcm_48x36x19_phynophys_para.e
+##mpirun gcm_160x240x15_phylmd_para.e
+
+mpirun gcm_80x120x15_phynophys_para.e
 
 ################