Changeset 1844 for trunk/LMDZ.TITAN/libf/phytitan

Timestamp:

Dec 7, 2017, 10:15:41 AM (7 years ago)

Author:

jvatant

Message:

Enable tracers input profiles in 1D
--JVO

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• dyn1d/rcm1d.F (modified) (3 diffs)
• physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.TITAN/libf/phytitan/dyn1d/rcm1d.F

@@ -6 +6 @@
       use mod_grid_phy_lmdz, only : regular_lonlat
       use infotrac, only: nqtot, tname
-      use tracer_h, only: initracer2
+      use tracer_h
       use surfdat_h, only: albedodat, phisfi,
      &                     zmea, zstd, zsig, zgam, zthe,
@@ -122 +122 @@
 !     added by JVO
       REAL tankCH4(1)
+
+      CHARACTER*20 file_prof
+      LOGICAL findprof
 
 !     added by AS to avoid the use of adv trac common
@@ -732 +735 @@
 ! ------------------
       if(tracer) then
-        call initracer2(nq,tname) ! We need tracers infos in physdem1
-      endif
+        call initracer2(nq,tname) ! We need tracers names for physdem1 and also to read profiles
+        
+        ! Inquire if we have tracers input profiles
+        
+        do iq=1,nq
+        
+          file_prof = "profile_"//trim(noms(iq))
+
+          INQUIRE(file=trim(file_prof),exist=findprof)
+        
+          IF (findprof) THEN
+            ! Read profile for tracer. First line is surface value.
+            PRINT *, "Input file found for tracer ", noms(iq)
+            OPEN(11,file=trim(file_prof),status='old',form='formatted')
+            PRINT *, "---> reading tracer profile from input file ..."
+            READ(11,*) qsurf(1)
+            DO ilayer=1,nlayer
+              READ (11,*) q(ilayer,iq)
+            ENDDO
+          ELSE
+            PRINT *, "No input profile found for ", trim(noms(iq))
+            PRINT *, "---> we set it to zero !!"
+          ENDIF
+
+        enddo ! nq loop
+        
+      endif ! end if tracer
 
 c  Write a "startfi" file

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -408 +408 @@
 ! --------------------------------
       if (firstcall) then
-        ! initracer:
-        !   initialize nmicro, 
-        call initracer2(nq,nametrac)
-
-        ! ----------------------------------------------------------------------------
+
+        ! Initialisation of nmicro as well as tracers names, indexes ...
+        if (ngrid.ne.1) then
+           call initracer2(nq,nametrac) ! Already done in rcm1d
+        endif
 
         ! Allocate saved arrays.