Changeset 1783 for trunk

Timestamp:

Sep 14, 2017, 8:31:56 PM (7 years ago)

Author:

mturbet

Message:

add new option in newstart

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/dynphy_lonlat/phystd/newstart.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/README

@@ -1354 +1354 @@
  - Added a tmpk variable
  - Save ~ 50% time on the RT, ~30% on the whole code on the tested simulations
+
+== 07/06/2017 == MT
+Add qs=x option to newstart.F (chose a specific value to a surface tracer)

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

@@ -588 +588 @@
       write(*,*) 'q=0 : ALL tracer =zero'
       write(*,*) 'q=x : give a specific uniform value to one tracer'
+      write(*,*) 'qs=x : give a uniform value to a surface tracer'
       write(*,*) 'q=profile    : specify a profile for a tracer'
 !      write(*,*) 'ini_q : tracers initialisation for chemistry, water an
@@ -951 +952 @@
                ENDDO
              ENDDO
+             write(*,*) 'SURFACE value of tracer ',trim(tname(iq)),
+     &                   ' ? (kg/m2)'
+             read(*,*) val
+             DO ig=1,ngridmx
+                 qsurf(ig,iq)=val
+             ENDDO
+             
+c       qs=x : initialise surface tracer manually 
+c       --------------------------------
+        else if (trim(modif).eq.'qs=x') then
+             write(*,*) 'Which tracer do you want to modify ?'
+             do iq=1,nqtot
+               write(*,*)iq,' : ',trim(tname(iq))
+             enddo
+             write(*,*) '(choose between 1 and ',nqtot,')'
+             read(*,*) iq 
              write(*,*) 'SURFACE value of tracer ',trim(tname(iq)),
      &                   ' ? (kg/m2)'