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

Timestamp:

Oct 1, 2018 2:37:10 AM (6 years ago)

Author:

raasch

Message:

modularization of the ocean code

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! changes concerning modularization of ocean option
+! 
+! 3289 2018-09-28 10:23:58Z suehring
 ! Introduce module parameter for number of inflow profiles
 ! 
@@ -549 +552 @@
     USE netcdf_data_input_mod,                                                 &
         ONLY:  init_3d, netcdf_data_input_init_3d
+
+    USE ocean_mod,                                                             &
+        ONLY:  ocean_init, ocean_init_arrays
     
     USE particle_attributes,                                                   &
@@ -772 +778 @@
                  s_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
                  s_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-#endif
-    ENDIF
-
-    IF ( ocean )  THEN
-#if defined( __nopointer )
-       ALLOCATE( prho(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
-                 rho_ocean(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
-                 sa(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                            &
-                 sa_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
-                 tsa_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-#else
-       ALLOCATE( prho_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                        &
-                 rho_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                         &
-                 sa_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
-                 sa_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
-                 sa_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-       prho => prho_1
-       rho_ocean  => rho_1  ! routines calc_mean_profile and diffusion_e require
-                      ! density to be apointer
 #endif
     ENDIF
@@ -1010 +997 @@
        s => s_1;  s_p => s_2;  ts_m => s_3
     ENDIF    
-
-    IF ( ocean )  THEN
-       sa => sa_1;  sa_p => sa_2;  tsa_m => sa_3
-    ENDIF
 #endif
+
 !
 !-- Initialize arrays for turbulence closure
@@ -1022 +1006 @@
     CALL init_surface_arrays
 !
-!-- Allocate microphysics module arrays
-    IF ( bulk_cloud_model )  THEN
-       CALL bcm_init_arrays
-    ENDIF
-!
-!-- Allocate land surface model arrays
-    IF ( land_surface )  THEN
-       CALL lsm_init_arrays
-    ENDIF
-
-!
-!-- Allocate wind turbine model arrays
-    IF ( wind_turbine )  THEN
-       CALL wtm_init_arrays
-    ENDIF
-!
-!-- Allocate gust module arrays
-    IF ( gust_module_enabled )  THEN
-       CALL gust_init_arrays
-    ENDIF
+!-- Allocate arrays for other modules
+    IF ( bulk_cloud_model    )  CALL bcm_init_arrays
+    IF ( gust_module_enabled )  CALL gust_init_arrays
+    IF ( land_surface        )  CALL lsm_init_arrays
+    IF ( ocean_mode          )  CALL ocean_init_arrays
+    IF ( wind_turbine        )  CALL wtm_init_arrays
+    IF ( uv_exposure         )  CALL uvem_init_arrays
 
 !
@@ -1053 +1024 @@
     IF ( uv_exposure )  THEN
        CALL uvem_init
-    ENDIF
-!
-!-- Allocate uv exposure arrays
-    IF ( uv_exposure )  THEN
-       CALL uvem_init_arrays
     ENDIF
 
@@ -1230 +1196 @@
           ENDIF
 
-          IF ( ocean )  THEN
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   sa(:,j,i) = sa_init
-                ENDDO
-             ENDDO
-          ENDIF
 !
 !--       Set velocity components at non-atmospheric / oceanic grid points to 
@@ -1380 +1339 @@
           ENDIF
 
-          IF ( ocean )  THEN
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   sa(:,j,i) = sa_init
-                ENDDO
-             ENDDO
-          ENDIF
 !
 !--       Compute initial temperature field and other constants used in case
@@ -1449 +1401 @@
        hom(:,1,7,:)  = SPREAD( pt(:,nys,nxl), 2, statistic_regions+1 )
 
-
-!
-!--    Store initial salinity profile
-       IF ( ocean )  THEN
-          hom(:,1,26,:)  = SPREAD( sa(:,nys,nxl), 2, statistic_regions+1 )
-       ENDIF
-
        IF ( humidity )  THEN
 !
@@ -1557 +1502 @@
        ENDIF       
 
-       IF ( ocean )  THEN
-          tsa_m = 0.0_wp
-          sa_p  = sa
-       ENDIF
-       
        CALL location_message( 'finished', .TRUE. )
 
@@ -1810 +1750 @@
        ENDIF
        IF ( passive_scalar )  s_p  = s
-       IF ( ocean          )  sa_p = sa
-
 !
 !--    Allthough tendency arrays are set in prognostic_equations, they have
@@ -1821 +1759 @@
        ENDIF
        IF ( passive_scalar )  ts_m  = 0.0_wp
-       IF ( ocean          )  tsa_m = 0.0_wp
 !
 !--    Initialize synthetic turbulence generator in case of restart.
@@ -2241 +2178 @@
     IF ( dt_dvrp /= 9999999.9_wp )  CALL init_dvrp
 
-    IF ( ocean )  THEN
-!
-!--    Initialize quantities needed for the ocean model
-       CALL init_ocean
-
-    ENDIF
-!
-!--    Initialize quantities for handling cloud physics
-!--    This routine must be called before lpm_init, because
-!--    otherwise, array d_exner, needed in data_output_dvrp (called by
-!--    lpm_init) is not defined.
-    IF ( .NOT. ocean )  THEN
+!
+!-- Initialize quantities needed for the ocean model
+    IF ( ocean_mode )  CALL ocean_init
+
+!
+!-- Initialize quantities for handling cloud physics.
+!-- This routine must be called before lpm_init, becaus otherwise,
+!-- array d_exner, needed in data_output_dvrp (called by lpm_init) is not defined.
+    IF ( .NOT. ocean_mode )  THEN
 
        ALLOCATE( hyp(nzb:nzt+1) )
@@ -2454 +2388 @@
     rdf_sc = 0.0_wp
     IF ( rayleigh_damping_factor /= 0.0_wp )  THEN
-       IF (  .NOT.  ocean )  THEN
+
+       IF (  .NOT.  ocean_mode )  THEN
           DO  k = nzb+1, nzt
              IF ( zu(k) >= rayleigh_damping_height )  THEN
@@ -2464 +2399 @@
           ENDDO
        ELSE
+!
+!--       In ocean mode, rayleigh damping is applied in the lower part of the
+!--       model domain
           DO  k = nzt, nzb+1, -1
              IF ( zu(k) <= rayleigh_damping_height )  THEN
@@ -2473 +2411 @@
           ENDDO
        ENDIF
+
     ENDIF
     IF ( scalar_rayleigh_damping )  rdf_sc = rdf