Changeset 5218


Ignore:
Timestamp:
Sep 23, 2024, 2:23:42 PM (3 months ago)
Author:
abarral
Message:

Rename lmdz_xer.f90 to slatec_xer.f90
Rename lmdz_libmath_pch.f90 to slatec_libmath_pch.f90
Merge r5138, r5149 from trunk which had been missed in previous merges

Location:
LMDZ6/branches/Amaury_dev
Files:
6 edited
2 moved

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev

  • LMDZ6/branches/Amaury_dev/DefLists/field_def_lmdz.xml

    r5132 r5218  
    588588            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
    589589            <!-- water oxygen H216O H217O H218O -->
     590     <field_group id="iso2D_H216O" enabled="False" > <!-- H216O -->
    590591        <field id="precipH216O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
    591592        <field id="plulH216O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
    … …  
    597598        <field id="evap_oceH216O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
    598599        <field id="evap_sicH216O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
     600     </field_group>
     601     <field_group id="iso2D_H217O" enabled="False" > <!-- H217O -->
    599602        <field id="precipH217O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
    600603        <field id="plulH217O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
    … …  
    606609        <field id="evap_oceH217O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
    607610        <field id="evap_sicH217O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
     611     </field_group>
     612     <field_group id="iso2D_H218O" enabled="False" > <!-- H218O -->
    608613        <field id="precipH218O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
    609614        <field id="plulH218O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
    … …  
    615620        <field id="evap_oceH218O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
    616621        <field id="evap_sicH218O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
     622    </field_group>
    617623            <!-- water hydrogen HDO HTO -->
     624    <field_group id="iso2D_HDO" enabled="False" > <!-- HDO -->
    618625        <field id="precipHDO" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
    619626        <field id="plulHDO"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
    … …  
    625632        <field id="evap_oceHDO"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
    626633        <field id="evap_sicHDO"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
     634    </field_group>
     635    <field_group id="iso2D_HTO" enabled="False" > <!-- HTO -->
    627636        <field id="precipHTO" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
    628637        <field id="plulHTO"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
    … …  
    634643        <field id="evap_oceHTO"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
    635644        <field id="evap_sicHTO"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
     645   </field_group>
    636646          <!-- end Added SN isotopes 2D fields 07 2023 -->
    637647   </field_group>
    … …  
    896906             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
    897907             <!-- water oxygen H216O H217O H218O -->
     908     <field_group id="iso3D_H216O" enabled="False" > <!-- H216O -->       
    898909        <field id="ovapH216O"   long_name="Specific humidity" unit="kg/kg" />
    899910        <field id="oliqH216O"   long_name="Liquid water"      unit="kg/kg" />
    … …  
    907918        <field id="dqajsH216O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    908919        <field id="dqtheH216O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
     920     </field_group>
     921     <field_group id="iso3D_H217O" enabled="False" > <!-- H217O -->
    909922        <field id="ovapH217O"   long_name="Specific humidity" unit="kg/kg" />
    910923        <field id="oliqH217O"   long_name="Liquid water"      unit="kg/kg" />
    … …  
    918931        <field id="dqajsH217O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    919932        <field id="dqtheH217O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
     933     </field_group>
     934     <field_group id="iso3D_H218O" enabled="False" > <!-- H218O -->
    920935        <field id="ovapH218O"   long_name="Specific humidity" unit="kg/kg" />
    921936        <field id="oliqH218O"   long_name="Liquid water"      unit="kg/kg" />
    … …  
    929944        <field id="dqajsH218O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    930945        <field id="dqtheH218O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
     946     </field_group>
    931947            <!-- water hydrogen HDO HTO -->
     948     <field_group id="iso3D_HDO" enabled="False" > <!-- HDO -->
    932949        <field id="ovapHDO"   long_name="Specific humidity" unit="kg/kg" />
    933950        <field id="oliqHDO"   long_name="Liquid water"      unit="kg/kg" />
    … …  
    941958        <field id="dqajsHDO"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    942959        <field id="dqtheHDO"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
     960     </field_group>
     961     <field_group id="iso3D_HTO" enabled="False" > <!-- HTO -->
    943962        <field id="ovapHTO"   long_name="Specific humidity" unit="kg/kg" />
    944963        <field id="oliqHTO"   long_name="Liquid water"      unit="kg/kg" />
    … …  
    952971        <field id="dqajsHTO"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
    953972        <field id="dqtheHTO"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
     973     </field_group>
    954974             <!-- End Added SN isotopes 3D fields 07 2023 -->
    955975   </field_group>

  • LMDZ6/branches/Amaury_dev/DefLists/file_def_histday_lmdz.xml

    r5132 r5218  
    453453            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
    454454            <!-- water oxygen H216O H217O H218O -->
    455               <field_group id="dd_2D_H216O" enabled="False" > <!-- H216O -->
    456455                <field field_ref="precipH216O"     level="1" />
    457456                <field field_ref="plulH216O"       level="1" />
    458457                <field field_ref="plucH216O"       level="1" />
     458                <field field_ref="runofflandH216O" level="1" />
    459459                <field field_ref="evapH216O"       level="1" />
    460460                <field field_ref="evap_terH216O"   level="1" />
    … …  
    462462                <field field_ref="evap_oceH216O"   level="1" />
    463463                <field field_ref="evap_sicH216O"   level="1" />
    464               </field_group>
    465               <field_group id="dd_2D_H217O" enabled="False" > <!-- H217O -->
    466464                <field field_ref="precipH217O"     level="1" />
    467465                <field field_ref="plulH217O"       level="1" />
    468466                <field field_ref="plucH217O"       level="1" />
     467                <field field_ref="runofflandH217O" level="1" />
    469468                <field field_ref="evapH217O"       level="1" />
    470469                <field field_ref="evap_terH217O"   level="1" />
    … …  
    472471                <field field_ref="evap_oceH217O"   level="1" />
    473472                <field field_ref="evap_sicH217O"   level="1" />
    474               </field_group>
    475               <field_group id="dd_2D_H218O" enabled="False" > <!-- H218O -->
    476473                <field field_ref="precipH218O"     level="1" />
    477474                <field field_ref="plulH218O"       level="1" />
    478475                <field field_ref="plucH218O"       level="1" />
     476                <field field_ref="runofflandH218O" level="1" />
    479477                <field field_ref="evapH218O"       level="1" />
    480478                <field field_ref="evap_terH218O"   level="1" />
    … …  
    482480                <field field_ref="evap_oceH218O"   level="1" />
    483481                <field field_ref="evap_sicH218O"   level="1" />
    484               </field_group>
    485482            <!-- water hydrogen HDO HTO -->
    486               <field_group id="dd_2D_HDO" enabled="False" > <!-- HDO -->
    487483                <field field_ref="precipHDO"    level="1" />
    488484                <field field_ref="plulHDO"      level="1" />
    489485                <field field_ref="plucHDO"      level="1" />
     486                <field field_ref="runofflandHDO" level="1" />
    490487                <field field_ref="evapHDO"      level="1" />
    491488                <field field_ref="evap_terHDO"  level="1" />
    … …  
    493490                <field field_ref="evap_oceHDO"  level="1" />
    494491                <field field_ref="evap_sicHDO"  level="1" />
    495               </field_group>
    496               <field_group id="dd_2D_HTO" enabled="False" > <!-- HTO -->
    497492                <field field_ref="plulHTO"      level="1" />
    498493                <field field_ref="precipHTO"    level="1" />
    499494                <field field_ref="plucHTO"      level="1" />
     495                <field field_ref="runofflandHTO" level="1" />
    500496                <field field_ref="evapHTO"      level="1" />
    501497                <field field_ref="evap_terHTO"  level="1" />
    … …  
    503499                <field field_ref="evap_oceHTO"  level="1" />
    504500                <field field_ref="evap_sicHTO"  level="1" />
    505               </field_group>
    506501             <!-- end Added SN isotopes 2D fields 07 2023 -->
    507502            </field_group>
    … …  
    682677             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
    683678             <!-- water oxygen H216O H217O H218O -->
    684              <field_group id="dd_3D_H216O" enabled="False"  > <!-- H216O -->
    685679                <field field_ref="ovapH216O"   level="1" />
    686680                <field field_ref="oliqH216O"   level="1" />
    … …  
    694688                <field field_ref="dqajsH216O"  level="1" />
    695689                <field field_ref="dqtheH216O"  level="1" />
    696               </field_group>
    697               <field_group id="dd_3D_H217O" enabled="False" > <!-- H217O -->
    698690                <field field_ref="ovapH217O"   level="1" />
    699691                <field field_ref="oliqH217O"   level="1" />
    … …  
    707699                <field field_ref="dqajsH217O"  level="1" />
    708700                <field field_ref="dqtheH217O"  level="1" />
    709               </field_group>
    710               <field_group id="dd_3D_H218O" enabled="False" > <!-- H218O -->
    711701                <field field_ref="ovapH218O"   level="1" />
    712702                <field field_ref="oliqH218O"   level="1" />
    … …  
    720710                <field field_ref="dqajsH218O"  level="1" />
    721711                <field field_ref="dqtheH218O"  level="1" />
    722               </field_group>
    723712            <!-- water hydrogen HDO HTO -->
    724               <field_group id="dd_3D_HDO" enabled="False" > <!-- HDO -->
    725713                <field field_ref="ovapHDO"   level="1" />
    726714                <field field_ref="oliqHDO"   level="1" />
    … …  
    734722                <field field_ref="dqajsHDO"  level="1" />
    735723                <field field_ref="dqtheHDO"  level="1" />
    736               </field_group>
    737               <field_group id="dd_3D_HTO" enabled="False" > <!-- HTO -->
    738724                <field field_ref="ovapHTO"   level="1" />
    739725                <field field_ref="oliqHTO"   level="1" />
    … …  
    747733                <field field_ref="dqajsHTO"  level="1" />
    748734                <field field_ref="dqtheHTO"  level="1" />
    749               </field_group>
    750735             <!-- End Added SN isotopes 3D fields 07 2023 -->
    751736            </field_group>

  • LMDZ6/branches/Amaury_dev/DefLists/file_def_histmth_lmdz.xml

    r5132 r5218  
    502502            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
    503503            <!-- water oxygen H216O H217O H218O -->
    504               <field_group id="mm_2D_H216O" enabled="False" > <!-- H216O -->
    505504                <field field_ref="precipH216O"     level="1" />
    506505                <field field_ref="plulH216O"       level="1" />
    507506                <field field_ref="plucH216O"       level="1" />
     507                <field field_ref="runofflandH216O" level="1" />
    508508                <field field_ref="evapH216O"       level="1" />
    509509                <field field_ref="evap_terH216O"   level="1" />
    … …  
    511511                <field field_ref="evap_oceH216O"   level="1" />
    512512                <field field_ref="evap_sicH216O"   level="1" />
    513               </field_group>
    514               <field_group id="mm_2D_H217O" enabled="False" > <!-- H217O -->
    515513                <field field_ref="precipH217O"     level="1" />
    516514                <field field_ref="plulH217O"       level="1" />
    517515                <field field_ref="plucH217O"       level="1" />
     516                <field field_ref="runofflandH217O" level="1" />
    518517                <field field_ref="evapH217O"       level="1" />
    519518                <field field_ref="evap_terH217O"   level="1" />
    … …  
    521520                <field field_ref="evap_oceH217O"   level="1" />
    522521                <field field_ref="evap_sicH217O"   level="1" />
    523               </field_group>
    524               <field_group id="mm_2D_H218O" enabled="False" > <!-- H218O -->
    525522                <field field_ref="precipH218O"     level="1" />
    526523                <field field_ref="plulH218O"       level="1" />
    527524                <field field_ref="plucH218O"       level="1" />
     525                <field field_ref="runofflandH218O" level="1" />
    528526                <field field_ref="evapH218O"       level="1" />
    529527                <field field_ref="evap_terH218O"   level="1" />
    … …  
    531529                <field field_ref="evap_oceH218O"   level="1" />
    532530                <field field_ref="evap_sicH218O"   level="1" />
    533               </field_group>
    534531            <!-- water hydrogen HDO HTO -->
    535               <field_group id="mm_2D_HDO" enabled="False" > <!-- HDO -->
    536532                <field field_ref="precipHDO"     level="1" />
    537533                <field field_ref="plulHDO"       level="1" />
    538534                <field field_ref="plucHDO"       level="1" />
     535                <field field_ref="runofflandHDO" level="1" />
    539536                <field field_ref="evapHDO"       level="1" />
    540537                <field field_ref="evap_terHDO"   level="1" />
    … …  
    542539                <field field_ref="evap_oceHDO"   level="1" />
    543540                <field field_ref="evap_sicHDO"   level="1" />
    544               </field_group>
    545               <field_group id="mm_2D_HTO" enabled="False" > <!-- HTO -->
    546541                <field field_ref="plulHTO"       level="1" />
    547542                <field field_ref="precipHTO"     level="1" />
    548543                <field field_ref="plucHTO"       level="1" />
     544                <field field_ref="runofflandHTO" level="1" />
    549545                <field field_ref="evapHTO"       level="1" />
    550546                <field field_ref="evap_terHTO"   level="1" />
    … …  
    552548                <field field_ref="evap_oceHTO"   level="1" />
    553549                <field field_ref="evap_sicHTO"   level="1" />
    554               </field_group>
    555550         <!-- end Added SN isotopes 2D fields 07 2023 -->
    556 
    557551            </field_group>
    558552
    … …  
    738732             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
    739733             <!-- water oxygen H216O H217O H218O -->
    740              <field_group id="mm_3D_H216O" enabled="False" > <!-- H216O -->
    741734                <field field_ref="ovapH216O"   level="1" />
    742735                <field field_ref="oliqH216O"   level="1" />
    … …  
    750743                <field field_ref="dqajsH216O"  level="1" />
    751744                <field field_ref="dqtheH216O"  level="1" />
    752               </field_group>
    753               <field_group id="mm_3D_H217O" enabled="False" > <!-- H217O -->
    754745                <field field_ref="ovapH217O"   level="1" />
    755746                <field field_ref="oliqH217O"   level="1" />
    … …  
    763754                <field field_ref="dqajsH217O"  level="1" />
    764755                <field field_ref="dqtheH217O"  level="1" />
    765               </field_group>
    766               <field_group id="mm_3D_H218O" enabled="False" > <!-- H218O -->
    767756                <field field_ref="ovapH218O"   level="1" />
    768757                <field field_ref="oliqH218O"   level="1" />
    … …  
    776765                <field field_ref="dqajsH218O"  level="1" />
    777766                <field field_ref="dqtheH218O"  level="1" />
    778               </field_group>
    779767            <!-- water hydrogen HDO HTO -->
    780               <field_group id="mm_3D_HDO" enabled="False" > <!-- HDO -->
    781768                <field field_ref="ovapHDO"   level="1" />
    782769                <field field_ref="oliqHDO"   level="1" />
    … …  
    790777                <field field_ref="dqajsHDO"  level="1" />
    791778                <field field_ref="dqtheHDO"  level="1" />
    792               </field_group>
    793               <field_group id="mm_3D_HTO" enabled="False" > <!-- HTO -->
    794779                <field field_ref="ovapHTO"   level="1" />
    795780                <field field_ref="oliqHTO"   level="1" />
    … …  
    803788                <field field_ref="dqajsHTO"  level="1" />
    804789                <field field_ref="dqtheHTO"  level="1" />
    805               </field_group>
    806790             <!-- End Added SN isotopes 3D fields 07 2023 -->
    807  
    808791            </field_group>
    809792

  • LMDZ6/branches/Amaury_dev/libf/misc/slatec_libmath_pch.f90

    r5217 r5218  
    22
    33MODULE lmdz_libmath_pch
    4   USE lmdz_xer, ONLY: xermsg
     4  USE slatec_xer, ONLY: xermsg
    55  IMPLICIT NONE; PRIVATE
    66  PUBLIC pchfe_95, pchsp_95

  • LMDZ6/branches/Amaury_dev/libf/phylmd/phys_output_write_mod.F90

    r5160 r5218  
    29212921            i = INDEX(outiso, '_', .TRUE.)
    29222922            outiso = outiso(1:i-1)//outiso(i+1:LEN_TRIM(outiso))
    2923             ! dayly
    2924             CALL xios_set_fieldgroup_attr("dd_2D_"//TRIM(outiso), enabled=.TRUE.)
    2925             CALL xios_set_fieldgroup_attr("dd_3D_"//TRIM(outiso), enabled=.TRUE.)
    2926             ! monthly
    2927             CALL xios_set_fieldgroup_attr("mm_2D_"//TRIM(outiso), enabled=.TRUE.)
    2928             CALL xios_set_fieldgroup_attr("mm_3D_"//TRIM(outiso), enabled=.TRUE.)
     2923
     2924            CALL xios_set_fieldgroup_attr("iso2D_"//TRIM(outiso), enabled=.TRUE.)
     2925            CALL xios_set_fieldgroup_attr("iso3D_"//TRIM(outiso), enabled=.TRUE.)
    29292926
    29302927          ENDDO

  • LMDZ6/branches/Amaury_dev/libf/phylmdiso/reevap.F90

    r5160 r5218  
    2828    REAL, DIMENSION(klon,klev,nqtot), INTENT(IN) ::     qx
    2929    REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_t_eva
    30     REAL, DIMENSION(klon,klev,nqtot), INTENT(OUT) ::        d_qx_eva
     30    REAL, DIMENSION(klon,klev,nqtot), INTENT(OUT) ::    d_qx_eva
    3131
    3232    REAL za,zb,zdelta,zlvdcp,zlsdcp
    … …  
    3535    !--------Stochastic Boundary Layer Triggering: ALE_BL--------
    3636
    37 DO ixt=1,1+ntiso
     37DO ixt = 1, 1+ntiso
    3838    ! Re-evaporer l'eau liquide nuageuse
    3939
    … …  
    4343!PRINT *,'rrevap ; fl_cor_ebil:',fl_cor_ebil,' iflag_ice_thermo:',iflag_ice_thermo,' RVTMP2',RVTMP2
    4444    DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
    45        DO i = 1, klon
     45      DO i = 1, klon
    4646
    47         IF (ixt==1) THEN
     47        IF (ixt == 1) THEN ! water
    4848         IF (fl_cor_ebil > 0) THEN
    4949          !zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*(q_seri(i,k)+ql_seri(i,k)+qs_seri(i,k)))
    … …  
    5151          zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*(qx(i,k,ivapcur)+qx(i,k,iliqcur)+qx(i,k,isolcur)))
    5252          zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*(qx(i,k,ivapcur)+qx(i,k,iliqcur)+qx(i,k,isolcur)))
    53          else
     53         ELSE
    5454           zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*qx(i,k,ivapcur))
    5555           !jyg<
    … …  
    5757           !                  A verifier !!!
    5858           zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*qx(i,k,ivapcur))
    59          end if
     59         END IF
    6060         IF (iflag_ice_thermo == 0) THEN
    6161             zlsdcp=zlvdcp
    62          ENDIF
     62         END IF
    6363          !>jyg
     64        END IF
     65        IF (iflag_ice_thermo == 0) THEN
     66          !pas necessaire a priori
     67
     68          zdelta = MAX(0., SIGN(1., RTT - t_seri(i, k)))
     69          zdelta = 0.
     70          zb = MAX(0.0, qx(i, k, iliqcur))
     71          IF (ixt == 1) THEN
     72            za = - MAX(0.0, qx(i, k, iliqcur)) &
     73                    * (zlvdcp * (1. - zdelta) + zlsdcp * zdelta)
     74            d_t_eva(i, k) = za
     75          END IF
     76          d_qx_eva(i, k, ivapcur) = zb
     77          d_qx_eva(i, k, iliqcur) = -qx(i, k, iliqcur)
     78          d_qx_eva(i, k, isolcur) = 0.
     79
     80        ELSE
     81
    6482             !CR: on r\'e-\'evapore eau liquide et glace
    6583
    … …  
    6886             !        za = - MAX(0.0,ql_seri(i,k)) &
    6987             !             * (zlvdcp*(1.-zdelta)+zlsdcp*zdelta)
    70              za = - MAX(0.0,qx(i,k,iliqcur))*zlvdcp &
    71                   - MAX(0.0,qx(i,k,iliqcur))*zlsdcp
    72              d_t_eva(i,k) = za
     88          IF (ixt == 1) THEN
     89            za = - MAX(0.0, qx(i, k, iliqcur)) * zlvdcp &
     90                    - MAX(0.0, qx(i, k, iliqcur)) * zlsdcp
     91            d_t_eva(i, k) = za
    7392
    74           endif !if (ixt.EQ.1) THEN
     93          END IF
    7594            !zb = MAX(0.0,ql_seri(i,k)+qs_seri(i,k))
    7695            !d_q_eva(i,k) = zb
    … …  
    82101            d_qx_eva(i,k,iliqcur) = -qx(i,k,iliqcur)
    83102            d_qx_eva(i,k,isolcur) = -qx(i,k,isolcur)
     103        END IF
    84104
    85105
    86        ENDDO
    87     ENDDO
     106      END DO
     107    END DO
    88108
    89     enddo ! do ixt=1,1+niso*(nzone +1)
    90 
     109    END DO ! DO ixt = 1, 1+niso*(nzone +1)
    91110
    92111
Note: See TracChangeset for help on using the changeset viewer.