Changeset 227 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/forceWCfn.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/forceWCfn.F (modified) (3 diffs)

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

@@ -1 @@
-      subroutine forceWCfn(ngrid,nq,pplev,pt,dq,dqs)
+      subroutine forceWCfn(ngrid,nlayer,nq,pplev,pt,dq,dqs)
 
       USE tracer_h
@@ -18 +18 @@
 !==================================================================
 
-#include "dimensions.h"
-#include "dimphys.h"
+!#include "dimensions.h"
+!#include "dimphys.h"
 #include "comcstfi.h"
 
-      INTEGER ngrid,nq
+      INTEGER ngrid,nlayer,nq
 
       real masse, Wtot, Wdiff
 
-      real pplev(ngrid,nlayermx+1)
+      real pplev(ngrid,nlayer+1)
       real pt(ngrid)
 
       real dqs(ngrid,nq) 
-      real dq(ngrid,nlayermx,nq)
+      real dq(ngrid,nlayer,nq)
 
       integer iq, ig, ilay
@@ -37 +37 @@
         do ig=1,ngrid
            Wtot = 0.0
-           do ilay=1,nlayermx
+           do ilay=1,nlayer
               masse = (pplev(ig,ilay) - pplev(ig,ilay+1))/g
               Wtot  = Wtot + masse*dq(ig,ilay,iq)