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

Timestamp:

Jul 27, 2007 9:09:17 AM (17 years ago)

Author:

raasch

Message:

preliminary version for coupled runs

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Flux initialization in case of coupled runs, +momentum fluxes at top boundary
 !
 ! Former revisions:
@@ -102 +102 @@
               ts(nys-1:nyn+1,nxl-1:nxr+1), tswst_1(nys-1:nyn+1,nxl-1:nxr+1),  &
               us(nys-1:nyn+1,nxl-1:nxr+1), usws_1(nys-1:nyn+1,nxl-1:nxr+1),   &
-              vsws_1(nys-1:nyn+1,nxl-1:nxr+1), z0(nys-1:nyn+1,nxl-1:nxr+1) )
+              uswst_1(nys-1:nyn+1,nxl-1:nxr+1),                               &
+              vsws_1(nys-1:nyn+1,nxl-1:nxr+1),                                &
+              vswst_1(nys-1:nyn+1,nxl-1:nxr+1), z0(nys-1:nyn+1,nxl-1:nxr+1) )
 
     IF ( timestep_scheme(1:5) /= 'runge' )  THEN
@@ -111 +113 @@
                  tswst_2(nys-1:nyn+1,nxl-1:nxr+1), &
                  usws_2(nys-1:nyn+1,nxl-1:nxr+1),  &
+                 uswst_2(nys-1:nyn+1,nxl-1:nxr+1), &
+                 vswst_2(nys-1:nyn+1,nxl-1:nxr+1), &
                  vsws_2(nys-1:nyn+1,nxl-1:nxr+1) )
     ENDIF
@@ -260 +264 @@
        tswst_m => tswst_1;  tswst => tswst_2
        usws_m  => usws_1;   usws  => usws_2
+       uswst_m => uswst_1;  uswst => uswst_2
        vsws_m  => vsws_1;   vsws  => vsws_2
+       vswst_m => vswst_1;  vswst => vswst_2
        e_m  => e_1;   e  => e_2;   e_p  => e_3;   te_m  => e_3
        kh_m => kh_1;  kh => kh_2
@@ -287 +293 @@
        tswst => tswst_1
        usws  => usws_1
+       uswst => uswst_1
        vsws  => vsws_1
+       vswst => vswst_1
        e     => e_1;   e_p  => e_2;   te_m  => e_3;   e_m  => e_3
        kh    => kh_1
@@ -381 +389 @@
              vsws = 0.0
           ENDIF
+          uswst = top_momentumflux_u
+          vswst = top_momentumflux_v
 
 !
@@ -494 +504 @@
                            ! e = (u*/cm)**2 at k=nzb+1
           ENDIF
-          e    = 0.0
-          rif  = 0.0
-          ts   = 0.0
-          us   = 0.0
-          usws = 0.0
-          vsws = 0.0
+          e     = 0.0
+          rif   = 0.0
+          ts    = 0.0
+          us    = 0.0
+          usws  = 0.0
+          uswst = top_momentumflux_u
+          vsws  = 0.0
+          vswst = top_momentumflux_v
           IF ( humidity  .OR.  passive_scalar ) qs = 0.0
 
@@ -617 +629 @@
                 saswst = top_salinityflux
              ENDIF
-         ENDIF
+          ENDIF
+
+!
+!--       Initialization in case of a coupled model run
+          IF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
+             tswst = 0.0
+             IF ( ASSOCIATED( tswst_m ) )  tswst_m = tswst
+          ENDIF
 
        ENDIF