Ignore:
Timestamp:
Sep 16, 2016, 12:09:13 PM (8 years ago)
Author:
mlefevre
Message:

Added the reading of planet_constant file for prescribed planet

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/share/module_model_constants.F

    r1599 r1604  
    255255     womeg        = 0.2992549E-06
    256256     wdaysec      = 0.1008707E08
     257    ELSE if ( planet == "prescribed") then
     258     open(17,file='planet_constant',form='formatted',status='old')
     259     rewind(17)
     260     read(17,*) g
     261     read(17,*) cp
     262     read(17,*) molmas
     263     read(17,*) t0
     264     read(17,*) p0
     265     read(17,*) rad
     266     close(17)
     267     r_d=8.314/(molmas*1e-3)
     268     p1000mb=p0
     269     reradius=1./rad
    257270    Else
    258271     write(*,*) "unknown planet type:", planet
     
    275288     rfcp         =0.25/cp
    276289     p608         =rvovrd-1.
    277      z_scale      =r_d*t0/g
     290     z_scale      =(r_d*t0/g)*1e-3
    278291    write(*,*) "planet : ", planet
    279292    write(*,*) "g = ", g
Note: See TracChangeset for help on using the changeset viewer.