Changeset 4853 for LMDZ6/trunk/libf/phylmd/ecrad/driver
- Timestamp:
- Mar 19, 2024, 3:34:21 PM (6 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd/ecrad/driver
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/ecrad/driver/Makefile
r4773 r4853 20 20 test_programs: $(TEST_PROGRAMS) 21 21 22 # Link ecrad executable; add "-lifs" if you want to use the "satur" 23 # routine in ecrad_driver.F90 22 24 $(EXECUTABLE): $(OBJECTS) ../lib/*.a ecrad_driver.o 23 25 $(FC) $(FCFLAGS) ecrad_driver.o $(OBJECTS) $(LIBS) -o $(EXECUTABLE) -
LMDZ6/trunk/libf/phylmd/ecrad/driver/ecrad_driver.F90
r4773 r4853 55 55 implicit none 56 56 57 ! Uncomment this if you want to use the "satur" routine below 58 !#include "satur.intfb.h" 59 57 60 ! The NetCDF file containing the input profiles 58 61 type(netcdf_file) :: file … … 269 272 270 273 ! Compute saturation with respect to liquid (needed for aerosol 271 ! hydration) call 274 ! hydration) call... 272 275 call thermodynamics%calc_saturation_wrt_liquid(driver_config%istartcol,driver_config%iendcol) 273 276 277 ! ...or alternatively use the "satur" function in the IFS (requires 278 ! adding -lifs to the linker command line) but note that this 279 ! computes saturation with respect to ice at colder temperatures, 280 ! which is almost certainly incorrect 281 !allocate(thermodynamics%h2o_sat_liq(ncol,nlev)) 282 !call satur(driver_config%istartcol, driver_config%iendcol, ncol, 1, nlev, .false., & 283 ! 0.5_jprb * (thermodynamics.pressure_hl(:,1:nlev)+thermodynamics.pressure_hl(:,2:nlev)), & 284 ! 0.5_jprb * (thermodynamics.temperature_hl(:,1:nlev)+thermodynamics.temperature_hl(:,2:nlev)), & 285 ! thermodynamics%h2o_sat_liq, 2) 286 274 287 ! Check inputs are within physical bounds, printing message if not 275 288 is_out_of_bounds = gas%out_of_physical_bounds(driver_config%istartcol, driver_config%iendcol, &
Note: See TracChangeset
for help on using the changeset viewer.