Changeset 2369

Timestamp:

Sep 23, 2015, 10:39:13 PM (9 years ago)

Author:

oboucher

Message:

A correction on recent changes to rmu0 calculation
radpas was hard-coded in one place

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1100 +1100 @@
      ! pour obtenir le meme resultat.
      dtime=pdtphys
-     radpas = NINT( 86400./dtime/nbapp_rad)
+     IF (MOD(INT(86400./dtime),nbapp_rad).EQ.0) THEN
+       radpas = NINT( 86400./dtime/nbapp_rad)
+     ELSE 
+       WRITE(lunout,*) 'le nombre de pas de temps physique doit etre un multiple de nbapp_rad'
+       WRITE(lunout,*) 'changer nbapp_rad ou alors commenter ce test mais 1+1<>2'
+       abort_message='nbre de pas de temps physique n est pas multiple de nbapp_rad'
+       call abort_physic(modname,abort_message,1)
+     ENDIF
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -1848 +1855 @@
      !  radpas-1 pas de temps suivant. Petite ruse avec MOD pour prendre en
      !  compte le premier pas de temps de la physique pendant lequel itaprad=0 
-        zdtime1=dtime*REAL(-MOD(itaprad,4)-1)      
-        zdtime2=dtime*REAL(radpas-MOD(itaprad,4)-1) 
+        zdtime1=dtime*REAL(-MOD(itaprad,radpas)-1)      
+        zdtime2=dtime*REAL(radpas-MOD(itaprad,radpas)-1) 
         CALL zenang(zlongi,jH_cur,zdtime1,zdtime2,rlat,rlon,rmu0,fract)
      !