Ignore:
Timestamp:
May 6, 2016, 12:30:29 PM (9 years ago)
Author:
emillour
Message:

All GCMs:
Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation (up to rev r2420 of LMDZ5)

  • all physics packages:
  • added module callphysiq_mod.F90 in dynphy_lonlat/phy* which contains the routine "call_physiq" which is called by calfis* and calls the physics. This way different "physiq" routine from different physics packages may be called: The calfis* routines now exposes all available fields that might be transmitted to physiq but which is actually send (ie: expected/needed by physiq) is decided in call_physiq.
  • turned "physiq.F[90]" into module "physiq_mod.F[90]" for better control of "physiq" arguments. for phyvenus/phytitan, extracted gr_fi_ecrit from physiq.F as gr_fi_ecrit.F90 (note that it can only work in serial).
  • misc:
  • updated wxios.F90 to keep up with LMDZ5 modifications.
  • dyn3d_common:
  • infotrac.F90 keep up with LMDZ5 modifications (cosmetics)
  • dyn3d:
  • gcm.F90: cosmetic cleanup.
  • leapfrog.F90: fix computation of date as function of itau.
  • dyn3dpar:
  • gcm.F: cosmetic cleanup.
  • leapfrog_p.F90: fix computation of date as function of itau.

NB: physics are given the date corresponding to the end of the
physics step.

  • dynphy_lonlat:
  • calfis.F : added computation of relative wind vorticity.
  • calfis_p.F: added computation of relative wind vorticity (input required by Earth physics)

EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified trunk/LMDZ.COMMON/libf/dyn3d_common/infotrac.F90

    r1508 r1549  
    728728
    729729    do iq=nqo+1,nqtot
    730        write(lunout,*) 'infotrac 569: iq,tnom_0(iq)=',iq,tnom_0(iq)
     730!       write(lunout,*) 'infotrac 569: iq,tnom_0(iq)=',iq,tnom_0(iq)
    731731       do phase=1,nqo   
    732732        do ixt= 1,niso_possibles   
    733733         tnom_trac=trim(tnom_0(phase))//'_'
    734734         tnom_trac=trim(tnom_trac)//trim(tnom_iso(ixt))
    735          write(*,*) 'phase,ixt,tnom_trac=',phase,ixt,tnom_trac     
     735!         write(*,*) 'phase,ixt,tnom_trac=',phase,ixt,tnom_trac     
    736736         IF (tnom_0(iq) == tnom_trac) then
    737           write(lunout,*) 'Ce traceur est un isotope'
     737!          write(lunout,*) 'Ce traceur est un isotope'
    738738          nb_iso(ixt,phase)=nb_iso(ixt,phase)+1   
    739739          nb_isoind(phase)=nb_isoind(phase)+1   
     
    742742          indnum_fn_num(ixt)=iso_indnum(iq)
    743743          phase_num(iq)=phase
    744           write(lunout,*) 'iso_num(iq)=',iso_num(iq)
    745           write(lunout,*) 'iso_indnum(iq)=',iso_indnum(iq)
    746           write(lunout,*) 'indnum_fn_num(ixt)=',indnum_fn_num(ixt)
    747           write(lunout,*) 'phase_num(iq)=',phase_num(iq)
     744!          write(lunout,*) 'iso_num(iq)=',iso_num(iq)
     745!          write(lunout,*) 'iso_indnum(iq)=',iso_indnum(iq)
     746!          write(lunout,*) 'indnum_fn_num(ixt)=',indnum_fn_num(ixt)
     747!          write(lunout,*) 'phase_num(iq)=',phase_num(iq)
    748748          goto 20
    749749         else if (iqpere(iq).gt.0) then         
    750750          if (tnom_0(iqpere(iq)) == tnom_trac) then
    751            write(lunout,*) 'Ce traceur est le fils d''un isotope'
     751!           write(lunout,*) 'Ce traceur est le fils d''un isotope'
    752752           ! c'est un traceur d'isotope
    753753           nb_traciso(ixt,phase)=nb_traciso(ixt,phase)+1
     
    756756           zone_num(iq)=nb_traciso(ixt,phase)
    757757           phase_num(iq)=phase
    758            write(lunout,*) 'iso_num(iq)=',iso_num(iq)
    759            write(lunout,*) 'phase_num(iq)=',phase_num(iq)
    760            write(lunout,*) 'zone_num(iq)=',zone_num(iq)
     758!           write(lunout,*) 'iso_num(iq)=',iso_num(iq)
     759!           write(lunout,*) 'phase_num(iq)=',phase_num(iq)
     760!           write(lunout,*) 'zone_num(iq)=',zone_num(iq)
    761761           goto 20
    762762          endif !if (tnom_0(iqpere(iq)) == trim(tnom_0(phase))//trim(tnom_iso(ixt))) then
     
    767767      enddo !do iq=1,nqtot
    768768
    769       write(lunout,*) 'iso_num=',iso_num
    770       write(lunout,*) 'iso_indnum=',iso_indnum
    771       write(lunout,*) 'zone_num=',zone_num 
    772       write(lunout,*) 'phase_num=',phase_num
    773       write(lunout,*) 'indnum_fn_num=',indnum_fn_num
     769!      write(lunout,*) 'iso_num=',iso_num
     770!      write(lunout,*) 'iso_indnum=',iso_indnum
     771!      write(lunout,*) 'zone_num=',zone_num 
     772!      write(lunout,*) 'phase_num=',phase_num
     773!      write(lunout,*) 'indnum_fn_num=',indnum_fn_num
    774774
    775775      do ixt= 1,niso_possibles 
     
    824824    ! dimensions isotopique:
    825825    ntraciso=niso*(ntraceurs_zone+1)
    826     WRITE(lunout,*) 'niso=',niso
    827     WRITE(lunout,*) 'ntraceurs_zone,ntraciso=',ntraceurs_zone,ntraciso   
     826!    WRITE(lunout,*) 'niso=',niso
     827!    WRITE(lunout,*) 'ntraceurs_zone,ntraciso=',ntraceurs_zone,ntraciso   
    828828 
    829829    ! flags isotopiques:
     
    833833        ok_isotopes=.false.
    834834    endif
    835     WRITE(lunout,*) 'ok_isotopes=',ok_isotopes
     835!    WRITE(lunout,*) 'ok_isotopes=',ok_isotopes
    836836 
    837837    if (ok_isotopes) then
     
    843843        alpha_ideal=(/1.0,1.01,1.006,1.003,1.0/)
    844844    endif !if (ok_isotopes) then 
    845     WRITE(lunout,*) 'ok_iso_verif=',ok_iso_verif
    846     WRITE(lunout,*) 'ok_init_iso=',ok_init_iso
     845!    WRITE(lunout,*) 'ok_iso_verif=',ok_iso_verif
     846!    WRITE(lunout,*) 'ok_init_iso=',ok_init_iso
    847847
    848848    if (ntraceurs_zone.gt.0) then
     
    851851        ok_isotrac=.false.
    852852    endif   
    853     WRITE(lunout,*) 'ok_isotrac=',ok_isotrac
     853!    WRITE(lunout,*) 'ok_isotrac=',ok_isotrac
    854854
    855855    ! remplissage du tableau iqiso(ntraciso,phase)
     
    875875        index_trac(:,:)=0.0
    876876    endif !if (ok_isotrac) then
    877     write(lunout,*) 'index_trac=',index_trac   
     877!    write(lunout,*) 'index_trac=',index_trac   
    878878
    879879! Finalize :
Note: See TracChangeset for help on using the changeset viewer.