Changeset 534 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Feb 17, 2012, 11:06:39 AM (13 years ago)

Author:

jleconte

Message:

Corrects a bug in rcm1d

File:

• trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F

@@ -29 +29 @@
 
 #include "dimensions.h"
+#include "paramet.h"
 #include "dimphys.h"
 #include "comgeomfi.h"
@@ -45 +46 @@
 #include "logic.h"
 #include "advtrac.h"
+#include "comgeom.h"
 
 c --------------------------------------------------------------
@@ -273 +275 @@
       write(*,*)" time = ",time
       time=time/24.E+0 ! convert time (hours) to fraction of sol
-      
-       ecritphy=day_step ! default value for ecritphy
-       PRINT *,'Nunber of steps between writediagfi ?'
-       call getin("ecritphy",ecritphy)
-       write(*,*) " ecritphy = ",ecritphy
 
 
@@ -291 +288 @@
       call getin("day_step",day_step)
       write(*,*) " day_step = ",day_step
+
+       
+      ecritphy=day_step ! default value for ecritphy
+      PRINT *,'Nunber of steps between writediagfi ?'
+      call getin("ecritphy",ecritphy)
+      write(*,*) " ecritphy = ",ecritphy
 
       ndt=10 ! default value for ndt
@@ -408 +411 @@
 c   la surface de chaque maille est inutile en 1D --->
       area(1)=1.E+0
+      aire(1)=area(1) !JL+EM to have access to the area in the diagfi.nc files. area in comgeomfi.h and aire in comgeom.h
 
 c   le geopotentiel au sol est inutile en 1D car tout est controle