Changeset 3700 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Jan 26, 2019 5:03:42 PM (5 years ago)

Author:

knoop

Message:

Moved user_define_netdf_grid into user_module.f90
Added module interface for the definition of additional timeseries

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -617 +617 @@
 
     USE netcdf_interface,                                                      &
-        ONLY:  dots_max, dots_num, dots_unit, dots_label
+        ONLY:  dots_max
 
     USE netcdf_data_input_mod,                                                 &
@@ -734 +734 @@
               sums_l(nzb:nzt+1,pr_palm+max_pr_user+max_pr_cs,0:threads_per_task-1),      &
               sums_l_l(nzb:nzt+1,0:statistic_regions,0:threads_per_task-1),    &
-              sums_wsts_bc_l(nzb:nzt+1,0:statistic_regions),                   &
-              ts_value(dots_max,0:statistic_regions) )
+              sums_wsts_bc_l(nzb:nzt+1,0:statistic_regions) )
+    ALLOCATE( ts_value(dots_max,0:statistic_regions) )
     ALLOCATE( ptdf_x(nxlg:nxrg), ptdf_y(nysg:nyng) )
 
@@ -2041 +2041 @@
     rmask(:,nxlg:nxl-1,:) = 0.0_wp;  rmask(:,nxr+1:nxrg,:) = 0.0_wp
     rmask(nysg:nys-1,:,:) = 0.0_wp;  rmask(nyn+1:nyng,:,:) = 0.0_wp
-
-!
-!-- Temporary solution to add LSM and radiation time series to the default
-!-- output
-    IF ( land_surface  .OR.  radiation )  THEN
-       IF ( TRIM( radiation_scheme ) == 'rrtmg' )  THEN
-          dots_num = dots_num + 15
-       ELSE
-          dots_num = dots_num + 11
-       ENDIF
-    ENDIF
 !
 !-- To do: New concept for these non-topography grid points!
@@ -2382 +2371 @@
        ENDDO
     ENDIF
-!
-!-- Check if maximum number of allowed timeseries is exceeded
-    IF ( dots_num > dots_max )  THEN
-       WRITE( message_string, * ) 'number of time series quantities exceeds',  &
-                                  ' its maximum of dots_max = ', dots_max,     &
-                                  '&Please increase dots_max in modules.f90.'
-       CALL message( 'init_3d_model', 'PA0194', 1, 2, 0, 6, 0 )    
-    ENDIF
 
 !