Changeset 1642


Ignore:
Timestamp:
Dec 13, 2016, 11:38:35 AM (8 years ago)
Author:
slebonnois
Message:

SL: implementation of XIOS outputs + removal of deltatemp (chemistry) + bug correction in radlwsw

Location:
trunk/LMDZ.VENUS
Files:
4 edited
3 copied

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml

    r1637 r1642  
    2626        <grid id="grid_3d">
    2727            <domain id="dom_glo" />
    28                 <!-- <axis id="presnivs" /> -->
    29                 <axis id="altitude" />
     28                <axis id="presnivs" />
     29                <!-- <axis id="altitude" /> -->
    3030        </grid>
    3131        <!-- to output scalars -->
    … …  
    4040                      enabled=".TRUE.">
    4141       
    42         <!-- 0D variables -->
     42        <!-- 0D variables
    4343        <field_group id="fields_0D" grid_ref="scalarpoint">
    4444            <field id="ls"
    4545                   long_name="Solar longitude"
    4646                   unit="degrees" />
    47         </field_group>   
     47        </field_group>    -->
    4848       
    4949        <!-- 2D variables -->
    5050        <field_group id="fields_2D" domain_ref="dom_glo">
    51             <field id="ps"
     51            <field id="psol"
    5252                   long_name="Surface Pressure"
    5353                   unit="Pa" />
    54             <field id="area"
     54            <field id="aire"
    5555                   long_name="Mesh area"
    5656                   unit="m2" />
    … …  
    5959        <!-- 3D variables -->
    6060        <field_group id="fields_3D" grid_ref="grid_3d">
    61             <field id="temperature"
     61            <field id="temp"
    6262                   long_name="Atmospheric temperature"
    6363                   unit="K" />
    64             <field id="u"
     64            <field id="vitu"
    6565                   long_name="Eastward Zonal Wind"
    6666                   unit="m/s" />
    67             <field id="v"
     67            <field id="vitv"
    6868                   long_name="Northward Meridional Wind"
    6969                   unit="m/s" />
    … …  
    8181              enabled=".true.">
    8282                   
    83             <!-- VARS 0D -->
     83            <!-- VARS 0D
    8484            <field_group operation="instant"
    8585                         freq_op="1ts">
    8686                <field field_ref="ls" />
    87             </field_group>
     87            </field_group> -->
    8888
    8989            <!-- VARS 2D -->
    9090            <field_group operation="instant"
    9191                         freq_op="1ts">
    92                 <field field_ref="area" operation="once" />
    93                 <field field_ref="ps" />
     92                <field field_ref="aire" operation="once" />
     93                <field field_ref="psol" />
    9494            </field_group>
    9595
    … …  
    9797            <field_group operation="instant"
    9898                         freq_op="1ts">
    99                 <field field_ref="temperature" />
    100                 <field field_ref="u" />
    101                 <field field_ref="v" />
     99                <field field_ref="temp" />
     100                <field field_ref="vitu" />
     101                <field field_ref="vitv" />
    102102            </field_group>
    103103        </file>
    … …  
    108108              output_freq="1d"
    109109              type="one_file"
     110              output_level="2"
    110111              enabled=".true.">
    111112                   
    112             <!-- VARS 0D -->
     113            <!-- VARS 0D
    113114            <field_group operation="average"
    114115                         freq_op="1ts">
    115116                <field field_ref="ls" />
    116             </field_group>
     117            </field_group> -->
    117118
    118119            <!-- VARS 2D -->
    119120            <field_group operation="average"
    120121                         freq_op="1ts">
    121                 <field field_ref="area" operation="once" />
    122                 <field field_ref="ps" />
     122                <field field_ref="aire" level="1" operation="once" />
     123                <field field_ref="psol" level="1" />
    123124            </field_group>
    124125
    … …  
    126127            <field_group operation="average"
    127128                         freq_op="1ts">
    128                 <field field_ref="temperature" />
    129                 <field field_ref="u" />
    130                 <field field_ref="v" />
     129                <field field_ref="temp" level="1" />
     130                <field field_ref="vitu" level="1" />
     131                <field field_ref="vitv" level="1" />
    131132            </field_group>
    132133        </file>

  • trunk/LMDZ.VENUS/deftank/iodef.xml

    r1637 r1642  
    1111                <variable_group id="parameters" >
    1212                    <variable id="using_server" type="bool">false</variable>
    13                     <variable id="info_level" type="int">100</variable>
     13                    <variable id="info_level" type="int">10</variable>
    1414                </variable_group>
    1515            </variable_definition>

  • trunk/LMDZ.VENUS/libf/phyvenus/clesphys.h

    r1442 r1642  
    1212       LOGICAL ok_kzmin
    1313       LOGICAL callnlte,callnirco2,callthermos
    14        LOGICAL ok_cloud, ok_chem, reinit_trac, ok_sedim, ok_deltatemp
     14       LOGICAL ok_cloud, ok_chem, reinit_trac, ok_sedim
    1515       INTEGER nbapp_rad, nbapp_chim, iflag_con, iflag_ajs
    1616       INTEGER lev_histins, lev_histday, lev_histmth
    … …  
    2727     &     ok_orodr, ok_orolf, ok_gw_nonoro, ok_kzmin,                  &
    2828     &     callnlte,callnirco2,callthermos,                             &
    29      &     ok_cloud, ok_chem, reinit_trac, ok_sedim, ok_deltatemp
     29     &     ok_cloud, ok_chem, reinit_trac, ok_sedim
    3030
    3131       COMMON/clesphys_i/ nbapp_rad, nbapp_chim,                        &

  • trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90

    r1442 r1642  
    116116!Config Key  = OK_mensuel
    117117!Config Desc = Pour des sorties mensuelles
    118 !Config Def  = .true.
     118!Config Def  = .false.
    119119!Config Help = Pour creer le fichier histmth contenant les sorties
    120120!              mensuelles
    121121!
    122   ok_mensuel = .true.
     122  ok_mensuel = .false.
    123123  call getin('OK_mensuel', ok_mensuel)
    124124!
    … …  
    361361  ok_sedim = .FALSE.
    362362  call getin('ok_sedim',ok_sedim)
    363 
    364 !
    365 !Config Key  = ok_deltatemp
    366 !Config Desc = 
    367 !Config Def  = .FALSE.
    368 !Config Help =
    369 !
    370   ok_deltatemp = .FALSE.
    371   call getin('ok_deltatemp',ok_deltatemp)
    372363
    373364!
    … …  
    502493  write(numout,*)' ok_chem = ',ok_chem
    503494  write(numout,*)' ok_sedim = ',ok_sedim
    504   write(numout,*)' ok_deltatemp = ',ok_deltatemp
    505495  write(numout,*)' nb_mode = ',nb_mode
    506496  write(numout,*)' callnlte = ',callnlte

  • trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

    r1621 r1642  
    7979      use mod_grid_phy_lmdz, only: nbp_lon
    8080      use infotrac_phy, only: iflag_trac, tname, ttext
     81      use vertical_layers_mod, only: pseudoalt
     82#ifdef CPP_XIOS     
     83      use xios_output_mod, only: initialize_xios_output,
     84     &                           update_xios_timestep,
     85     &                           send_xios_field
     86#endif
    8187      IMPLICIT none
    8288c======================================================================
    … …  
    333339      REAL :: d_tr(klon,klev,nqmax)
    334340
    335 c Champ de modification de la temperature par rapport a VIRAII
    336       REAL delta_temp(klon,klev)
    337 c      SAVE delta_temp
    338       REAL mat_dtemp(33,50)
    339       SAVE mat_dtemp
    340 
    341341c Variables tendance sedimentation
    342342
    … …  
    469469
    470470         CALL printflag( ok_journe,ok_instan )
     471
     472#ifdef CPP_XIOS
     473
     474         write(*,*) "physiq: call initialize_xios_output"
     475         call initialize_xios_output(rjourvrai,gmtime,pdtphys,RDAY,
     476     &                               presnivs,pseudoalt)
     477#endif
     478
    471479c
    472480c---------
    … …  
    654662      endif
    655663
    656 c======================================================================
    657 c      Lecture du fichier DeltaT
    658 c======================================================================
    659 
    660 c     ATTENTION tout ce qui suit est pour un 48*32*50
    661 
    662       if (ok_deltatemp) then
    663 
    664       print*,'lecture de VenusDeltaT.txt '
    665       open(99, form = 'formatted', status = 'old', file =
    666      &     'VenusDeltaT.dat')
    667       print*,'Ouverture de VenusDeltaT.txt '
    668 
    669       DO ilev = 1, klev
    670          read(99,'(33(1x,e13.6))') (mat_dtemp(ilat,ilev),ilat=1,33)
    671          print*,'lecture de VenusDeltaT.txt ligne:',ilev
    672       ENDDO
    673      
    674       close(99)
    675       print*,'FIN lecture de VenusDeltaT.txt ok.'
    676 
    677       DO k = 1, klev
    678       DO i = 1, klon     
    679          ilat=(latitude_deg(i)/5.625) + 17.
    680          delta_temp(i,k)=mat_dtemp(INT(ilat),k)
    681       ENDDO
    682       ENDDO
    683      
    684       endif
    685        
    686664      ENDIF ! debut
    687665c======================================================================
    688666c======================================================================
     667! ------------------------------------------------------
     668!   Initializations done at every physical timestep:
     669! ------------------------------------------------------
    689670
    690671c Mettre a zero des variables de sortie (pour securite)
    … …  
    744725     s      , fs_bound, fq_bound )
    745726      END IF
     727
     728c====================================================================
     729c XIOS outputs
     730
     731#ifdef CPP_XIOS     
     732      ! update XIOS time/calendar
     733      call update_xios_timestep
     734#endif     
    746735
    747736c====================================================================
    … …  
    878867! Case 3: Full chemistry and/or clouds
    879868
    880          if (ok_deltatemp) then
    881 !           PRINT*,'Def de delta_temp'
    882            DO k = 1, klev
    883            DO i = 1, klon     
    884              ilat=(latitude_deg(i)/5.625) + 17.
    885 !             PRINT*,INT(ilat),latitude_deg(i),mat_dtemp(INT(ilat),k)
    886              delta_temp(i,k)=mat_dtemp(INT(ilat),k)
    887            ENDDO
    888            ENDDO
    889      
    890          endif
    891 
    892 !         Pas appel concentrations car BUG ... comprend pas
    893 !         valeur de mmean trop hautes 5,25E+29
    894 !         call concentrations2(pplay,t_seri,d_t,tr_seri, nqmax)
    895 
    896          if (ok_deltatemp) then
    897 ! Utilisation du champ de temperature modifie         
    898            call phytrac_chimie(                             
    899      I             debut,
    900      I             gmtime,
    901      I             nqmax,
    902      I             klon,
    903      I             latitude_deg,
    904      I             longitude_deg,
    905      I             nlev,
    906      I             dtime,
    907      I             t_seri+delta_temp,
    908      I             pplay,
    909      O             tr_seri)
    910          else
    911          
    912869           call phytrac_chimie(                             
    913870     I             debut,
    … …  
    922879     I             pplay,
    923880     O             tr_seri)
    924          endif
    925881
    926882c        CALL WriteField_phy('Pression',pplay,nlev)
    … …  
    935891         if (ok_sedim) then
    936892
    937          if (ok_deltatemp) then
    938 ! Utilisation du champ de temperature modifie 
    939            CALL new_cloud_sedim(
    940      I                 klon,
    941      I                     nlev,
    942      I                     dtime,
    943      I                 pplay,
    944      I                     paprs,
    945      I                     t_seri+delta_temp,
    946      I                 tr_seri,
    947      O                     d_tr_sed,
    948      O                     d_tr_ssed,
    949      I                     nqmax,
    950      O                 Fsedim)
    951          else
    952          
    953893           CALL new_cloud_sedim(
    954894     I                 klon,
    … …  
    963903     I                     nqmax,
    964904     O                 Fsedim)
    965 
    966          endif
    967905
    968906        DO k = 1, klev
    … …  
    17421680#endif
    17431681
     1682#endif
     1683
     1684! XIOS outputs
     1685! This can be done ANYWHERE in the physics routines !
     1686
     1687#ifdef CPP_XIOS     
     1688! Send fields to XIOS: (NB these fields must also be defined as
     1689! <field id="..." /> in context_lmdz_physics.xml to be correctly used)
     1690     
     1691      CALL send_xios_field("phis",phis)
     1692      cell_area_out(:)=cell_area(:)
     1693      if (is_north_pole_phy) cell_area_out(1)=cell_area(1)/nbp_lon
     1694      if (is_south_pole_phy) cell_area_out(klon)=cell_area(klon)/nbp_lon
     1695      CALL send_xios_field("aire",cell_area_out)
     1696      CALL send_xios_field("tsol",ftsol)
     1697      CALL send_xios_field("psol",paprs(:,1))
     1698      CALL send_xios_field("ue",ue)
     1699c VENUS: regardee a l envers!!!!!!!!!!!!!!!
     1700      CALL send_xios_field("ve",-1.*ve)
     1701
     1702      CALL send_xios_field("temp",t_seri)
     1703      CALL send_xios_field("vitu",u_seri)
     1704c VENUS: regardee a l envers!!!!!!!!!!!!!!!
     1705      CALL send_xios_field("vitv",-1.*v_seri)
     1706     
    17441707#endif
    17451708

  • trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F

    r1639 r1642  
    401401       if (isza0.eq.-99) then
    402402         mat0  = indexve(lat) +ips0
    403          mat0A = indexve(latA)+ips0A
    404          mat0B = indexve(latB)+ips0B
     403         if (latA.gt.0) then
     404           mat0A = indexve(latA)+ips0A
     405           mat0B = indexve(latB)+ips0B
     406         endif
    405407       else
    406408         mat0  = indexve(lat) +(isza0 -1)*nbpsve(lat) +ips0
    407          mat0A = indexve(latA)+(isza0A-1)*nbpsve(latA)+ips0A
    408          mat0B = indexve(latB)+(isza0B-1)*nbpsve(latB)+ips0B
     409         if (latA.gt.0) then
     410           mat0A = indexve(latA)+(isza0A-1)*nbpsve(latA)+ips0A
     411           mat0B = indexve(latB)+(isza0B-1)*nbpsve(latB)+ips0B
     412         endif
    409413       endif
    410414

  • trunk/LMDZ.VENUS/libf/phyvenus/xios_output_mod.F90

    r1641 r1642  
    8080    endif
    8181
    82     ! Initialize the XIOS domain coreesponding to this process:
     82    ! Initialize the XIOS domain corresponding to this process:
    8383    if (prt_level>=10) write(lunout,*) "initialize_xios_output: call wxios_domain_param"
    8484    CALL wxios_domain_param("dom_glo", is_sequential, nbp_lon, jj_nb, nbp_lon, nbp_lat, &
    8585                            1, nbp_lon, ii_begin, ii_end, jj_begin, jj_end,             &
    8686                            klon_mpi+2*(nbp_lon-1), data_ibegin, data_iend,             &
    87                             lat_reg*(180./pi), lon_reg*(180./pi),                       &
     87! VENUS IS SEEN UPSIDE DOWN, SO CENTRAL SYMMETRY TO PUT NORTH UP AGAIN
     88                           -1.*lat_reg*(180./pi), -1.*lon_reg*(180./pi),                &
    8889                            is_south_pole_dyn,mpi_rank)
    8990
    … …  
    100101                              start_date=start_date, &
    101102                              time_origin=time_origin, &
    102                               month_lengths=[30,30,30,30,30,30,30,30,30,30,30,30])
     103                              month_lengths=[2])
    103104   
    104105    ! 4. Finalize the context:
Note: See TracChangeset for help on using the changeset viewer.