Changeset 377 for codes/icosagcm/trunk/src/etat0_dcmip3.f90

Timestamp:

04/14/16 23:12:41 (8 years ago)

Author:

dubos

Message:

New : positive advection option for theta

File:

• codes/icosagcm/trunk/src/etat0_dcmip3.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/etat0_dcmip3.f90

@@ -8 +8 @@
   
   SUBROUTINE compute_etat0(ngrid,lon,lat, phis,ps,temp,ulon,ulat,geopot,q)
-    USE genmod
     USE dcmip_initial_conditions_test_1_2_3
     USE disvert_mod
@@ -27 +26 @@
     pp=peq
     DO ij=1,ngrid
-       CALL test3_gravity_wave(lon(ij),lat(ij),pp,dummy,0, &
+       CALL test3_gravity_wave(scale_factor, lon(ij),lat(ij),pp,dummy,0, &
             dummy,dummy,dummy,dummy,phis(ij),ps(ij),dummy,dummy)
     END DO
@@ -33 +32 @@
        DO ij=1,ngrid
           pp = ap(l) + bp(l)*ps(ij) ! half-layer pressure
-          CALL test3_gravity_wave(lon(ij),lat(ij),pp,zz,0, &
+          CALL test3_gravity_wave(scale_factor, lon(ij),lat(ij),pp,zz,0, &
                dummy,dummy,dummy,dummy,dummy,dummy,dummy,dummy)
           geopot(ij,l) = g*zz ! initialize geopotential for NH
@@ -41 +40 @@
        DO ij=1,ngrid
           pp = .5*(ap(l)+ap(l+1)) + .5*(bp(l)+bp(l+1))*ps(ij) ! full-layer pressure
-          CALL test3_gravity_wave(lon(ij),lat(ij),pp,dummy,0, &
+          CALL test3_gravity_wave(scale_factor, lon(ij),lat(ij),pp,dummy,0, &
                ulon(ij,l),ulat(ij,l),dummy,Temp(ij,l),dummy,dummy,dummy,dummy)
        END DO