Ignore:
Timestamp:
Oct 9, 2024, 2:06:39 PM (6 weeks ago)
Author:
afalco
Message:

Pluto PCM: added conduction & molvis.
AF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.PLUTO/libf/phypluto/tabfi_mod.F90

    r3453 r3455  
    1212!   C. Hourdin 15/11/96
    1313!
    14 !   Object:        Lecture du tab_cntrl physique dans un fichier 
     14!   Object:        Lecture du tab_cntrl physique dans un fichier
    1515!   ------            et initialisation des constantes physiques
    1616!
    … …  
    2121!     ------
    2222!
    23 !      - nid:    unitne logique du fichier ou on va lire le tab_cntrl   
    24 !                      (ouvert dans le programme appellant) 
     23!      - nid:    unitne logique du fichier ou on va lire le tab_cntrl
     24!                      (ouvert dans le programme appellant)
    2525!
    2626!                 si nid=0:
    2727!                       pas de lecture du tab_cntrl mais
    2828!                       Valeurs par default des constantes physiques
    29 !       
    30 !      - tab0:    Offset de tab_cntrl a partir duquel sont ranges 
     29!
     30!      - tab0:    Offset de tab_cntrl a partir duquel sont ranges
    3131!                  les parametres physiques (50 pour start_archive)
    3232!
    … …  
    4444!      - p_rad
    4545!      - p_omeg   !
    46 !      - p_g      ! Constantes physiques ayant des 
     46!      - p_g      ! Constantes physiques ayant des
    4747!      - p_mugaz  ! homonymes dynamiques
    4848!      - p_daysec !
    … …  
    5858      use mod_phys_lmdz_para, only: is_parallel
    5959      use planete_mod, only: year_day, periastr, apoastr, peri_day, &
    60                              obliquit, z0, lmixmin, emin_turb
     60                             obliquit, z0, lmixmin, emin_turb, &
     61                             tpal, adjust
    6162      use comcstfi_mod, only: rad, omeg, g, mugaz, rcp, cpp, r
    6263      use time_phylmdz_mod, only: dtphys, daysec
    6364      use callkeys_mod, only: cpp_mugaz_mode
    6465      implicit none
    65  
     66
    6667      include "netcdf.inc"
    6768
    … …  
    8788      LOGICAL :: found
    8889      CHARACTER(len=5) :: modname="tabfi"
    89      
     90
    9091      write(*,*)"tabfi: nid=",nid," tab0=",tab0," Lmodif=",Lmodif
    9192
    … …  
    9798        lmax=0 ! not used anyways
    9899        !day_ini already set via inifis
    99         time=0 
     100        time=0
    100101! Informations about planet for dynamics and physics
    101102        ! rad,cpp,g,r,rcp already initialized by inifis
    … …  
    219220      p_rad=rad
    220221
    221       ENDIF    ! end of (nid = 0) 
     222      ENDIF    ! end of (nid = 0)
    222223
    223224!-----------------------------------------------------------------------
    224225!       Write physical constants to output before modifying them
    225226!-----------------------------------------------------------------------
    226  
     227
    227228   6  FORMAT(a20,e15.6,e15.6)
    228229   5  FORMAT(a20,f12.2,f12.2)
    229  
     230
    230231      write(*,*) '*****************************************************'
    231232      write(*,*) 'Reading tab_cntrl when calling tabfi before changes'
    … …  
    306307      write(*,*) '(10)     daysec   : length of a sol (s)'
    307308      write(*,*)
    308  
    309  
     309
     310
    310311      do while(modif(1:1).ne.'hello')
    311312        write(*,*)
    … …  
    316317        read(*,fmt='(a20)') modif
    317318        if (modif(1:1) .eq. ' ') goto 999
    318  
     319
    319320        write(*,*)
    320321        write(*,*) modif(1:len_trim(modif)) , ' : '
    … …  
    365366 106      read(*,*,iostat=ierr) emisice(1)
    366367          if(ierr.ne.0) goto 106
    367           write(*,*) 
     368          write(*,*)
    368369          write(*,*) ' emisice(1) (new value):',emisice(1)
    369370          write(*,*)
    … …  
    373374 107      read(*,*,iostat=ierr) emisice(2)
    374375          if(ierr.ne.0) goto 107
    375           write(*,*) 
     376          write(*,*)
    376377          write(*,*) ' emisice(2) (new value):',emisice(2)
    377378
    … …  
    381382 110      read(*,*,iostat=ierr) iceradius(1)
    382383          if(ierr.ne.0) goto 110
    383           write(*,*) 
     384          write(*,*)
    384385          write(*,*) ' iceradius(1) (new value):',iceradius(1)
    385386          write(*,*)
    … …  
    389390 111      read(*,*,iostat=ierr) iceradius(2)
    390391          if(ierr.ne.0) goto 111
    391           write(*,*) 
     392          write(*,*)
    392393          write(*,*) ' iceradius(2) (new value):',iceradius(2)
    393394
    … …  
    397398 112      read(*,*,iostat=ierr) dtemisice(1)
    398399          if(ierr.ne.0) goto 112
    399           write(*,*) 
     400          write(*,*)
    400401          write(*,*) ' dtemisice(1) (new value):',dtemisice(1)
    401402          write(*,*)
    … …  
    405406 113      read(*,*,iostat=ierr) dtemisice(2)
    406407          if(ierr.ne.0) goto 113
    407           write(*,*) 
     408          write(*,*)
    408409          write(*,*) ' dtemisice(2) (new value):',dtemisice(2)
    409410
    … …  
    414415 115      read(*,*,iostat=ierr) obliquit
    415416          if(ierr.ne.0) goto 115
    416           write(*,*) 
     417          write(*,*)
    417418          write(*,*) ' obliquit (new value):',obliquit
    418419
    … …  
    423424 116      read(*,*,iostat=ierr) peri_day
    424425          if(ierr.ne.0) goto 116
    425           write(*,*) 
     426          write(*,*)
    426427          write(*,*) ' peri_day (new value):',peri_day
    427428
    … …  
    432433 117      read(*,*,iostat=ierr) periastr
    433434          if(ierr.ne.0) goto 117
    434           write(*,*) 
     435          write(*,*)
    435436          write(*,*) ' periastr (new value):',periastr
    436  
     437
    437438        else if (modif(1:len_trim(modif)) .eq. 'apoastr') then
    438439          write(*,*) 'current value:',apoastr
    … …  
    441442 118      read(*,*,iostat=ierr) apoastr
    442443          if(ierr.ne.0) goto 118
    443           write(*,*) 
     444          write(*,*)
    444445          write(*,*) ' apoastr (new value):',apoastr
    445  
     446
    446447        else if (modif(1:len_trim(modif)) .eq. 'volcapa') then
    447448          write(*,*) 'current value:',volcapa
    … …  
    449450 119      read(*,*,iostat=ierr) volcapa
    450451          if(ierr.ne.0) goto 119
    451           write(*,*) 
     452          write(*,*)
    452453          write(*,*) ' volcapa (new value):',volcapa
    453        
     454
    454455        else if (modif(1:len_trim(modif)).eq.'rad') then
    455456          write(*,*) 'current value:',rad
    … …  
    457458 120      read(*,*,iostat=ierr) rad
    458459          if(ierr.ne.0) goto 120
    459           write(*,*) 
     460          write(*,*)
    460461          write(*,*) ' rad (new value):',rad
    461462
    … …  
    465466 121      read(*,*,iostat=ierr) omeg
    466467          if(ierr.ne.0) goto 121
    467           write(*,*) 
     468          write(*,*)
    468469          write(*,*) ' omeg (new value):',omeg
    469        
     470
    470471        else if (modif(1:len_trim(modif)).eq.'g') then
    471472          write(*,*) 'current value:',g
    … …  
    473474 122      read(*,*,iostat=ierr) g
    474475          if(ierr.ne.0) goto 122
    475           write(*,*) 
     476          write(*,*)
    476477          write(*,*) ' g (new value):',g
    477478
    … …  
    481482 123      read(*,*,iostat=ierr) mugaz
    482483          if(ierr.ne.0) goto 123
    483           write(*,*) 
     484          write(*,*)
    484485          write(*,*) ' mugaz (new value):',mugaz
    485486          r=8.314511/(mugaz/1000.0)
    … …  
    491492 124      read(*,*,iostat=ierr) rcp
    492493          if(ierr.ne.0) goto 124
    493           write(*,*) 
     494          write(*,*)
    494495          write(*,*) ' rcp (new value):',rcp
    495496          r=8.314511/(mugaz/1000.0)
    … …  
    502503          call su_gases
    503504          call calc_cpp_mugaz
    504           write(*,*) 
     505          write(*,*)
    505506          write(*,*) ' cpp (new value):',cpp
    506507          write(*,*) ' mugaz (new value):',mugaz
    … …  
    508509          rcp=r/cpp
    509510          write(*,*) ' rcp (new value):',rcp
    510          
     511
    511512        else if (modif(1:len_trim(modif)).eq.'daysec') then
    512513          write(*,*) 'current value:',daysec
    … …  
    514515 125      read(*,*,iostat=ierr) daysec
    515516          if(ierr.ne.0) goto 125
    516           write(*,*) 
     517          write(*,*)
    517518          write(*,*) ' daysec (new value):',daysec
    518519
    … …  
    520521        else if (modif(1:len_trim(modif)).eq.'year_day') then
    521522          write(*,*) 'current value:',year_day
    522           write(*,*) 'enter new value:' 
     523          write(*,*) 'enter new value:'
    523524 126      read(*,*,iostat=ierr) year_day
    524525          if(ierr.ne.0) goto 126
    … …  
    534535!       Write values of physical constants after modifications
    535536!-----------------------------------------------------------------------
    536  
     537
    537538      write(*,*) '*****************************************************'
    538539      write(*,*) 'Reading tab_cntrl when calling tabfi AFTER changes'
    … …  
    548549      write(*,5) '(9)             rcp',tab_cntrl(tab0+9),rcp
    549550      write(*,6) '(11)        dtphys?',tab_cntrl(tab0+11),dtphys
    550  
     551
    551552      write(*,5) '(14)       year_day',tab_cntrl(tab0+14),year_day
    552553      write(*,5) '(15)       periastr',tab_cntrl(tab0+15),periastr
    … …  
    554555      write(*,5) '(17)       peri_day',tab_cntrl(tab0+17),peri_day
    555556      write(*,5) '(18)       obliquit',tab_cntrl(tab0+18),obliquit
    556  
     557
    557558      write(*,6) '(19)             z0',tab_cntrl(tab0+19),z0
    558559      write(*,6) '(21)      emin_turb',tab_cntrl(tab0+21),emin_turb
    559560      write(*,5) '(20)        lmixmin',tab_cntrl(tab0+20),lmixmin
    560  
     561
    561562      write(*,5) '(26)        emissiv',tab_cntrl(tab0+26),emissiv
    562563      write(*,5) '(24)     emisice(1)',tab_cntrl(tab0+24),emisice(1)
    … …  
    566567      write(*,5) '(33)   dtemisice(1)',tab_cntrl(tab0+33),dtemisice(1)
    567568      write(*,5) '(34)   dtemisice(2)',tab_cntrl(tab0+34),dtemisice(2)
    568  
     569
    569570      write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
    570571
    571       write(*,*) 
    572       write(*,*) 
     572      write(*,*)
     573      write(*,*)
    573574
    574575      ENDIF ! of if (Lmodif == 1)
Note: See TracChangeset for help on using the changeset viewer.