Changeset 4444 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_save.F90

Timestamp:

Feb 21, 2023, 3:26:41 PM (20 months ago)

Author:

idelkadi

Message:

Update of the ECRAD radiative code version implemented in the LMDZ model.
Upgrade to the :
​https://github.com/ecmwf/ecrad/trunk
Version svn : 749
UUID du dépôt : 44b0ca93-0ed8-356e-d663-ce57b7db7bff

File:

• LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_save.F90 (modified) (5 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_save.F90

@@ -33 +33 @@
   ! thermodynamics object
   subroutine save_fluxes(file_name, config, thermodynamics, flux, &
-       &                 iverbose, is_hdf5_file, experiment_name)
+       &                 iverbose, is_hdf5_file, experiment_name, &
+       &                 is_double_precision)
 
     use yomhook,                  only : lhook, dr_hook
@@ -50 +51 @@
     integer,          optional, intent(in) :: iverbose
     logical,          optional, intent(in) :: is_hdf5_file
+    logical,          optional, intent(in) :: is_double_precision
     character(len=*), optional, intent(in) :: experiment_name
 
@@ -96 +98 @@
     ! output file column varies most slowly so need to transpose
     call out_file%transpose_matrices(.true.)
+
+    ! Set default precision for file, if specified
+    if (present(is_double_precision)) then
+      call out_file%double_precision(is_double_precision)
+    end if
 
     ! Spectral fluxes in memory are dimensioned (nband,ncol,nlev), but
@@ -885 +892 @@
            &   dim2_name="column", dim1_name="level", &
            &   units_str="m", long_name="Ice effective radius")
-      if (allocated(cloud%re_ice)) then
+      if (associated(cloud%re_ice)) then
         call out_file%define_variable("re_ice", &
              &   dim2_name="column", dim1_name="level", &
@@ -966 +973 @@
       call out_file%put("q_ice", cloud%q_ice)
       call out_file%put("re_liquid", cloud%re_liq)
-      if (allocated(cloud%re_ice)) then
+      if (associated(cloud%re_ice)) then
         call out_file%put("re_ice", cloud%re_ice)
       end if