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

Timestamp:

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

Author:

aslmd

Message:

added loading of profile for few condensable species

File:

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

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 
+            !Fe_vap 
+            if (txt .eq. "Fe_vap") then 
                   q(:,iq)=0. !kg/kg of atmosphere
                   qsurf(iq) = 0. !kg/kg of atmosphere
@@ -697 +697 @@
                   endif
                   close(91)
-            endif ! (txt .eq. "generic_Fe_vap") 
-            !generic Fe_ice 
-            if (txt .eq. "generic_Fe_ice") then 
+            endif ! (txt .eq. "Fe_vap") 
+            !Fe_ice 
+            if (txt .eq. "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
-            !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
+
+            ! Cr_vap 
+            if (txt .eq. "Cr_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Cr_vap" input file
                   open(91,file='profile_Cr_vap',status='old',
      &            form="formatted",iostat=ierr)
@@ -719 +720 @@
                   endif
                   close(91)
-            endif ! (txt .eq. "generic_Cr_vap") 
-            !generic Cr_ice 
-            if (txt .eq. "generic_Cr_ice") then 
+            endif ! (txt .eq. "Cr_vap") 
+            !Cr_ice 
+            if (txt .eq. "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
+
+            ! Mg_vap 
+            if (txt .eq. "Mg_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Mg_vap" input file
+                  open(91,file='profile_Mg_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_Mg_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "Mg_vap") 
+            !Mg_ice 
+            if (txt .eq. "Mg_ice") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+            endif ! we only initialize the solid at 0
+
+            ! Na_vap 
+            if (txt .eq. "Na_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Na_vap" input file
+                  open(91,file='profile_Na_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_Na_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. Na_vap") 
+            !Cr_ice 
+            if (txt .eq. "Na_ice") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+            endif ! we only initialize the solid at 0
+            ! KCl_vap 
+            if (txt .eq. "KCl_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_KCL_vap" input file
+                  open(91,file='profile_KCl_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_KCl_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "KCl_vap") 
+            !KCl_ice 
+            if (txt .eq. "KCl_ice") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+            endif ! we only initialize the solid at 0
+
+            ! Mn_vap 
+            if (txt .eq. "Mn_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Mn_vap" input file
+                  open(91,file='profile_Mn_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_Mn_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "Mn_vap") 
+            !Mn_ice 
+            if (txt .eq. "Mn_ice") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+            endif ! we only initialize the solid at 0
+
+            ! Zn_vap 
+            if (txt .eq. "Zn_vap") then 
+                  q(:,iq)=0. !kg/kg of atmosphere
+                  qsurf(iq) = 0. !kg/kg of atmosphere
+                  ! Look for a "profile_Zn_vap" input file
+                  open(91,file='profile_Zn_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_Zn_vap file !"
+                  endif
+                  close(91)
+            endif ! (txt .eq. "Zn_vap") 
+            !Zn_ice 
+            if (txt .eq. "Zn_ice") then 
                   q(:,iq)=0. !kg/kg of atmosphere
                   qsurf(iq) = 0. !kg/kg of atmosphere