Changeset 1533 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Timestamp:

Apr 7, 2016, 5:25:33 PM (9 years ago)

Author:

mturbet

Message:

Tracers profiles now available for input

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -109 +109 @@
       character*2 str2
       character (len=7) :: str7
+      character(len=44) :: txt
 
       logical oldcompare, earthhack,saveprofile
@@ -122 +123 @@
       REAL cloudfrac(1,llm)
       REAL hice(1),totcloudfrac(1)
-      REAL qzero1D   !initial water amount on the ground
 
 !     added by BC for ocean
@@ -550 +550 @@
         qsurf(iq) = 0.
       ENDDO
-
-      if (water) then
-         qzero1D               = 0.0
-         qsurf(i_h2o_vap) = qzero1D
-      endif
+      
+      if (tracer) then ! Initialize tracers here. 
+             
+         write(*,*) "rcm1d : initializing tracers profiles"
+
+         do iq=1,nq
+         
+            txt=""
+            write(txt,"(a)") tname(iq)
+            write(*,*)"  tracer:",trim(txt)
+              
+            ! CO2
+            if (txt.eq."co2_ice") then
+               q(:,iq)=0.   ! kg/kg of atmosphere
+               qsurf(iq)=0. ! kg/m2 at the surface               
+               ! Look for a "profile_co2_ice" input file
+               open(91,file='profile_co2_ice',status='old',
+     &         form='formatted',iostat=ierr)
+               if (ierr.eq.0) then
+                  read(91,*) qsurf(iq)
+                  do ilayer=1,nlayer
+                     read(91,*) q(ilayer,iq)
+                  enddo
+               else
+                  write(*,*) "No profile_co2_ice file!"
+               endif
+               close(91)
+            endif ! of if (txt.eq."co2")
+          
+            ! WATER VAPOUR
+            if (txt.eq."h2o_vap") then
+               q(:,iq)=0.   ! kg/kg of atmosphere
+               qsurf(iq)=0. ! kg/m2 at the surface
+               ! Look for a "profile_h2o_vap" input file   
+               open(91,file='profile_h2o_vap',status='old',
+     &         form='formatted',iostat=ierr)
+               if (ierr.eq.0) then
+                  read(91,*) qsurf(iq)
+                  do ilayer=1,nlayer
+                     read(91,*) q(ilayer,iq)
+                  enddo
+               else
+                  write(*,*) "No profile_h2o_vap file!"
+               endif
+               close(91)
+            endif ! of if (txt.eq."h2o_vap")
+            
+            ! WATER ICE
+            if (txt.eq."h2o_ice") then
+               q(:,iq)=0.   ! kg/kg of atmosphere
+               qsurf(iq)=0. ! kg/m2 at the surface
+               ! Look for a "profile_h2o_ice" input file
+               open(91,file='profile_h2o_ice',status='old',
+     &         form='formatted',iostat=ierr)
+               if (ierr.eq.0) then
+                  read(91,*) qsurf(iq)
+                  do ilayer=1,nlayer
+                     read(91,*) q(ilayer,iq)
+                  enddo
+               else
+                  write(*,*) "No profile_h2o_ice file!"
+               endif
+               close(91)
+            endif ! of if (txt.eq."h2o_ice")
+
+         enddo ! of do iq=1,nq
+         
+      endif ! of tracer
 
 c   Initialisation pour prendre en compte les vents en 1-D