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

Timestamp:

May 30, 2018 5:43:55 PM (6 years ago)

Author:

suehring

Message:

Speed-up NetCDF input; Revise NetCDF-input routines and remove input via io-blocks; Temporarily revoke renaming of input variables in dynamic driver; More detailed error messages created; Bugfix in mapping 3D buildings; Bugfix in land-surface model at pavement surfaces; Bugfix in initialization with inifor

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (7 diffs, 1 prop)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Move initialization call for nudging and 1D/3D offline nesting. 
+! Revise initialization with inifor data.
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Error messages revised
 ! 
@@ -1116 +1120 @@
 
 !
-!-- Initialize nudging if required
-    IF ( nudging )  THEN
-       CALL nudge_init
-    ENDIF
-
-!
-!-- Initialize reading of large scale forcing from external file - if required
-    IF ( large_scale_forcing  .OR.  forcing )  THEN
-       CALL lsf_init
-    ENDIF
-
-!
 !-- Allocate arrays containing the RK coefficient for calculation of 
 !-- perturbation pressure and turbulent fluxes. At this point values are
@@ -1166 +1158 @@
           CALL location_message( 'initializing with INIFOR', .FALSE. )
 !
-!--       Read initial 1D profiles from NetCDF file if available. 
+!--       Read initial 1D profiles or 3D data from NetCDF file, depending 
+!--       on the provided level-of-detail. 
 !--       At the moment, only u, v, w, pt and q are provided. 
           CALL netcdf_data_input_init_3d
@@ -1177 +1170 @@
           IF ( ANY( zu(1:nzt+1) /= init_3d%zu_atmos(1:init_3d%nzu) ) )  THEN
 
-             CALL netcdf_data_input_interpolate( init_3d%u_init(nzb+1:nzt+1),  &
+             IF( init_3d%lod_u == 1 )                                          &
+                CALL netcdf_data_input_interpolate(                            &
+                                                 init_3d%u_init(nzb+1:nzt+1),  &
                                                  zu(nzb+1:nzt+1),              &
                                                  init_3d%zu_atmos )
-             CALL netcdf_data_input_interpolate( init_3d%v_init(nzb+1:nzt+1),  &
+                                                 
+             IF( init_3d%lod_v == 1 )                                          &
+                CALL netcdf_data_input_interpolate(                            &
+                                                 init_3d%v_init(nzb+1:nzt+1),  &
                                                  zu(nzb+1:nzt+1),              &
                                                  init_3d%zu_atmos )
+                                                 
 !              CALL netcdf_data_input_interpolate( init_3d%w_init(nzb+1:nzt),    &
 !                                                  zw(nzb+1:nzt),                &
 !                                                  init_3d%zw_atmos )
-             IF ( .NOT. neutral )                                              &
+
+             IF ( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                  &
                 CALL netcdf_data_input_interpolate(                            &
                                              init_3d%pt_init(nzb+1:nzt+1),     &
                                              zu(nzb+1:nzt+1),                  &
                                              init_3d%zu_atmos )
-             IF ( humidity )                                                   &
+                                             
+             IF ( humidity  .AND.  init_3d%lod_q == 1 )                        &
                 CALL netcdf_data_input_interpolate(                            &
                                              init_3d%q_init(nzb+1:nzt+1),      &
@@ -1197 +1198 @@
                                              init_3d%zu_atmos )
           ENDIF
-
-          u_init = init_3d%u_init
-          v_init = init_3d%v_init    
-          IF( .NOT. neutral )  pt_init = init_3d%pt_init
-          IF( humidity      )  q_init  = init_3d%q_init
+!
+!--       In case of LOD=1, initialize 1D profiles and 3D data.  
+          IF( init_3d%lod_u == 1 )  u_init = init_3d%u_init
+          IF( init_3d%lod_v == 1 )  v_init = init_3d%v_init    
+          IF( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                      &
+             pt_init = init_3d%pt_init
+          IF( humidity  .AND.  init_3d%lod_q == 1 )                            &
+             q_init  = init_3d%q_init
 
 !
@@ -1210 +1214 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                u(:,j,i) = u_init(:)
-                v(:,j,i) = v_init(:)
-                IF( .NOT. neutral )  pt(:,j,i) = pt_init(:)
-                IF( humidity      )  q(:,j,i)  = q_init(:)
+                IF( init_3d%lod_u == 1 )  u(:,j,i) = u_init(:)
+                IF( init_3d%lod_v == 1 )  v(:,j,i) = v_init(:)
+                IF( .NOT. neutral  .AND.  init_3d%lod_pt == 1 )                &
+                   pt(:,j,i) = pt_init(:)
+                IF( humidity  .AND.  init_3d%lod_q == 1 )                      &
+                   q(:,j,i)  = q_init(:)
              ENDDO
           ENDDO
@@ -2287 +2293 @@
 
 !
+!-- Initialize nudging if required
+    IF ( nudging )  CALL nudge_init
+
+!
+!-- Initialize 1D/3D offline-nesting with COSMO model and read data from 
+!-- external file.
+    IF ( large_scale_forcing  .OR.  forcing )  CALL lsf_init
+    
+!
 !-- Initialize surface forcing corresponding to large-scale forcing. Therein, 
 !-- initialize heat-fluxes, etc. via datatype. Revise it later!