Changeset 778
- Timestamp:
- Sep 12, 2012, 3:45:33 PM (12 years ago)
- Location:
- trunk
- Files:
-
- 1 added
- 1 deleted
- 3 edited
-
LMDZ.TITAN/libf/phytitan/readstartphy.F (modified) (10 diffs)
-
LMDZ.VENUS/libf/phyvenus/YOEGWD.h (modified) (1 diff)
-
LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F (deleted)
-
LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90 (added)
-
LMDZ.VENUS/libf/phyvenus/readstartphy.F (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/readstartphy.F
r102 r778 20 20 #include "temps.h" 21 21 c====================================================================== 22 integer ngridmx 23 parameter (ngridmx=(2+(jjm-1)*iim - 1/jjm)) 24 22 25 CHARACTER*(*) fichnom 23 REAL rlat( klon), rlon(klon)24 REAL tsol( klon)25 REAL tsoil( klon,nsoilmx)26 REAL albe( klon)26 REAL rlat(ngridmx), rlon(ngridmx) 27 REAL tsol(ngridmx) 28 REAL tsoil(ngridmx,nsoilmx) 29 REAL albe(ngridmx) 27 30 cIM BEG alblw 28 REAL alblw( klon)31 REAL alblw(ngridmx) 29 32 cIM END alblw 30 REAL radsol( klon)31 REAL sollw( klon)32 real solsw( klon)33 real fder( klon)33 REAL radsol(ngridmx) 34 REAL sollw(ngridmx) 35 real solsw(ngridmx) 36 real fder(ngridmx) 34 37 INTEGER length 35 38 PARAMETER (length=100) … … 114 117 PRINT*, 'phyetat0: Le champ <TS> est present' 115 118 #ifdef NC_DOUBLE 116 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol (1))117 #else 118 ierr = NF_GET_VAR_REAL(nid, nvarid, tsol (1))119 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol) 120 #else 121 ierr = NF_GET_VAR_REAL(nid, nvarid, tsol) 119 122 #endif 120 123 IF (ierr.NE.NF_NOERR) THEN … … 124 127 xmin = 1.0E+20 125 128 xmax = -1.0E+20 126 DO i = 1, klon129 DO i = 1, ngridmx 127 130 xmin = MIN(tsol(i),xmin) 128 131 xmax = MAX(tsol(i),xmax) … … 143 146 PRINT*, "phyetat0: Le champ <Tsoil"//str2//"> est absent" 144 147 PRINT*, " Il prend donc la valeur de surface" 145 DO i=1, klon148 DO i=1, ngridmx 146 149 tsoil(i,isoil)=tsol(i) 147 150 ENDDO … … 158 161 ENDIF 159 162 ENDDO 163 160 164 c 161 165 c Lecture de albedo au sol: … … 179 183 xmin = 1.0E+20 180 184 xmax = -1.0E+20 181 DO i = 1, klon185 DO i = 1, ngridmx 182 186 xmin = MIN(albe(i),xmin) 183 187 xmax = MAX(albe(i),xmax) … … 207 211 xmin = 1.0E+20 208 212 xmax = -1.0E+20 209 DO i = 1, klon213 DO i = 1, ngridmx 210 214 xmin = MIN(solsw(i),xmin) 211 215 xmax = MAX(solsw(i),xmax) … … 233 237 xmin = 1.0E+20 234 238 xmax = -1.0E+20 235 DO i = 1, klon239 DO i = 1, ngridmx 236 240 xmin = MIN(sollw(i),xmin) 237 241 xmax = MAX(sollw(i),xmax) … … 260 264 xmin = 1.0E+20 261 265 xmax = -1.0E+20 262 DO i = 1, klon266 DO i = 1, ngridmx 263 267 xmin = MIN(fder(i),xmin) 264 268 xmax = MAX(fder(i),xmax) … … 285 289 xmin = 1.0E+20 286 290 xmax = -1.0E+20 287 DO i = 1, klon291 DO i = 1, ngridmx 288 292 xmin = MIN(radsol(i),xmin) 289 293 xmax = MAX(radsol(i),xmax) -
trunk/LMDZ.VENUS/libf/phyvenus/YOEGWD.h
r3 r778 1 C-----------------------------------------------------------------2 C* *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS3 C-----------------------------------------------------------------4 C 5 integer NKTOPG,NTOP6 real GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT7 real GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC8 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 9 *,GHMAX,GRAHILO,GSIGCR,NKTOPG,NTOP,GSSEC,GTSEC,GVSEC1 ! ----------------------------------------------------------------- 2 !* *COMMON* *YOEGWD* - PARAMETERS FOR GRAVITY WAVE DRAG CALCULATIONS 3 ! ----------------------------------------------------------------- 4 ! 5 integer :: NKTOPG,NTOP 6 real :: GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT 7 real :: GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC 8 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT & 9 ,GHMAX,GRAHILO,GSIGCR,NKTOPG,NTOP,GSSEC,GTSEC,GVSEC 10 10 -
trunk/LMDZ.VENUS/libf/phyvenus/readstartphy.F
r101 r778 20 20 #include "temps.h" 21 21 c====================================================================== 22 integer ngridmx 23 parameter (ngridmx=(2+(jjm-1)*iim - 1/jjm)) 24 22 25 CHARACTER*(*) fichnom 23 REAL rlat( klon), rlon(klon)24 REAL tsol( klon)25 REAL tsoil( klon,nsoilmx)26 REAL albe( klon)26 REAL rlat(ngridmx), rlon(ngridmx) 27 REAL tsol(ngridmx) 28 REAL tsoil(ngridmx,nsoilmx) 29 REAL albe(ngridmx) 27 30 cIM BEG alblw 28 REAL alblw( klon)31 REAL alblw(ngridmx) 29 32 cIM END alblw 30 REAL radsol( klon)31 REAL sollw( klon)32 real solsw( klon)33 real fder( klon)33 REAL radsol(ngridmx) 34 REAL sollw(ngridmx) 35 real solsw(ngridmx) 36 real fder(ngridmx) 34 37 INTEGER length 35 38 PARAMETER (length=100) … … 114 117 PRINT*, 'phyetat0: Le champ <TS> est present' 115 118 #ifdef NC_DOUBLE 116 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol (1))117 #else 118 ierr = NF_GET_VAR_REAL(nid, nvarid, tsol (1))119 ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol) 120 #else 121 ierr = NF_GET_VAR_REAL(nid, nvarid, tsol) 119 122 #endif 120 123 IF (ierr.NE.NF_NOERR) THEN … … 124 127 xmin = 1.0E+20 125 128 xmax = -1.0E+20 126 DO i = 1, klon129 DO i = 1, ngridmx 127 130 xmin = MIN(tsol(i),xmin) 128 131 xmax = MAX(tsol(i),xmax) … … 143 146 PRINT*, "phyetat0: Le champ <Tsoil"//str2//"> est absent" 144 147 PRINT*, " Il prend donc la valeur de surface" 145 DO i=1, klon148 DO i=1, ngridmx 146 149 tsoil(i,isoil)=tsol(i) 147 150 ENDDO … … 158 161 ENDIF 159 162 ENDDO 163 160 164 c 161 165 c Lecture de albedo au sol: … … 179 183 xmin = 1.0E+20 180 184 xmax = -1.0E+20 181 DO i = 1, klon185 DO i = 1, ngridmx 182 186 xmin = MIN(albe(i),xmin) 183 187 xmax = MAX(albe(i),xmax) … … 207 211 xmin = 1.0E+20 208 212 xmax = -1.0E+20 209 DO i = 1, klon213 DO i = 1, ngridmx 210 214 xmin = MIN(solsw(i),xmin) 211 215 xmax = MAX(solsw(i),xmax) … … 233 237 xmin = 1.0E+20 234 238 xmax = -1.0E+20 235 DO i = 1, klon239 DO i = 1, ngridmx 236 240 xmin = MIN(sollw(i),xmin) 237 241 xmax = MAX(sollw(i),xmax) … … 260 264 xmin = 1.0E+20 261 265 xmax = -1.0E+20 262 DO i = 1, klon266 DO i = 1, ngridmx 263 267 xmin = MIN(fder(i),xmin) 264 268 xmax = MAX(fder(i),xmax) … … 285 289 xmin = 1.0E+20 286 290 xmax = -1.0E+20 287 DO i = 1, klon291 DO i = 1, ngridmx 288 292 xmin = MIN(radsol(i),xmin) 289 293 xmax = MAX(radsol(i),xmax)
Note: See TracChangeset
for help on using the changeset viewer.
