Changeset 2232


Ignore:
Timestamp:
Mar 13, 2015, 9:21:25 AM (9 years ago)
Author:
fhourdin
Message:

1) Fusion des procédures clcdrag.F90 et coefcdrag.F90

dans une procédure unique cdrag.F90.
Les deux anciennes sont obsolètes mais maintenues dans le
code quelques mois pour d'éventuels tests.

2) modification du nom du module arth.F90 en arth_m.F90 pour

avoir le même nom de fichier et de module (règle adoptée
pour LMDZ et utilisée par makelmdz).

1) merging clcdrag.F90 and coefcdrag.F90 in one procedure cdrag.F90
2) renaming arth.F90 in arth_m.F90

  1. Wang/FH
Location:
LMDZ5/trunk/libf
Files:
2 added
1 deleted
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/clcdrag.F90

    r2126 r2232  
    6161  REAL, DIMENSION(klon) :: zgeop1       ! geopotentiel au 1er niveau du modele
    6262  LOGICAL, PARAMETER    :: zxli=.FALSE. ! calcul des cdrags selon Laurent Li
     63
     64  CHARACTER (LEN=80) :: abort_message
     65  CHARACTER (LEN=20) :: modname = 'clcdrag'
     66
     67
    6368!
    6469! Fonctions thermodynamiques et fonctions d'instabilite
    … …  
    6671  fsta(x) = 1.0 / (1.0+10.0*x*(1+8.0*x))
    6772  fins(x) = SQRT(1.0-18.0*x)
     73
     74  abort_message='obsolete, remplace par cdrag, use at you own risk'
     75  CALL abort_gcm(modname,abort_message,1)
     76
     77
    6878
    6979! ================================================================= c

  • LMDZ5/trunk/libf/phylmd/coefcdrag.F90

    r2158 r2232  
    6464      REAL, dimension(klon) :: trm0, trm1
    6565
     66  CHARACTER (LEN=80) :: abort_message
     67  CHARACTER (LEN=20) :: modname = 'coefcdra'
     68
     69
     70!
     71
     72
    6673!-------------------------------------------------------------------------
    6774      REAL :: fsta, fins, x
    … …  
    6976      fins(x) = SQRT(1.0-18.0*x)
    7077!-------------------------------------------------------------------------
     78
     79  abort_message='obsolete, remplace par cdrag, use at you own risk'
     80  CALL abort_gcm(modname,abort_message,1)
     81
    7182!
    7283      DO i = 1, knon

  • LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

    r2227 r2232  
    553553    REAL, DIMENSION(klon)              :: tair1, qair1, tairsol
    554554    REAL, DIMENSION(klon)              :: psfce, patm
    555     REAL, DIMENSION(klon)              :: qairsol, zgeo1
     555    REAL, DIMENSION(klon)              :: qairsol, zgeo1, speed, zri1, pref !speed, zri1, pref, added by Fuxing WANG, 04/03/2015
    556556    REAL, DIMENSION(klon)              :: rugo1
    557557    REAL, DIMENSION(klon)              :: yfluxsens
    … …  
    714714    REAL, DIMENSION(klon)       :: ytrmb3_w
    715715!
    716     REAL, DIMENSION(klon)              :: uzon_x, vmer_x
     716    REAL, DIMENSION(klon)              :: uzon_x, vmer_x, speed_x, zri1_x, pref_x !speed_x, zri1_x, pref_x, added by Fuxing WANG, 04/03/2015
    717717    REAL, DIMENSION(klon)              :: zgeo1_x, tair1_x, qair1_x, tairsol_x
    718718!
    719     REAL, DIMENSION(klon)              :: uzon_w, vmer_w
     719    REAL, DIMENSION(klon)              :: uzon_w, vmer_w, speed_w, zri1_w, pref_w !speed_w, zri1_w, pref_w, added by Fuxing WANG, 04/03/2015
    720720    REAL, DIMENSION(klon)              :: zgeo1_w, tair1_w, qair1_w, tairsol_w
    721721
    … …  
    13621362!!!
    13631363!!! nrlmd & jyg les 02/05/2011, 13/06/2011, 05/02/2012
    1364         CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
    1365             yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
     1364! Faire disparaitre les lignes commentees fin 2015 (le temps des tests)
     1365!       CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
     1366!           yu(:,1), yv(:,1), yt(:,1), yq(:,1), &
     1367!           yts, yqsurf, yrugos, &
     1368!           ycdragm, ycdragh )
     1369! Fuxing WANG, 04/03/2015, replace the clcdrag by the merged version: cdrag
     1370        DO i = 1, klon
     1371           zgeo1(i) = RD * yt(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
     1372                * (ypaprs(i,1)-ypplay(i,1))
     1373           speed(i) = SQRT(yu(i,1)**2+yv(i,1)**2)
     1374        END DO
     1375        CALL cdrag(knon, nsrf, &
     1376            speed, yt(:,1), yq(:,1), zgeo1, ypaprs(:,1),&
    13661377            yts, yqsurf, yrugos, &
    1367             ycdragm, ycdragh )
     1378            ycdragm, ycdragh, zri1, pref )
     1379
    13681380! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013
    13691381     IF (ok_prescr_ust) then
    … …  
    13801392        IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh
    13811393       ELSE  !(iflag_split .eq.0)
    1382         CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
    1383             yu_x(:,1), yv_x(:,1), yt_x(:,1), yq_x(:,1), &
     1394
     1395! Faire disparaitre les lignes commentees fin 2015 (le temps des tests)
     1396!       CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
     1397!           yu_x(:,1), yv_x(:,1), yt_x(:,1), yq_x(:,1), &
     1398!           yts_x, yqsurf, yrugos, &
     1399!           ycdragm_x, ycdragh_x )
     1400! Fuxing WANG, 04/03/2015, replace the clcdrag by the merged version: cdrag
     1401        DO i = 1, knon
     1402           zgeo1_x(i) = RD * yt_x(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
     1403                * (ypaprs(i,1)-ypplay(i,1))
     1404           speed_x(i) = SQRT(yu_x(i,1)**2+yv_x(i,1)**2)
     1405        END DO
     1406        CALL cdrag(knon, nsrf, &
     1407            speed_x, yt_x(:,1), yq_x(:,1), zgeo1_x, ypaprs(:,1),&
    13841408            yts_x, yqsurf, yrugos, &
    1385             ycdragm_x, ycdragh_x )
     1409            ycdragm_x, ycdragh_x, zri1_x, pref_x )
     1410
    13861411! --- special Dice. JYG+MPL 25112013
    13871412        IF (ok_prescr_ust) then

  • LMDZ5/trunk/libf/phylmd/screenc.F90

    r1907 r2232  
    6666! Richardson at reference level
    6767!
    68       CALL coefcdrag (klon, knon, nsrf, zxli, &
     68!      CALL coefcdrag (klon, knon, nsrf, zxli, &
     69!                    speed, temp, q_zref, gref, &
     70!                    psol, ts, qsurf, rugos, &
     71!                    okri, ri1, &
     72!                    cdram, cdrah, cdran, zri1, &
     73!                    pref)
     74! Fuxing WANG, 04/03/2015, replace the coefcdrag by the merged version: cdrag
     75      CALL cdrag (knon, nsrf, &
    6976                    speed, temp, q_zref, gref, &
    7077                    psol, ts, qsurf, rugos, &
    71                     okri, ri1, &
    72                     cdram, cdrah, cdran, zri1, &
    73                     pref)
    74 !
     78                    cdram, cdrah, zri1, pref)
    7579      DO i = 1, knon
    7680        delu(i) = ustar(i)/sqrt(cdram(i))

  • LMDZ5/trunk/libf/phylmd/stdlevvar.F90

    r2126 r2232  
    9898!
    9999      okri=.FALSE.
    100       CALL coefcdrag(klon, knon, nsrf, zxli, &
    101  &                   speed, t1, q1, z1, psol, &
    102  &                   ts1, qsurf, rugos, okri, ri1,  &         
    103  &                   cdram, cdrah, cdran, zri1, pref)           
     100!      CALL coefcdrag(klon, knon, nsrf, zxli, &
     101! &                   speed, t1, q1, z1, psol, &
     102! &                   ts1, qsurf, rugos, okri, ri1,  &         
     103! &                   cdram, cdrah, cdran, zri1, pref)           
     104! Fuxing WANG, 04/03/2015, replace the coefcdrag by the merged version: cdrag
     105      CALL cdrag(knon, nsrf, &
     106 &                   speed, t1, q1, z1, &
     107 &                   psol, ts1, qsurf, rugos, &
     108 &                   cdram, cdrah, zri1, pref)
     109
    104110! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013
    105111     IF (ok_prescr_ust) then
Note: See TracChangeset for help on using the changeset viewer.