Ignore:
Timestamp:
09/18/14 14:54:39 (10 years ago)
Author:
millour
Message:

Bug correction in the initialization of specific gas constant in Saturn physics.
EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inifis.F

    r227 r273  
    8888      r=pr 
    8989      rcp=r/cpp 
     90      ! Ehouarn debug: 
     91      write(*,*) "inifis: rad=",rad 
     92      write(*,*) "        daysec=",daysec 
     93      write(*,*) "        dtphys=",dtphys 
     94      write(*,*) "        cpp=",cpp 
     95      write(*,*) "        g=",g 
     96      write(*,*) "        r=",r 
     97      write(*,*) "        rcp=",rcp 
    9098 
    9199      avocado = 6.02214179e23   ! added by RW 
    … …  
    645653               STOP 
    646654           ELSE 
    647                write(*,*) "mugaz=",mugaz 
     655               write(*,*) "inifis: mugaz=",mugaz 
    648656           ENDIF 
     657           !Ehouarn: once mugaz has been set, r, the specific 
     658           ! gas constant must be computed 
     659           r=8.3144622/(mugaz*1.e-3) 
     660           write(*,*) "inifis: r=",r 
     661            
    649662           cpp = -99999. 
    650663           PRINT *,'SPECIFIC HEAT CAPACITY in J K-1 kg-1 ?' 
    … …  
    654667               STOP 
    655668           ELSE 
    656                write(*,*) "cpp=",cpp 
     669               write(*,*) "inifis: cpp=",cpp 
    657670           ENDIF 
    658671!         else 
Note: See TracChangeset for help on using the changeset viewer.