Changeset 1311

Timestamp:

Jul 11, 2014, 8:37:21 AM (10 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in testphys1d: routine ini_fillgeom expects arrays and not scalars as arguments, so define longitude and latitude as arrays (of 1 element).
EM

File:

• trunk/LMDZ.MARS/libf/phymars/testphys1d.F (modified) (3 diffs)

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

@@ -104 +104 @@
       Logical  tracerdyn
       integer :: nq=1 ! number of tracers
-      real :: latitude, longitude
+      real :: latitude(1), longitude(1)
 
       character*2 str2
@@ -464 +464 @@
 c  latitude/longitude
 c  ------------------
-      latitude=0 ! default value for latitude
+      latitude(1)=0 ! default value for latitude
       PRINT *,'latitude (in degrees) ?'
-      call getin("latitude",latitude)
+      call getin("latitude",latitude(1))
       write(*,*) " latitude = ",latitude
       latitude=latitude*pi/180.E+0
@@ -478 +478 @@
       call phys_state_var_init(1,llm,nq,
      .          daysec,dtphys,rad,g,r,cpp)
-      call ini_fillgeom(1,latitude,longitude,1.0)
+      call ini_fillgeom(1,latitude,longitude,(/1.0/))
       call conf_phys(1,llm,nq)