Changeset 227 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/evap.F

Timestamp:

07/16/14 18:05:01 (10 years ago)

Author:

milmd

Message:

Last LMDZ version (1315) with OpenMP directives and other stuff

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/evap.F (modified) (6 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/evap.F

@@ -1 @@
-      subroutine evap(ngrid,nq,dtime,pt, pq, pdq, pdt,  
+      subroutine evap(ngrid,nlayer,nq,dtime,pt, pq, pdq, pdt,  
      $     dqevap,dtevap, qevap, tevap) 
         
@@ -7 +7 @@
       implicit none
 
-#include "dimensions.h"
-#include "dimphys.h"
+!#include "dimensions.h"
+!#include "dimphys.h"
 #include "comcstfi.h"
 
@@ -24 +24 @@
 !==================================================================
 
-      INTEGER ngrid,nq
+      INTEGER ngrid,nlayer,nq
 
 ! Arguments:
-      REAL pt(ngrid,nlayermx) 
-      REAL pq(ngrid,nlayermx,nq) 
-      REAL pdt(ngrid,nlayermx) 
-      REAL pdq(ngrid,nlayermx,nq) 
-      REAL dqevap(ngrid,nlayermx)
-      REAL dtevap(ngrid,nlayermx)
-      REAL qevap(ngrid,nlayermx,nq)
+      REAL pt(ngrid,nlayer) 
+      REAL pq(ngrid,nlayer,nq) 
+      REAL pdt(ngrid,nlayer) 
+      REAL pdq(ngrid,nlayer,nq) 
+      REAL dqevap(ngrid,nlayer)
+      REAL dtevap(ngrid,nlayer)
+      REAL qevap(ngrid,nlayer,nq)
       REAL dtime
  
 ! Local:
-      REAL tevap(ngrid,nlayermx) 
+      REAL tevap(ngrid,nlayer) 
       REAL zlvdcp
       REAL zlsdcp
@@ -47 +47 @@
 !
 
-      DO l=1,nlayermx
+      DO l=1,nlayer
         DO ig=1,ngrid
          qevap(ig,l,igcm_h2o_vap)=pq(ig,l,igcm_h2o_vap) 
@@ -57 +57 @@
       ENDDO
 
-      DO l = 1, nlayermx  
+      DO l = 1, nlayer
         DO ig = 1, ngrid
          zlvdcp=RLVTT/RCPD!/(1.0+RVTMP2*qevap(ig,l,igcm_h2o_vap))
@@ -77 +77 @@
       ENDDO
 
-      RETURN
       END