Changeset 1060 for LMDZ4/trunk/libf


Ignore:
Timestamp:
Dec 16, 2008, 7:08:30 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Correction bug cas physique Newtonnienne (iflag_phys=2) pour ne plus
passer par la physique
EM/FH/IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r1007 r1060  
    1 !
    2 ! $Header$
    31!
    42c
     
    271269         IF( MOD(itau,idissip ).EQ.0.AND..NOT.forward ) apdiss = .TRUE.
    272270         IF( MOD(itau,iphysiq ).EQ.0.AND..NOT.forward
    273      s          .and. iflag_phys.NE.0                 ) apphys = .TRUE.
     271     s          .and. iflag_phys.EQ.1                 ) apphys = .TRUE.
    274272      ELSE
    275273         IF( MOD(itau   ,iconser) .EQ. 0              ) conser = .TRUE.
    276274         IF( MOD(itau+1,idissip)  .EQ. 0              ) apdiss = .TRUE.
    277          IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.NE.0) apphys=.TRUE.
     275         IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.EQ.1) apphys=.TRUE.
    278276      END IF
    279277
     
    396394     e  , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
    397395      ENDIF
    398 #else
    399 
    400 c   Calcul academique de la physique = Rappel Newtonien + fritcion
     396#endif
     397       ENDIF ! of IF( apphys )
     398
     399      IF(iflag_phys.EQ.2) THEN ! "Newtonian physics" case
     400c   Calcul academique de la physique = Rappel Newtonien + friction
    401401c   --------------------------------------------------------------
    402402       teta(:,:)=teta(:,:)
    403403     s  -iphysiq*dtvr*(teta(:,:)-tetarappel(:,:))/taurappel
    404404       call friction(ucov,vcov,iphysiq*dtvr)
    405 
    406 #endif
     405      ENDIF
     406
    407407
    408408c-jld
    409        ENDIF
    410409
    411410        CALL pression ( ip1jmp1, ap, bp, ps, p                  )
Note: See TracChangeset for help on using the changeset viewer.