Changeset 563 for trunk/LMDZ.MARS/libf/dyn3d/newstart.F

Timestamp:

Mar 6, 2012, 11:36:52 AM (13 years ago)

Author:

emillour

Message:

Mars GCM:

Added option "q=profile" in newstart to initialize a given tracer with

a profile specified in an ASCII file.

Minor fix on lect_start_archive : when tracers were not found, a default

value was asked to user, but twice! Once is enough.

EM

File:

• trunk/LMDZ.MARS/libf/dyn3d/newstart.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/dyn3d/newstart.F

@@ -155 +155 @@
       
       INTEGER :: nq,numvanle
-      character(len=20) :: txt ! to store some text
+      character(len=50) :: txt ! to store some text
       integer :: count
+      real :: profile(llm+1) ! to store an atmospheric profile + surface value
 
 ! MONS data:
@@ -482 +483 @@
       write(*,*) 'q=0 : ALL tracer =zero'
       write(*,*) 'q=x : give a specific uniform value to one tracer'
+      write(*,*) 'q=profile : specify a profile for a tracer'
       write(*,*) 'ini_q : tracers initialisation for chemistry, water an
      $d ice   '
@@ -779 +781 @@
                  qsurf(ig,iq)=val
              ENDDO
+
+c       q=profile : initialize tracer with a given profile
+c       --------------------------------------------------
+        else if (trim(modif) .eq. 'q=profile') then
+             write(*,*) 'Tracer profile will be sought in ASCII file'
+             write(*,*) "'profile_tracer' where 'tracer' is tracer name"
+             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)"
+             write(*,*) 'Which tracer do you want to set?'
+             do iq=1,nqmx
+               write(*,*)iq,' : ',trim(tnom(iq))
+             enddo
+             write(*,*) '(choose between 1 and ',nqmx,')'
+             read(*,*) iq 
+             if ((iq.lt.1).or.(iq.gt.nqmx)) then
+               ! wrong value for iq, go back to menu
+               write(*,*) "wrong input value:",iq
+               cycle
+             endif
+             ! look for input file 'profile_tracer'
+             txt="profile_"//trim(tnom(iq))
+             open(41,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(41,*,iostat=ierr) profile(l)
+                 if (ierr.ne.0) then ! something went wrong
+                   exit ! quit loop
+                 endif
+               enddo
+               if (ierr.eq.0) then
+                 ! initialize tracer values
+                 qsurf(:,iq)=profile(1)
+                 do l=1,llm
+                   q(:,:,l,iq)=profile(l+1)
+                 enddo
+                 write(*,*)'OK, tracer ',trim(tnom(iq)),' initialized ',
+     &                     'using values from file ',trim(txt)
+               else
+                 write(*,*)'problem reading file ',trim(txt),' !'
+                 write(*,*)'No modifications to tracer ',trim(tnom(iq))
+               endif
+             else
+               write(*,*)'Could not find file ',trim(txt),' !'
+               write(*,*)'No modifications to tracer ',trim(tnom(iq))
+             endif
+             
 
 c       ini_q : Initialize tracers for chemistry