Changeset 123 for trunk

Timestamp:

May 19, 2011, 4:54:40 PM (14 years ago)

Author:

slebonnois

Message:

SL: encore un bug oublie dans physiq.F

File:

• trunk/libf/phytitan/physiq.F (modified) (7 diffs)

trunk/libf/phytitan/physiq.F

@@ -33 +33 @@
 c debut-----input-L-variable logique indiquant le premier passage
 c lafin-----input-L-variable logique indiquant le dernier passage
-c rjourvrai-input-R-valeur réelle du jour (NBjours.gmtime)
+c rjourvrai-input-R-NBjours
 c gmtime----input-R-temps universel dans la journee (fraction de jour)
 c pdtphys---input-R-pas d'integration pour la physique (seconde)
@@ -814 +814 @@
 c Pour TITAN:
 c  calcul de la longitude solaire
-          CALL solarlong(rjourvrai,zls)
+          CALL solarlong(rjourvrai+gmtime,zls)
           print*,'Ls',zls*180./RPI      ! zls est en radians !!
       CALL orbite(zls,dist,pdecli) 
@@ -826 +826 @@
       ENDIF
 
-
-
-      
 c====================================================================
 c COUCHE LIMITE 
@@ -905 +902 @@
       ENDDO
       ENDDO
-c
+
 c        print*,"d_t_vdf1=",d_t_vdf(1,:)*dtime
 c        print*,"d_t_vdf2=",d_t_vdf(klon/2,:)*dtime
@@ -926 +923 @@
       ENDDO
       endif
-      
-            
+
       IF (if_ebil.ge.2) THEN 
         ztit='after clmain'
@@ -1056 +1052 @@
       itapchim = itapchim + 1
 
-
       if (iflag_trac.eq.1) then
          call phytrac (debut,lafin,
@@ -1074 +1069 @@
      .                        + d_tr_mph(:,:,nmicro+1:nqmax)*dtime
         endif
+
       endif