Changeset 605 for trunk

Timestamp:

Mar 29, 2012, 10:12:30 AM (13 years ago)

Author:

acolaitis

Message:

Added a number of corrections to pbl_parameter. The routine can now be called either without thermals (it will then use log law for wind and first level temp for t_out) and with thermals and richardson surface layer (it will conduct a Monin Obukhov interpolation). WARNING: The routine default output is now T instead of TETA. (not much change at this altitude though)

File:

• trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/phymars/pbl_parameters.F

@@ -1 +1 @@
       SUBROUTINE pbl_parameters(ngrid,nlay,ps,pplay,pz0, 
      & pg,zzlay,zzlev,pu,pv,wstar_in,hfmax,zmax,pts,ph,z_out,n_out,
-     & Teta_out,u_out,ustar,tstar,L_mo,vhf,vvv)
+     & T_out,u_out,ustar,tstar,L_mo,vhf,vvv)
       IMPLICIT NONE
 !=======================================================================
@@ -67 +67 @@
 !    --------
 
-      REAL, INTENT(OUT) :: Teta_out(ngrid,n_out),u_out(ngrid,n_out)
-      REAL T_out(ngrid,n_out)
+      REAL, INTENT(OUT) :: T_out(ngrid,n_out),u_out(ngrid,n_out)
+      REAL Teta_out(ngrid,n_out)
       REAL, INTENT(OUT) :: ustar(ngrid), tstar(ngrid)
       REAL wstar(ngrid)
@@ -156 +156 @@
          IF (zu2(ig) .ne. 0.) THEN
             ! Richardson number formulation proposed by D.E. England et al. (1995)
-            rib(ig) = (pg/ph(ig,1))
+          rib(ig) = (pg/pts(ig))
      &        *sqrt(zzlev(ig,2)*pz0(ig))
      &        *(((log(zzlev(ig,2)/pz0(ig)))**2)/(log(zzlev(ig,2)/pz0t)))
@@ -203 +203 @@
         residual = abs(pz0t-pz0tcomp(ig))
         ite = ite+1
- 
+
        ENDDO  ! of while
        pz0t=0.
@@ -273 +273 @@
      &                                  )
 
-           Teta_out(ig,n)=pts(ig)+(tstar(ig)/(prandtl(ig)*karman))*(
+           Teta_out(ig,n)=MAX(pts(ig)+(tstar(ig)/(prandtl(ig)*karman))*(
      &  2.*log(1.+sqrt(1.-16.*pz0tcomp(ig)/L_mo(ig)))
      & -2.*log(1.+sqrt(1.-16.*zout/L_mo(ig)))
      & +   log(zout/pz0tcomp(ig))
-     &                                                        )
+     &                             ),ph(ig,1))
 
           ELSE
@@ -295 +295 @@
      &                                   )
 
-           Teta_out(ig,n)=pts(ig)+(tstar(ig)/(prandtl(ig)*karman))*(
+           Teta_out(ig,n)=MAX(pts(ig)+(tstar(ig)/(prandtl(ig)*karman))*(
      &     (8./L_mo(ig))*(zout-pz0tcomp(ig))
      &   + (48./(L_mo(ig))**2)*(zout**2-pz0tcomp(ig)**2)
      &   + (1280./(3.*(L_mo(ig))**3))*(zout**3-pz0tcomp(ig)**3)
      &   + log(zout/pz0tcomp(ig))
-     &                                                           )
+     &                             ),ph(ig,1))
 
           ENDIF
@@ -317 +317 @@
 
       IF ((.not.calltherm).and.(calladj)) THEN
-           Teta_out(:,n)=ph(:,1)
+       Teta_out(:,n)=ph(:,1)
+       u_out(:,n)=(sqrt(cdn(:))*sqrt(pu(:,1)*pu(:,1)+pv(:,1)*pv(:,1))
+     &                                /karman)*log(zout/pz0(:))
       ENDIF
-
               T_out(:,n) = Teta_out(:,n)*(exp(
      &   (zout/zzlay(:,1))*(log(pplay(:,1)/ps))