Changeset 3310 for trunk/LMDZ.GENERIC
- Timestamp:
- Apr 22, 2024, 4:34:50 PM (8 months ago)
- Location:
- trunk/LMDZ.GENERIC/libf/phystd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90
r2870 r3310 133 133 if (is_master) print*,'Initialising H2-H2 continuum from HITRAN database...' 134 134 135 ! 1.1 Open the ASCII files and set nS according to version135 ! 1.1 Open the ASCII files and set nS according to version 136 136 ! Only two possible versions for now : 2011 or 2018 (sanity check in inifis_mod) 137 137 if (versH2H2cia.eq.2011) then 138 nS = 2428139 if (H2orthopara_mixture.eq."normal") then140 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2011.cia'141 else if (H2orthopara_mixture.eq."equilibrium") then142 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2011.cia'143 else if (H2orthopara_mixture.eq."hot") then144 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_2011_extended.cia'145 ns = 800146 endif147 else if (versH2H2cia.eq.2018 ) then148 nS = 9600149 if (H2orthopara_mixture.eq."normal") then150 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2018.cia'151 else if (H2orthopara_mixture.eq."equilibrium") then152 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2018.cia'153 154 endif138 nS = 2428 139 if (H2orthopara_mixture.eq."normal") then 140 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2011.cia' 141 else if (H2orthopara_mixture.eq."equilibrium") then 142 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2011.cia' 143 else if (H2orthopara_mixture.eq."hot") then 144 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_2011_extended.cia' 145 ns = 800 146 endif 147 else if (versH2H2cia.eq.2018 .and. H2orthopara_mixture.eq."normal") then 148 nS = 9600 149 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2018.cia' 150 else if (versH2H2cia.eq.2018 .and. H2orthopara_mixture.eq."equilibrium") then 151 nS = 10860 152 dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_0-15000cm-1_40-400K.cia' 153 endif 154 155 155 156 156 if(.not.allocated(wn_arr)) allocate(wn_arr(nS)) -
trunk/LMDZ.GENERIC/libf/phystd/interpolateH2He.F90
r2870 r3310 46 46 logical firstcall 47 47 48 save wn_arr, temp_arr, abs_arr !read by master48 save ns, wn_arr, temp_arr, abs_arr !read by master 49 49 50 50 character*100 dt_file … … 64 64 if (H2orthopara_mixture .eq. "hot") then ! .and. (temp .gt. 3000.0)) then 65 65 ! print*,"We're in the Hot Jupiter case " 66 nS = 1998167 66 if (ztemp .gt. 9900) then 68 67 if (strictboundcia) then … … 86 85 endif ! temp .gt. 9900 87 86 else !not in Hot Jupiter case 88 nS = 242889 87 if(ztemp.gt.400)then 90 88 if (strictboundcia) then … … 128 126 endif ! of (H2orthopara_mixture .eq. "hot") 129 127 130 if (.not. allocated(wn_arr)) allocate(wn_arr(nS))131 if (.not. allocated(abs_arr)) allocate(abs_arr(nS,nT))132 128 133 129 amagatH2 = (273.15/temp)*(presH2/101325.0) … … 140 136 ! 1.1 Open the ASCII files 141 137 if (H2orthopara_mixture.eq."normal") then 142 dt_file=TRIM(datadir)//'/continuum_data/H2-He_norm_2011.cia' 138 nS = 2428 139 dt_file=TRIM(datadir)//'/continuum_data/H2-He_norm_2011.cia' 143 140 else if (H2orthopara_mixture.eq."equilibrium") then 144 dt_file=TRIM(datadir)//'/continuum_data/H2-He_eq_2011.cia' 141 nS = 1501 142 dt_file=TRIM(datadir)//'/continuum_data/H2-He_eq_0-15000cm-1_40-400K.cia' 145 143 else if (H2orthopara_mixture .eq. "hot") then !we use a dataset than can go as high as 9900 K 146 dt_file=TRIM(datadir)//'/continuum_data/H2-He_2011.cia' 144 nS = 19981 145 dt_file=TRIM(datadir)//'/continuum_data/H2-He_2011.cia' 147 146 endif 147 148 if (.not. allocated(wn_arr)) allocate(wn_arr(nS)) 149 if (.not. allocated(abs_arr)) allocate(abs_arr(nS,nT)) 150 148 151 149 152 !$OMP MASTER … … 190 193 print*,' and He partial pressure ',presHe,' Pa' 191 194 endif 192 endif 195 endif !! if (firstcall) 193 196 194 197 call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,ztemp,abcoef,ind)
Note: See TracChangeset
for help on using the changeset viewer.