- Timestamp:
- 09/18/14 14:54:39 (10 years ago)
- Location:
- codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inifis.F
r227 r273 88 88 r=pr 89 89 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 90 98 91 99 avocado = 6.02214179e23 ! added by RW … … 645 653 STOP 646 654 ELSE 647 write(*,*) " mugaz=",mugaz655 write(*,*) "inifis: mugaz=",mugaz 648 656 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 649 662 cpp = -99999. 650 663 PRINT *,'SPECIFIC HEAT CAPACITY in J K-1 kg-1 ?' … … 654 667 STOP 655 668 ELSE 656 write(*,*) " cpp=",cpp669 write(*,*) "inifis: cpp=",cpp 657 670 ENDIF 658 671 ! else -
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/tabfi.F
r270 r273 104 104 g=10.44 105 105 call getin_p("g",g) 106 !mugaz=2.34 !EM: does not give cpp=11500 106 !mugaz=2.34 !EM: does not give cpp=11500) 107 107 mugaz=2.53 ! with this value of mugaz, cpp=11500 108 108 call getin_p("mugaz",mugaz) 109 ! recompute r=R/mugaz to be consitent 110 r=8.3144621/(mugaz*1.e-3) 109 111 ! kappa 110 112 rcp=0.2857143 111 113 call getin_p("kappa",rcp) 112 cpp=(8.314 511/(mugaz/1000.0))/rcp114 cpp=(8.3144621/(mugaz/1000.0))/rcp 113 115 call getin_p("cpp",cpp) 114 116 ! write(*,*) "tabfi: cpp=",cpp … … 122 124 ! Orbital parameters 123 125 periastr=9.02151966094971 126 call getin_p("periastron",periastr) 124 127 apoastr=10.054479598999 128 call getin_p("apoastron",apoastr) 125 129 peri_day=19280 130 call getin_p("periastron_day",peri_day) 126 131 obliquit=26.7299995422363 127 ! Other parameters some physical paréametrizations need 132 call getin_p("obliquity",obliquit) 133 ! Other parameters some physical parametrizations need 128 134 z0=1e-2 129 135 lmixmin=30
Note: See TracChangeset
for help on using the changeset viewer.