Changeset 1403 for trunk/LMDZ.MARS/libf/dynlonlat_phylonlat/phymars/datareadnc.F

Timestamp:

Mar 31, 2015, 3:49:07 PM (10 years ago)

Author:

emillour

Message:

All models: Reorganizing the physics/dynamics interface.

• makelmdz and makelmdz_fcm scripts adapted to handle the new directory settings

• misc: (replaces what was the "bibio" directory)
• Should only contain extremely generic (and non physics or dynamics-specific) routines
• Therefore moved initdynav.F90, initfluxsto.F, inithist.F, writedynav.F90, write_field.F90, writehist.F to "dyn3d_common"

• dynlonlat_phylonlat: (new interface directory)
• This directory contains routines relevent to physics/dynamics grid interactions, e.g. routines gr_dyn_fi or gr_fi_dyn and calfis
• Moreover the dynlonlat_phylonlat contains directories "phy*" corresponding to each physics package "phy*" to be used. These subdirectories should only contain specific interfaces (e.g. iniphysiq) or main programs (e.g. newstart)

• phy*/dyn1d: this subdirectory contains the 1D model using physics from phy*

EM

Location:

trunk/LMDZ.MARS/libf/dynlonlat_phylonlat

Files:

• . (added)
• phymars (added)
• phymars/datareadnc.F (moved) (moved from trunk/LMDZ.MARS/libf/phymars/datareadnc.F) (2 diffs)

trunk/LMDZ.MARS/libf/dynlonlat_phylonlat/phymars/datareadnc.F

@@ -228 +228 @@
       if (k.eq.4) then
 
-          call multscal(imd*jmdp1,zdata,1000.,zdata)
-          call multscal(imd,longitude,pi/180.,longitude)
-          call multscal(jmdp1,latitude,pi/180.,latitude)
+          zdata(:)=1000.*zdata(:)
+          longitude(:)=(pi/180.)*longitude(:)
+          latitude(:)=(pi/180.)*latitude(:)
 
           call grid_noro1(360, 180, longitude, latitude, zdata,
@@ -299 +299 @@
 c-----------------------------------------------------------------------
 
-      DO i=1,iimp1*jjp1
-            phisinit(i)=1000.*phisinit(i)
-      ENDDO
+      phisinit(1:iimp1*jjp1)=1000.*phisinit(1:iimp1*jjp1)
       CALL dump2d(iimp1,jjp1,phisinit,'Altitude in m')
-      CALL multscal(iimp1*jjp1,phisinit,g,phisinit)
+      phisinit(:)=g*phisinit(:)
 
 c-----------------------------------------------------------------------