Changeset 1496 for palm/trunk/SOURCE/swap_timelevel.f90

Timestamp:

Dec 2, 2014 5:25:50 PM (9 years ago)

Author:

maronga

Message:

added beta version of a land surface model and a simple radiation model for clear sky conditions

File:

• palm/trunk/SOURCE/swap_timelevel.f90 (modified) (6 diffs)

palm/trunk/SOURCE/swap_timelevel.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added swapping of land surface model quantities
 ! 
 ! Former revisions:
@@ -69 +69 @@
         ONLY:  e, e_p, nr, nr_p, pt, pt_p, q, q_p, qr, qr_p, sa, sa_p, u, u_p, &
                v, v_p, w, w_p
+
+    USE land_surface_model_mod,                                                &
+        ONLY: land_surface, m_liq, m_liq_p, m_soil, m_soil_p, T_0, T_0_p,      &
+              T_soil, T_soil_p
 #else
     USE arrays_3d,                                                             &
@@ -74 +78 @@
                q_1, q_2, q_p, qr, qr_1, qr_2, qr_p, sa, sa_1, sa_2, sa_p, u,   &
                u_1, u_2, u_p, v, v_1, v_2, v_p, w, w_1, w_2, w_p
+
+    USE land_surface_model_mod,                                                &
+        ONLY: land_surface, m_liq, m_liq_1, m_liq_2, m_liq_p, m_soil,          &
+              m_soil_1, m_soil_2, m_soil_p, T_0, T_0_1, T_0_2, T_0_p, T_soil,  &
+              T_soil_1, T_soil_2, T_soil_p
 #endif
 
@@ -116 +125 @@
     ENDIF
 
+    IF ( land_surface )  THEN
+       T_0    = T_0_p
+       T_soil = T_soil_p
+       IF ( humidity )  THEN
+          m_soil = m_soil_p
+          m_liq  = m_liq_p
+       ENDIF
+    ENDIF
+
+
     CALL cpu_log( log_point(28), 'swap_timelevel (nop)', 'stop' )
 #else
@@ -145 +164 @@
           ENDIF
 
+          IF ( land_surface )  THEN
+             T_0    => T_0_1;    T_0_p    => T_0_2
+             T_soil => T_soil_1; T_soil_p => T_soil_2
+             IF ( humidity )  THEN
+                m_soil => m_soil_1; m_soil_p  => m_soil_2
+                m_liq  => m_liq_1;  m_liq_p   => m_liq_2
+             ENDIF
+          ENDIF
+
 
        CASE ( 1 )
@@ -169 +197 @@
           ENDIF
 
+          IF ( land_surface )  THEN
+             T_0    => T_0_2;    T_0_p    => T_0_1
+             T_soil => T_soil_2; T_soil_p => T_soil_1
+             IF ( humidity )  THEN
+                m_soil => m_soil_2; m_soil_p  => m_soil_1
+                m_liq  => m_liq_2;  m_liq_p   => m_liq_1
+             ENDIF
+          ENDIF
+
 
     END SELECT