Changeset 2232
- Timestamp:
- Mar 13, 2015, 9:21:25 AM (10 years ago)
- Location:
- LMDZ5/trunk/libf
- Files:
-
- 2 added
- 1 deleted
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/clcdrag.F90
r2126 r2232 61 61 REAL, DIMENSION(klon) :: zgeop1 ! geopotentiel au 1er niveau du modele 62 62 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 63 68 ! 64 69 ! Fonctions thermodynamiques et fonctions d'instabilite … … 66 71 fsta(x) = 1.0 / (1.0+10.0*x*(1+8.0*x)) 67 72 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 68 78 69 79 ! ================================================================= c -
LMDZ5/trunk/libf/phylmd/coefcdrag.F90
r2158 r2232 64 64 REAL, dimension(klon) :: trm0, trm1 65 65 66 CHARACTER (LEN=80) :: abort_message 67 CHARACTER (LEN=20) :: modname = 'coefcdra' 68 69 70 ! 71 72 66 73 !------------------------------------------------------------------------- 67 74 REAL :: fsta, fins, x … … 69 76 fins(x) = SQRT(1.0-18.0*x) 70 77 !------------------------------------------------------------------------- 78 79 abort_message='obsolete, remplace par cdrag, use at you own risk' 80 CALL abort_gcm(modname,abort_message,1) 81 71 82 ! 72 83 DO i = 1, knon -
LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90
r2227 r2232 553 553 REAL, DIMENSION(klon) :: tair1, qair1, tairsol 554 554 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 556 556 REAL, DIMENSION(klon) :: rugo1 557 557 REAL, DIMENSION(klon) :: yfluxsens … … 714 714 REAL, DIMENSION(klon) :: ytrmb3_w 715 715 ! 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 717 717 REAL, DIMENSION(klon) :: zgeo1_x, tair1_x, qair1_x, tairsol_x 718 718 ! 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 720 720 REAL, DIMENSION(klon) :: zgeo1_w, tair1_w, qair1_w, tairsol_w 721 721 … … 1362 1362 !!! 1363 1363 !!! 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),& 1366 1377 yts, yqsurf, yrugos, & 1367 ycdragm, ycdragh ) 1378 ycdragm, ycdragh, zri1, pref ) 1379 1368 1380 ! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013 1369 1381 IF (ok_prescr_ust) then … … 1380 1392 IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh 1381 1393 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),& 1384 1408 yts_x, yqsurf, yrugos, & 1385 ycdragm_x, ycdragh_x ) 1409 ycdragm_x, ycdragh_x, zri1_x, pref_x ) 1410 1386 1411 ! --- special Dice. JYG+MPL 25112013 1387 1412 IF (ok_prescr_ust) then -
LMDZ5/trunk/libf/phylmd/screenc.F90
r1907 r2232 66 66 ! Richardson at reference level 67 67 ! 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, & 69 76 speed, temp, q_zref, gref, & 70 77 psol, ts, qsurf, rugos, & 71 okri, ri1, & 72 cdram, cdrah, cdran, zri1, & 73 pref) 74 ! 78 cdram, cdrah, zri1, pref) 75 79 DO i = 1, knon 76 80 delu(i) = ustar(i)/sqrt(cdram(i)) -
LMDZ5/trunk/libf/phylmd/stdlevvar.F90
r2126 r2232 98 98 ! 99 99 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 104 110 ! --- special Dice: on force cdragm ( a defaut de forcer ustar) MPL 05082013 105 111 IF (ok_prescr_ust) then
Note: See TracChangeset
for help on using the changeset viewer.