Changeset 901 for trunk/LMDZ.VENUS/libf
- Timestamp:
- Mar 12, 2013, 7:36:18 PM (12 years ago)
- Location:
- trunk/LMDZ.VENUS/libf/phyvenus
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/libf/phyvenus/iostart.F90
r892 r901 25 25 26 26 PUBLIC get_field,get_var,put_field,put_var 27 PUBLIC Open_startphy,close_startphy,open_restartphy,close_restartphy27 PUBLIC open_startphy,close_startphy,open_restartphy,close_restartphy 28 28 29 29 CONTAINS 30 30 31 SUBROUTINE Open_startphy(filename)31 SUBROUTINE open_startphy(filename) 32 32 USE netcdf 33 33 USE mod_phys_lmdz_para … … 45 45 ENDIF 46 46 47 END SUBROUTINE Open_startphy47 END SUBROUTINE open_startphy 48 48 49 49 SUBROUTINE Close_startphy … … 326 326 ENDIF 327 327 328 ierr = NF90_PUT_ATT (nid_restart, NF90_GLOBAL, "title","Fichier redem marage physique")328 ierr = NF90_PUT_ATT (nid_restart, NF90_GLOBAL, "title","Fichier redemarrage physique") 329 329 330 330 ierr = NF90_DEF_DIM (nid_restart, "index", length, idim1) -
trunk/LMDZ.VENUS/libf/phyvenus/phys_state_var_mod.F90
r892 r901 137 137 IMPLICIT NONE 138 138 139 deallocate(ftsol, f albe)139 deallocate(ftsol, ftsoil, falbe) 140 140 deallocate(zmea, zstd, zsig, zgam) 141 141 deallocate(zthe, zpic, zval) -
trunk/LMDZ.VENUS/libf/phyvenus/rcm1d.F
r894 r901 5 5 use control_mod 6 6 use comgeomphy 7 USE phys_state_var_mod 7 8 IMPLICIT NONE 8 9 … … 48 49 REAL play(llm) ! Pressure at the middle of the layers (Pa) 49 50 REAL plev(llm+1) ! intermediate pressure levels (pa) 50 REAL psurf ,tsurf(1)51 REAL psurf 51 52 REAL u(llm),v(llm) ! zonal, meridional wind 52 53 REAL gru,grv ! prescribed "geostrophic" background wind 53 54 REAL temp(llm) ! temperature at the middle of the layers 54 55 REAL,allocatable :: q(:,:) ! tracer mixing ratio (e.g. kg/kg) 55 REAL tsoil(nsoilmx) ! subsurface soik temperature (K)56 56 REAL zlay(llm) ! altitude estimee dans les couches (km) 57 57 REAL long(1),lati(1),area(1) 58 58 REAL cufi(1),cvfi(1) 59 REAL phisfi(1),albedo(1) 60 REAL solsw(1),sollwdown(1),dlw(1),radsol(1) 61 REAL zmea(1), zstd(1) 62 REAL zsig(1), zgam(1), zthe(1) 63 REAL zpic(1), zval(1) 59 REAL phisfi(1) 64 60 65 61 c Physical and dynamical tandencies (e.g. m.s-2, K/s, Pa/s) … … 196 192 long(1)=0.E+0 197 193 198 c Initialisation albedo199 c ----------------------200 c ne sert pas ici...201 albedo(1)=0.1202 c PRINT *,'Albedo du sol nu ?'203 c READ(unit,*) albedo(1)204 c PRINT *,albedo(1)205 206 194 c Initialisation speciales "physiq" 207 195 c --------------------------------- … … 271 259 pks=cpp*(psurf/preff)**rcp 272 260 261 c init des variables pour phyredem 262 c -------------------------------- 263 call phys_state_var_init 264 273 265 c profil de temperature et altitude au premier appel 274 266 c -------------------------------------------------- … … 286 278 print*," Pression Altitude Temperature" 287 279 ilayer=1 288 tsurf(1)=tmp2(0)280 ftsol(1)=tmp2(0) 289 281 temp(1)=tmp2(1) 290 282 zlay(1)=tmp2(1)*tmp1(1) 291 print*," 0", tsurf(1)283 print*," 0",ftsol(1) 292 284 print*,ilayer,play(ilayer),zlay(ilayer),temp(ilayer) 293 285 DO ilayer=2,nlayer … … 300 292 c ~~~~~~~~~~~~~~~~~~~~~~~ 301 293 DO isoil=1,nsoil 302 tsoil(isoil)=tsurf(1)294 ftsoil(1,isoil)=ftsol(1) 303 295 ENDDO 304 296 … … 323 315 zval(1) = 0. 324 316 317 c Initialisation albedo 318 c ---------------------- 319 320 falbe(1)=0.1 321 325 322 c Ecriture de "startphy.nc" 326 323 c ------------------------- … … 330 327 331 328 solsw(1) = 0. 332 sollw down(1)= 0.333 dlw(1) = 0.329 sollw(1) = 0. 330 fder(1) = 0. 334 331 radsol(1) = 0. 335 332 … … 337 334 soil_model = .true. 338 335 339 call phyredem("startphy.nc ", 340 . lati,long, 341 . tsurf,tsoil,albedo, 342 . solsw,sollwdown,dlw,radsol, 343 . zmea, zstd, zsig, zgam, zthe, zpic, zval, 344 . temp) 336 call phyredem("startphy.nc") 337 338 c deallocation des variables phyredem 339 c ----------------------------------- 340 call phys_state_var_end 345 341 346 342 c======================================================================= … … 457 453 458 454 OPEN(11,file='profile.new') 459 write (11,*) tsurf460 455 DO ilayer=1,nlayer 461 456 write (11,*) zlay(ilayer),temp(ilayer),tmp1(ilayer)
Note: See TracChangeset
for help on using the changeset viewer.