Changeset 1142
- Timestamp:
- Apr 7, 2009, 5:53:39 PM (16 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev
- Files:
-
- 2 added
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F
r1114 r1142 369 369 CALL gr_dyn_fi_p(llm,iip1,jjp1,klon,pphi,zphi) 370 370 371 c$OMP MASTER372 371 CALL gr_dyn_fi_p(1,iip1,jjp1,klon,pphis,zphis) 373 c$OMP END MASTER 372 374 373 c$OMP BARRIER 375 374 … … 895 894 c tendance sur la pression : 896 895 c ----------------------------------- 897 c$OMP MASTER898 896 CALL gr_fi_dyn_p(1,klon,iip1,jjp1,zdpsrf,pdpsfi) 899 c$OMP END MASTER900 897 c 901 898 c 62. enthalpie potentielle -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/times.F90
r1000 r1142 207 207 endif 208 208 209 ENDIF ! using_mp î209 ENDIF ! using_mp� 210 210 end subroutine allgather_timer_average 211 211 … … 222 222 end subroutine InitTime 223 223 224 function DiffTime 224 function DiffTime() 225 225 implicit none 226 226 double precision :: DiffTime … … 236 236 end function DiffTime 237 237 238 function DiffCpuTime 238 function DiffCpuTime() 239 239 implicit none 240 240 real :: DiffCpuTime -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/calltherm.F90
r1139 r1142 16 16 #include "thermcell.h" 17 17 #include "iniprint.h" 18 19 ! A inclure eventuellement dans les fichiers de configuration20 data r_aspect_thermals,l_mix_thermals/2.,30./21 data w2di_thermals/1/22 18 23 19 !IM 140508 … … 176 172 & ,tau_thermals,3) 177 173 else if (iflag_thermals.eq.11) then 178 stop 'cas non prevu dans calltherm'174 stop 'cas non prevu dans calltherm' 179 175 ! CALL thermcell_pluie(klon,klev,zdt & 180 176 ! & ,pplay,paprs,pphi,zlev & -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90
r1069 r1142 460 460 !**************************************************************************************** 461 461 ! Declarations specifiques pour le 1D. A reprendre 462 REAL :: fsens,flat 463 LOGICAL ok_flux_surf 464 DATA ok_flux_surf/.FALSE./ 465 !ym pas glop !! 466 COMMON /flux_arp/fsens,flat,ok_flux_surf 467 !$OMP THREADPRIVATE(/flux_arp/) 462 REAL,SAVE :: fsens,flat 463 LOGICAL,SAVE :: ok_flux_surf=.FALSE. 464 !$OMP THREADPRIVATE(fsens,flat,ok_flux_surf) 468 465 469 466 !**************************************************************************************** -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/thermcell.h
r1026 r1142 1 integer iflag_thermals,nsplit_thermals 2 real r_aspect_thermals,l_mix_thermals,tau_thermals 3 integer w2di_thermals,isplit 4 integer iflag_coupl,iflag_clos,iflag_wake 5 integer iflag_thermals_ed,iflag_thermals_optflux 1 integer :: iflag_thermals,nsplit_thermals 2 real,parameter :: r_aspect_thermals=2.,l_mix_thermals=30. 3 real :: tau_thermals 4 integer,parameter :: w2di_thermals=1 5 integer :: isplit 6 7 integer :: iflag_coupl,iflag_clos,iflag_wake 8 integer :: iflag_thermals_ed,iflag_thermals_optflux 6 9 7 10 common/ctherm1/iflag_thermals,nsplit_thermals 8 common/ctherm2/r_aspect_thermals,l_mix_thermals,tau_thermals 9 common/ctherm3/w2di_thermals 11 common/ctherm2/tau_thermals 10 12 common/ctherm4/iflag_coupl,iflag_clos,iflag_wake 11 13 common/ctherm5/iflag_thermals_ed,iflag_thermals_optflux 12 14 13 !$OMP THREADPRIVATE(/ctherm1/,/ctherm2/,/ctherm 3/,/ctherm4/)15 !$OMP THREADPRIVATE(/ctherm1/,/ctherm2/,/ctherm4/,/ctherm5/) -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/thermcell_flux.F90
r987 r1142 164 164 ! FH Version en cours de test; 165 165 ! par rapport a thermcell_flux, on fait une grande boucle sur "l" 166 ! et on modifie le flux avec tous les contr ôles appliques d'affilee166 ! et on modifie le flux avec tous les contr�les appliques d'affilee 167 167 ! pour la meme couche 168 168 ! Momentanement, on duplique le calcule du flux pour pouvoir comparer … … 264 264 if (entr(ig,l)<0.) then 265 265 print*,'N1 ig,l,entr',ig,l,entr(ig,l) 266 stop 'entr negatif'266 stop 'entr negatif' 267 267 endif 268 268 if (detr(ig,l).gt.fm(ig,l)) then … … 292 292 print*,'entr(ig,l)',entr(ig,l) 293 293 print*,'fm(ig,l)',fm(ig,l) 294 stop 'probleme dans thermcell flux'294 stop 'probleme dans thermcell flux' 295 295 endif 296 296 enddo … … 319 319 print*,'detr(ig,l)',detr(ig,l) 320 320 print*,'fm(ig,l)',fm(ig,l) 321 stop 'probleme dans thermcell flux'321 stop 'probleme dans thermcell flux' 322 322 endif 323 323 enddo … … 420 420 print*,'fm(ig,l+1)',fm(ig,l+1) 421 421 print*,'fm(ig,l)',fm(ig,l) 422 stop 'probleme dans thermcell_flux'422 stop 'probleme dans thermcell_flux' 423 423 endif 424 424 entr(ig,l+1)=entr(ig,l+1)-ddd -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/thermcell_flux2.F90
r1026 r1142 160 160 ! FH Version en cours de test; 161 161 ! par rapport a thermcell_flux, on fait une grande boucle sur "l" 162 ! et on modifie le flux avec tous les contr ôles appliques d'affilee162 ! et on modifie le flux avec tous les contr�les appliques d'affilee 163 163 ! pour la meme couche 164 164 ! Momentanement, on duplique le calcule du flux pour pouvoir comparer … … 256 256 if (entr(ig,l)<0.) then 257 257 print*,'N1 ig,l,entr',ig,l,entr(ig,l) 258 stop 'entr negatif'258 stop 'entr negatif' 259 259 endif 260 260 if (detr(ig,l).gt.fm(ig,l)) then … … 285 285 print*,'entr(ig,l)',entr(ig,l) 286 286 print*,'fm(ig,l)',fm(ig,l) 287 stop 'probleme dans thermcell flux'287 stop 'probleme dans thermcell flux' 288 288 endif 289 289 enddo … … 312 312 print*,'detr(ig,l)',detr(ig,l) 313 313 print*,'fm(ig,l)',fm(ig,l) 314 stop 'probleme dans thermcell flux'314 stop 'probleme dans thermcell flux' 315 315 endif 316 316 enddo … … 413 413 print*,'fm(ig,l+1)',fm(ig,l+1) 414 414 print*,'fm(ig,l)',fm(ig,l) 415 stop 'probleme dans thermcell_flux'415 stop 'probleme dans thermcell_flux' 416 416 endif 417 417 entr(ig,l+1)=entr(ig,l+1)-ddd -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/thermcell_main.F90
r1134 r1142 484 484 ! Test valable seulement en 1D mais pas genant 485 485 if (.not. (f0(1).ge.0.) ) then 486 stop 'Dans thermcell_main'486 stop 'Dans thermcell_main' 487 487 endif 488 488
Note: See TracChangeset
for help on using the changeset viewer.