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

Timestamp:

Aug 12, 2016 3:14:59 PM (8 years ago)

Author:

suehring

Message:

Prescribing scalar flux at model top; several bugfixes concering data output of scalars and output of flight data

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Initializaton of scalarflux at model top
+! Bugfixes in initialization of surface and top salinity flux, top scalar and 
+! humidity fluxes
 ! 
 ! Former revisions:
@@ -1205 +1207 @@
                 ENDIF
              ENDIF
-          ENDIF          
+          ENDIF    
+!
+!--       Initialize near-surface salinity flux
+          IF ( ocean )  saswsb = bottom_salinityflux
 
        ENDIF
@@ -1214 +1219 @@
 !--    The latent flux is zero in this case!
        IF ( use_top_fluxes )  THEN
-
-          IF ( constant_top_heatflux )  THEN
-!
-!--          Heat flux is prescribed
-             tswst = top_heatflux
-
-             IF ( humidity )  THEN
-                qswst = 0.0_wp
-                IF ( cloud_physics  .AND.  microphysics_seifert ) THEN
-                   nrswst = 0.0_wp
-                   qrswst = 0.0_wp
-                ENDIF
+!
+!--       Prescribe to heat flux
+          IF ( constant_top_heatflux )  tswst = top_heatflux
+!
+!--       Prescribe zero latent flux at the top      
+          IF ( humidity )  THEN
+             qswst = 0.0_wp
+             IF ( cloud_physics  .AND.  microphysics_seifert ) THEN
+                nrswst = 0.0_wp
+                qrswst = 0.0_wp
              ENDIF
-             IF ( passive_scalar )  sswst = 0.0_wp
-
-             IF ( ocean )  THEN
-                saswsb = bottom_salinityflux
-                saswst = top_salinityflux
-             ENDIF
-          ENDIF
-
+          ENDIF
+!
+!--       Prescribe top scalar flux
+          IF ( passive_scalar .AND. constant_top_scalarflux )                  &
+             sswst = top_scalarflux
+!
+!--       Prescribe top salinity flux
+          IF ( ocean .AND. constant_top_salinityflux)                          &
+             saswst = top_salinityflux
 !
 !--       Initialization in case of a coupled model run