Changeset 4226 for palm/trunk/SOURCE/radiation_model_mod.f90

Timestamp:

Sep 10, 2019 5:03:24 PM (5 years ago)

Author:

suehring

Message:

Offline nesting: data input modularized; time variable is defined relative to time_utc_init, so that input data is correctly mapped to actual model time; checks rephrased and new checks for the time dimension added; Netcdf input: routine to retrieve dimension length renamed

File:

• palm/trunk/SOURCE/radiation_model_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! - Netcdf input routine for dimension length renamed
+! - Define time variable for external radiation input relative to time_utc_init 
+! 
+! 4210 2019-09-02 13:07:09Z suehring
 ! - Revise steering of splitting diffuse and direct radiation
 ! - Bugfixes in checks
@@ -34 +38 @@
 ! 
 ! 4208 2019-09-02 09:01:07Z suehring
-! Bugfix in accessing albedo_pars in the clear-sky branch (merge from branch) 
+! Bugfix in accessing albedo_pars in the clear-sky branch
+! (merge from branch resler) 
 ! 
 ! 4198 2019-08-29 15:17:48Z gronemeier
@@ -259 +264 @@
                close_input_file,                                               &
                get_attribute,                                                  &
+               get_dimension_length,                                           &
                get_variable,                                                   &
                inquire_num_variables,                                          &
@@ -264 +270 @@
                input_file_dynamic,                                             &
                input_pids_dynamic,                                             &
-               netcdf_data_input_get_dimension_length,                         &
                num_var_pids,                                                   &
                pids_id,                                                        &
@@ -2510 +2515 @@
 !--       Input time dimension.
           IF ( check_existence( vars_pids, 'time_rad' ) )  THEN
-             CALL netcdf_data_input_get_dimension_length( pids_id,             &
-                                                          ntime,               &
-                                                          'time_rad' )
+             CALL get_dimension_length( pids_id, ntime, 'time_rad' )
           
              ALLOCATE( time_rad_f%var1d(0:ntime-1) )
@@ -2637 +2640 @@
           ENDIF
           
-          IF ( time_rad_f%var1d(0) /= 0.0_wp )  THEN
+          IF ( time_rad_f%var1d(0) /= time_utc_init )  THEN
              message_string = 'External radiation forcing: first point in ' // &
-                              'time is /= 0.0.'
+                              'time is /= time_utc_init.'
              CALL message( 'radiation_init', 'PA0313', 1, 2, 0, 6, 0 )
           ENDIF
           
-          IF ( end_time - spinup_time > time_rad_f%var1d(ntime-1) )  THEN
+          IF ( end_time - spinup_time > time_rad_f%var1d(ntime-1)              &
+                                      - time_utc_init )  THEN
              message_string = 'External radiation forcing does not cover ' //  &
                               'the entire simulation time.'
@@ -2832 +2836 @@
        ELSE
           t = 0
-          DO WHILE ( time_rad_f%var1d(t) <= time_since_reference_point )
+          DO WHILE ( time_rad_f%var1d(t) <=                                    &
+                     time_since_reference_point + time_utc_init )
              t = t + 1
           ENDDO
@@ -2838 +2843 @@
           tm = MAX( t-1, 0 )
           
-          fac_dt = ( time_since_reference_point - time_rad_f%var1d(tm) + dt_3d ) &
+          fac_dt = ( time_since_reference_point + time_utc_init                &
+                   - time_rad_f%var1d(tm) + dt_3d )                            &
                  / ( time_rad_f%var1d(t)  - time_rad_f%var1d(tm) )
           fac_dt = MIN( 1.0_wp, fac_dt )