Changeset 2691 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jun 21, 2022, 11:05:31 AM (3 years ago)

Author:

aslmd

Message:

added loading of tracer profile for Fe in rcm1d.F

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F (modified) (1 diff)

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

@@ -681 +681 @@
             endif ! of if (txt.eq."h2o_ice")
 
+            !generic Fe_vap (will need to add more species after)
+            if (txt .eq. "generic_Fe_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Fe_vap" input file
+                  open(91,file='profile_Fe_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_Fe_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "generic_Fe_vap") 
+            !generic Fe_ice 
+            if (txt .eq. "generic_Fe_ice") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+            endif ! we only initialize the solid at 0
          enddo ! of do iq=1,nq