Changeset 96 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Jun 4, 2007 8:07:41 AM (17 years ago)

Author:

raasch

Message:

more preliminary uncomplete changes for ocean version

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (12 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! Calculation of Deardorff velocity scale adjusted to be used with the ocean
+! Statistics for ocean version (salinity, density) added,
+! calculation of Deardorff velocity scale adjusted to be used with the ocean
 ! version (HAS STILL TO BE COMPLETED!!!)
 !
@@ -123 +124 @@
           ENDDO
        ENDDO
+
+!
+!--    Horizontally averaged profile of salinity
+       IF ( ocean )  THEN
+          !$OMP DO
+          DO  i = nxl, nxr
+             DO  j =  nys, nyn
+                DO  k = nzb_s_outer(j,i), nzt+1
+                   sums_l(k,23,tn)  = sums_l(k,23,tn) + &
+                                      sa(k,j,i) * rmask(j,i,sr)
+                ENDDO
+             ENDDO
+          ENDDO
+       ENDIF
 
 !
@@ -177 +192 @@
              sums_l(:,2,0) = sums_l(:,2,0) + sums_l(:,2,i)
              sums_l(:,4,0) = sums_l(:,4,0) + sums_l(:,4,i)
+             IF ( ocean )  THEN
+                sums_l(:,23,0) = sums_l(:,23,0) + sums_l(:,23,i)
+             ENDIF
              IF ( humidity )  THEN
                 sums_l(:,41,0) = sums_l(:,41,0) + sums_l(:,41,i)
@@ -200 +218 @@
        CALL MPI_ALLREDUCE( sums_l(nzb,4,0), sums(nzb,4), nzt+2-nzb, MPI_REAL, &
                            MPI_SUM, comm2d, ierr )
+       IF ( ocean )  THEN
+          CALL MPI_ALLREDUCE( sums_l(nzb,23,0), sums(nzb,23), nzt+2-nzb, &
+                              MPI_REAL, MPI_SUM, comm2d, ierr )
+       ENDIF
        IF ( humidity ) THEN
           CALL MPI_ALLREDUCE( sums_l(nzb,44,0), sums(nzb,44), nzt+2-nzb, &
@@ -221 +243 @@
        sums(:,2) = sums_l(:,2,0)
        sums(:,4) = sums_l(:,4,0)
+       IF ( ocean )  sums(:,23) = sums_l(:,23,0)
        IF ( humidity ) THEN
           sums(:,44) = sums_l(:,44,0)
@@ -241 +264 @@
        hom(:,1,2,sr) = sums(:,2)             ! v
        hom(:,1,4,sr) = sums(:,4)             ! pt
+
+!
+!--    Salinity
+       IF ( ocean )  THEN
+          sums(:,23) = sums(:,23) / ngp_2dh_outer(:,sr)
+          hom(:,1,23,sr) = sums(:,23)             ! sa
+       ENDIF
 
 !
@@ -377 +407 @@
 
 !
+!--             Salinity flux w"sa"
+                IF ( ocean )  THEN
+                   sums_l(k,65,tn) = sums_l(k,65,tn)                           &
+                                         - 0.5 * ( kh(k,j,i) + kh(k+1,j,i) )   &
+                                               * ( sa(k+1,j,i) - sa(k,j,i) )   &
+                                               * ddzu(k+1) * rmask(j,i,sr)
+                ENDIF
+
+!
 !--             Buoyancy flux, water flux (humidity flux) w"q"
                 IF ( humidity ) THEN
@@ -420 +459 @@
                 sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
                                     0.0 * rmask(j,i,sr)           ! v"pt"
+                IF ( ocean )  THEN
+                   sums_l(nzb,65,tn) = sums_l(nzb,65,tn) + &
+                                       saswsb(j,i) * rmask(j,i,sr)  ! w"sa"
+                ENDIF
                 IF ( humidity )  THEN
                    sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
@@ -450 +493 @@
                 sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
                                     0.0 * rmask(j,i,sr)           ! v"pt"
+                IF ( ocean )  THEN
+                   sums_l(nzt,65,tn) = sums_l(nzt,65,tn) + &
+                                       saswst(j,i) * rmask(j,i,sr)  ! w"sa"
+                ENDIF
                 IF ( humidity )  THEN
                    sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
@@ -508 +555 @@
                 sums_l(k,36,tn) = sums_l(k,36,tn) + pts**2 * w(k,j,i) * &
                                                     rmask(j,i,sr)
+
+!
+!--             Salinity flux and density (density does not belong to here,
+!--             but so far there is no suitable place to calculate)
+                IF ( ocean )  THEN
+                   pts = 0.5 * ( sa(k,j,i)   - hom(k,1,23,sr) + &
+                                 sa(k+1,j,i) - hom(k+1,1,23,sr) )
+                   sums_l(k,66,tn) = sums_l(k,66,tn) + pts * w(k,j,i) * &
+                                                       rmask(j,i,sr)
+                   sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) * &
+                                                       rmask(j,i,sr)
+                ENDIF
 
 !
@@ -739 +798 @@
        hom(:,1,21,sr) = sums(:,21)     ! w*pt*BC
        hom(:,1,22,sr) = sums(:,16) + sums(:,21)    ! wptBC
-                                       ! profiles 23-29 left empty for initial 
+                                       ! profile 24 is initial profile (sa)
+                                       ! profiles 25-29 left empty for initial 
                                        ! profiles
        hom(:,1,30,sr) = sums(:,30)     ! u*2
@@ -771 +831 @@
        hom(:,1,62,sr) = sums(:,62)     ! v*pt*
        hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
+       hom(:,1,64,sr) = sums(:,64)     ! rho
+       hom(:,1,65,sr) = sums(:,65)     ! w"sa"
+       hom(:,1,66,sr) = sums(:,66)     ! w*sa*
+       hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
 
        hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)