Changeset 3864 for palm/trunk/SOURCE/netcdf_data_input_mod.f90

Timestamp:

Apr 5, 2019 9:01:56 AM (5 years ago)

Author:

monakurppa

Message:

major changes in salsa: data input, format and performance

• Time-dependent emissions enabled: lod=1 for yearly PM emissions that are normalised depending on the time, and lod=2 for preprocessed emissions (similar to the chemistry module).
• Additionally, 'uniform' emissions allowed. This emission is set constant on all horisontal upward facing surfaces and it is created based on parameters surface_aerosol_flux, aerosol_flux_dpg/sigmag/mass_fracs_a/mass_fracs_b.
• All emissions are now implemented as surface fluxes! No 3D sources anymore.
• Update the emission information by calling salsa_emission_update if skip_time_do_salsa >= time_since_reference_point and next_aero_emission_update <= time_since_reference_point
• Aerosol background concentrations read from PIDS_DYNAMIC. The vertical grid must match the one applied in the model.
• Gas emissions and background concentrations can be also read in in salsa_mod if the chemistry module is not applied.
• In deposition, information on the land use type can be now imported from the land use model
• Use SI units in PARIN, i.e. n_lognorm given in #/m3 and dpg in metres.
• Apply 100 character line limit
• Change all variable names from capital to lowercase letter
• Change real exponents to integer if possible. If not, precalculate the value of exponent
• Rename in1a to start_subrange_1a, fn2a to end_subrange_1a etc.
• Rename nbins --> nbins_aerosol, ncc_tot --> ncomponents_mass and ngast --> ngases_salsa
• Rename ibc to index_bc, idu to index_du etc.
• Renamed loop indices b, c and sg to ib, ic and ig
• run_salsa subroutine removed
• Corrected a bud in salsa_driver: falsely applied ino instead of inh
• Call salsa_tendency within salsa_prognostic_equations which is called in module_interface_mod instead of prognostic_equations_mod
• Removed tailing white spaces and unused variables
• Change error message to start by PA instead of SA

File:

• palm/trunk/SOURCE/netcdf_data_input_mod.f90 (modified) (19 diffs)

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! get_variable_4d_to_3d_real modified to enable read in data of type
+! data(t,y,x,n) one timestep at a time + some routines made public
+! 
+! 3855 2019-04-03 10:00:59Z suehring
 ! Typo removed
 ! 
@@ -913 +917 @@
            netcdf_data_input_interpolate, netcdf_data_input_offline_nesting,   &
            netcdf_data_input_surface_data, netcdf_data_input_topo,             &
-           netcdf_data_input_var, get_attribute, get_variable, open_read_file
+           netcdf_data_input_var, get_attribute, get_variable, open_read_file, &
+           check_existence, inquire_num_variables, inquire_variable_names
 
 
@@ -5557 +5562 @@
        INTEGER(iwp)                  ::  ke              !< start index of 3rd dimension
        INTEGER(iwp)                  ::  ks              !< end index of 3rd dimension
-       
+
        INTEGER(KIND=1), DIMENSION(:,:,:), ALLOCATABLE   ::  tmp  !< temporary variable to read data from file according 
                                                                  !< to its reverse memory access
@@ -5566 +5571 @@
 !
 !--    Inquire variable id
-       nc_stat = NF90_INQ_VARID( id, TRIM( variable_name ), id_var )   
+       nc_stat = NF90_INQ_VARID( id, TRIM( variable_name ), id_var )
 !
 !--    Check for collective read-operation and set respective NetCDF flags if
@@ -5573 +5578 @@
 #if defined( __netcdf4_parallel )
           nc_stat = NF90_VAR_PAR_ACCESS (id, id_var, NF90_COLLECTIVE)
-#endif          
+#endif
        ENDIF
 !
@@ -5582 +5587 @@
        nc_stat = NF90_GET_VAR( id, id_var, tmp,                                &
                                start = (/ is+1,    js+1,    ks+1 /),           &
-                               count = (/ ie-is+1, je-js+1, ke-ks+1 /) )                              
-
-       CALL handle_error( 'get_variable_3d_int8', 533, variable_name )                               
+                               count = (/ ie-is+1, je-js+1, ke-ks+1 /) )
+
+       CALL handle_error( 'get_variable_3d_int8', 533, variable_name )
 !
 !--    Resort data. Please note, dimension subscripts of var all start at 1. 
@@ -5594 +5599 @@
           ENDDO
        ENDDO
-       
+
        DEALLOCATE( tmp )
 
@@ -5627 +5632 @@
        INTEGER(iwp)                  ::  ke              !< start index of 3rd dimension
        INTEGER(iwp)                  ::  ks              !< end index of 3rd dimension
