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

Timestamp:

Jun 21, 2022, 11:05:34 AM (2 years ago)

Author:

aslmd

Message:

added reading of Cr profile in rcm1d.def

File:

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

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

@@ -703 +703 @@
                   qsurf(iq) = 0. !kg/kg of atmosphere
             endif ! we only initialize the solid at 0
+            !generic Cr_vap (will need to add more species after)
+            if (txt .eq. "generic_Cr_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_Cr_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_Cr_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "generic_Cr_vap") 
+            !generic Cr_ice 
+            if (txt .eq. "generic_Cr_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