-       
+
        REAL(wp), DIMENSION(:,:,:), ALLOCATABLE   ::  tmp !< temporary variable to read data from file according 
                                                          !< to its reverse memory access
@@ -5641 +5646 @@
 !--    required. 
        IF ( collective_read )  THEN
-#if defined( __netcdf4_parallel )       
+#if defined( __netcdf4_parallel )
           nc_stat = NF90_VAR_PAR_ACCESS (id, id_var, NF90_COLLECTIVE)
 #endif
@@ -5653 +5658 @@
                                start = (/ is+1,    js+1,    ks+1 /),           &
                                count = (/ ie-is+1, je-js+1, ke-ks+1 /) )   
-                               
+
        CALL handle_error( 'get_variable_3d_real', 534, variable_name )
 !
@@ -5664 +5669 @@
           ENDDO
        ENDDO
-       
+
        DEALLOCATE( tmp )
 
@@ -5737 +5742 @@
        ELSE
 !
-!--    Allocate temporary variable according to memory access on file. 
-       ALLOCATE( tmp(is:ie,js:je,k1s:k1e,k2s:k2e) )
-!
-!--    Get variable
-          nc_stat = NF90_GET_VAR( id, id_var, tmp,                                &
+!--       Allocate temporary variable according to memory access on file. 
+          ALLOCATE( tmp(is:ie,js:je,k1s:k1e,k2s:k2e) )
+!
+!--       Get variable
+          nc_stat = NF90_GET_VAR( id, id_var, tmp,                             &
                                start = (/ is+1,    js+1,   k1s+1, k2s+1 /),    &
                                count = (/ ie-is+1, je-js+1,                    &
@@ -5748 +5753 @@
           CALL handle_error( 'get_variable_4d_real', 535, variable_name )
 !
-!--    Resort data. Please note, dimension subscripts of var all start at 1. 
+!--       Resort data. Please note, dimension subscripts of var all start at 1. 
           DO  i = is, ie 
              DO  j = js, je
@@ -5758 +5763 @@
              ENDDO
           ENDDO
-       
+
           DEALLOCATE( tmp )
        ENDIF
@@ -5828 +5833 @@
        ELSE
 !
-!--    Allocate temporary variable according to memory access on file. 
+!--       Allocate temporary variable according to memory access on file. 
           ALLOCATE( tmp(is:ie,js:je,ks:ke) )
 !
-!--    Get variable
-          nc_stat = NF90_GET_VAR( id, id_var, tmp,                                &
-                                  start = (/ ns+1, is+1,    js+1,    ks+1 /),           &
-                                  count = (/ 1, ie-is+1, je-js+1, ke-ks+1 /) )   
-                               
+!--       Get variable
+          nc_stat = NF90_GET_VAR( id, id_var, tmp,                             &
+                                  start = (/ is+1,    js+1,    ks+1,   ns+1 /),&
+                                  count = (/ ie-is+1, je-js+1, ke-ks+1, 1   /) )
+
           CALL handle_error( 'get_variable_4d_to_3d_real', 536, variable_name )
 !
-!--    Resort data. Please note, dimension subscripts of var all start at 1. 
-          DO  i = is, ie 
+!--       Resort data. Please note, dimension subscripts of var all start at 1.
+          DO  i = is, ie
              DO  j = js, je
                 DO  k = ks, ke
@@ -5846 +5851 @@
              ENDDO
           ENDDO
-       
+
          DEALLOCATE( tmp )
 
@@ -5962 +5967 @@
        CHARACTER(LEN=*)              ::  variable_name   !< variable name
 
-       INTEGER(iwp)                  ::  ns              !< start index for subdomain input along n dimension: ns coincides here with ne, since, we select only one value along the 1st dimension n
+       INTEGER(iwp)                  ::  ns              !< start index for subdomain input along n dimension: 
+                                                         !< ns coincides here with ne, since, we select only one 
+                                                         !< value along the 1st dimension n
 
        INTEGER(iwp)                  ::  t               !< index along t direction
@@ -5986 +5993 @@
 !
 !--    Inquire variable id
-       nc_stat = NF90_INQ_VARID( id, TRIM( variable_name ), id_var )  
+       nc_stat = NF90_INQ_VARID( id, TRIM( variable_name ), id_var )
 !
 !--    Check for collective read-operation and set respective NetCDF flags if
-!--    required. 
+!--    required.
        IF ( collective_read )  THEN
           nc_stat = NF90_VAR_PAR_ACCESS (id, id_var, NF90_COLLECTIVE)
@@ -6036 +6043 @@
           ENDDO 
 
-         DEALLOCATE( tmp )
+          DEALLOCATE( tmp )
 
        ENDIF