Changeset 5258 for LMDZ6/trunk/libf
- Timestamp:
- Oct 22, 2024, 4:17:24 PM (6 weeks ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 13 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/advtrac.F90
r4143 r5258 1 1 ! $Id$ 2 2 3 #define DEBUG_IO4 #undef DEBUG_IO5 3 SUBROUTINE advtrac(pbaru, pbarv, p, masse,q,iapptrac,teta, flxw, pk) 6 4 ! Auteur : F. Hourdin … … 14 12 USE control_mod, ONLY: iapp_tracvl, day_step 15 13 USE comconst_mod, ONLY: dtvr 16 14 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 15 USE strings_mod, ONLY: int2str 17 16 IMPLICIT NONE 18 17 ! … … 202 201 CALL massbar(massem,massebx,masseby) 203 202 204 #ifdef DEBUG_IO 203 IF (CPPKEY_DEBUGIO) THEN 205 204 CALL WriteField_u('massem',massem) 206 205 CALL WriteField_u('wg',wg) … … 213 212 CALL WriteField_u('q_adv'//trim(int2str(iq)),q(:,:,iq)) 214 213 END DO 215 #endif 214 END IF 216 215 217 216 IF(isoCheck) WRITE(*,*) 'advtrac 227' -
LMDZ6/trunk/libf/dyn3dmem/advect_new_loc.F90
r5246 r5258 9 9 USE comconst_mod, ONLY: daysec 10 10 USE logic_mod, ONLY: conser 11 11 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 12 12 IMPLICIT NONE 13 13 !======================================================================= … … 250 250 if (pole_nord) ijb=ijb+iip1 251 251 if (pole_sud) ije=ije-iip1 252 #ifdef DEBUG_IO 252 IF (CPPKEY_DEBUGIO) THEN 253 253 CALL WriteField_u('du_bis',du) 254 #endif 254 END IF 255 255 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 256 256 DO l=1,llm … … 264 264 ENDDO 265 265 !$OMP END DO NOWAIT 266 #ifdef DEBUG_IO 266 IF (CPPKEY_DEBUGIO) THEN 267 267 CALL WriteField_u('du1',du1) 268 268 CALL WriteField_u('du2',du2) 269 269 CALL WriteField_u('du_bis',du) 270 #endif 270 END IF 271 271 ijb=ij_begin 272 272 ije=ij_end -
LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.F90
r4469 r5258 1 1 2 #define DEBUG_IO3 #undef DEBUG_IO4 2 SUBROUTINE advtrac_loc(pbarug, pbarvg, wg, p, massem, q, teta, pk) 5 3 ! Auteur : F. Hourdin … … 21 19 USE times 22 20 USE advtrac_mod, ONLY: finmasse 23 21 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 22 USE strings_mod, ONLY: int2str 24 23 IMPLICIT NONE 25 24 ! … … 109 108 !ym ----> Reverifier lors de la parallelisation des autres schemas 110 109 111 #ifdef DEBUG_IO 110 IF (CPPKEY_DEBUGIO) THEN 112 111 CALL WriteField_u('massem',massem) 113 112 CALL WriteField_u('wg',wg) … … 120 119 CALL WriteField_u('q_adv'//trim(int2str(iq)),q(:,:,iq)) 121 120 END DO 122 #endif 121 END IF 123 122 124 123 ! … … 132 131 CALL vlspltgen_loc(q, 2., massem, wg, pbarug, pbarvg, dtvr, p, pk, teta ) 133 132 134 #ifdef DEBUG_IO 133 IF (CPPKEY_DEBUGIO) THEN 135 134 DO iq = 1, nqtot 136 135 CALL WriteField_u('q_adv'//trim(int2str(iq)),q(:,:,iq)) 137 136 END DO 138 #endif 137 END IF 139 138 140 139 GOTO 1234 -
LMDZ6/trunk/libf/dyn3dmem/caladvtrac_loc.F90
r5246 r5258 16 16 USE Vampir 17 17 USE write_field_loc 18 !18 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 19 19 IMPLICIT NONE 20 20 ! … … 131 131 !$OMP ENDDO NOWAIT 132 132 133 #ifdef DEBUG_IO 133 IF (CPPKEY_DEBUGIO) THEN 134 134 CALL WriteField_u('pbarug1',pbarug) 135 135 CALL WriteField_v('pbarvg1',pbarvg) 136 136 CALL WriteField_u('wg1',wg) 137 #endif 137 END IF 138 138 139 139 !$OMP BARRIER … … 178 178 179 179 180 #ifdef DEBUG_IO 180 IF (CPPKEY_DEBUGIO) THEN 181 181 CALL WriteField_u('pbarug1',pbarug_adv) 182 182 CALL WriteField_v('pbarvg1',pbarvg_adv) 183 183 CALL WriteField_u('wg1',wg_adv) 184 #endif 184 END IF 185 185 ! !write(*,*) 'caladvtrac 185' 186 186 CALL advtrac_loc( pbarug_adv,pbarvg_adv,wg_adv, & -
LMDZ6/trunk/libf/dyn3dmem/caldyn_loc.F90
r5246 r5258 1 !2 ! $Id: $3 !4 #undef DEBUG_IO5 !#define DEBUG_IO6 7 1 SUBROUTINE caldyn_loc & 8 2 (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis , & … … 13 7 vorpot, ecin, bern, massebxy, convm 14 8 USE comvert_mod, ONLY: ap, bp 15 9 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 16 10 IMPLICIT NONE 17 11 … … 90 84 CALL convmas2_loc ( convm ) 91 85 !$OMP BARRIER 92 #ifdef DEBUG_IO 86 IF (CPPKEY_DEBUGIO) THEN 93 87 call WriteField_u('ucont',ucont) 94 88 call WriteField_v('vcont',vcont) … … 102 96 call WriteField_u('dteta',dteta) 103 97 call WriteField_u('convm',convm) 104 #endif 98 END IF 105 99 106 100 !$OMP BARRIER … … 122 116 CALL dudv1_loc ( vorpot , pbaru , pbarv , du , dv ) 123 117 124 #ifdef DEBUG_IO 118 IF (CPPKEY_DEBUGIO) THEN 125 119 call WriteField_u('w',w) 126 120 call WriteField_v('vorpot',vorpot) 127 121 call WriteField_u('du',du) 128 122 call WriteField_v('dv',dv) 129 #endif 123 END IF 130 124 131 125 ! ! compute kinetic energy ecin() … … 136 130 CALL dudv2_loc ( teta , pkf , bern , du , dv ) 137 131 138 #ifdef DEBUG_IO 132 IF (CPPKEY_DEBUGIO) THEN 139 133 call WriteField_u('ecin',ecin) 140 134 call WriteField_u('bern',bern) … … 142 136 call WriteField_v('dv',dv) 143 137 call WriteField_u('pkf',pkf) 144 #endif 138 END IF 145 139 146 140 ijb=ij_begin-iip1 -
LMDZ6/trunk/libf/dyn3dmem/call_calfis_mod.F90
r5250 r5258 88 88 USE comvert_mod, ONLY: ap, bp, pressure_exner 89 89 USE temps_mod, ONLY: day_ini, day_ref, jd_ref, jh_ref, start_time 90 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS 91 90 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS, CPPKEY_DEBUGIO 91 USE strings_mod, ONLY: int2str 92 92 93 93 IMPLICIT NONE … … 118 118 119 119 120 #ifdef DEBUG_IO 120 IF (CPPKEY_DEBUGIO) THEN 121 121 CALL WriteField_u('ucovfi',ucov) 122 122 CALL WriteField_v('vcovfi',vcov) … … 127 127 CALL WriteField_u('qfi'//trim(int2str(iq)),q(:,:,iq)) 128 128 ENDDO 129 #endif 129 END IF 130 130 131 131 ! … … 169 169 !$OMP END MASTER 170 170 171 #ifdef DEBUG_IO 171 IF (CPPKEY_DEBUGIO) THEN 172 172 CALL WriteField_u('ucovfi',ucov) 173 173 CALL WriteField_v('vcovfi',vcov) … … 175 175 CALL WriteField_u('pfi',p) 176 176 CALL WriteField_u('pkfi',pk) 177 #endif 177 END IF 178 178 179 179 CALL SetTag(Request_physic,800) … … 219 219 220 220 221 #ifdef DEBUG_IO 221 IF (CPPKEY_DEBUGIO) THEN 222 222 CALL WriteField_u('ucovfi',ucov) 223 223 CALL WriteField_v('vcovfi',vcov) … … 228 228 CALL WriteField_u('qfi'//trim(int2str(iq)),q(:,:,iq)) 229 229 ENDDO 230 #endif 230 END IF 231 231 232 232 !$OMP BARRIER … … 302 302 303 303 304 #ifdef DEBUG_IO 304 IF (CPPKEY_DEBUGIO) THEN 305 305 CALL WriteField_u('dufi',dufi) 306 CALL WriteField_v('dvfi',dvfi) 306 CALL WriteField_v('dvfi',dvfi) 307 307 CALL WriteField_u('dtetafi',dtetafi) 308 308 CALL WriteField_u('dpfi',dpfi) … … 310 310 CALL WriteField_u('dqfi'//trim(int2str(iq)),dqfi(:,:,iq)) 311 311 ENDDO 312 #endif 312 END IF 313 313 314 314 !$OMP BARRIER … … 316 316 ! ajout des tendances physiques: 317 317 ! ------------------------------ 318 #ifdef DEBUG_IO 318 IF (CPPKEY_DEBUGIO) THEN 319 319 CALL WriteField_u('ucovfi',ucov) 320 320 CALL WriteField_v('vcovfi',vcov) … … 324 324 CALL WriteField_u('qfi'//trim(int2str(iq)),q(:,:,iq)) 325 325 ENDDO 326 #endif 327 328 #ifdef DEBUG_IO 326 END IF 327 328 IF (CPPKEY_DEBUGIO) THEN 329 329 CALL WriteField_u('ucovfi',ucov) 330 330 CALL WriteField_v('vcovfi',vcov) … … 334 334 CALL WriteField_u('qfi'//trim(int2str(iq)),q(:,:,iq)) 335 335 ENDDO 336 #endif 336 END IF 337 337 338 338 CALL addfi_loc( dtphys, leapf, forward , & … … 351 351 !$OMP BARRIER 352 352 353 #ifdef DEBUG_IO 353 IF (CPPKEY_DEBUGIO) THEN 354 354 CALL WriteField_u('ucovfi',ucov) 355 355 CALL WriteField_v('vcovfi',vcov) … … 359 359 CALL WriteField_u('qfi'//trim(int2str(iq)),q(:,:,iq)) 360 360 ENDDO 361 #endif 361 END IF 362 362 363 363 IF (ok_strato) THEN … … 398 398 ENDIF 399 399 400 #ifdef DEBUG_IO 400 IF (CPPKEY_DEBUGIO) THEN 401 401 CALL WriteField_u('ucovfi',ucov_dyn) 402 402 CALL WriteField_v('vcovfi',vcov_dyn) … … 406 406 CALL WriteField_u('qfi'//trim(int2str(iq)),q_dyn(:,:,iq)) 407 407 ENDDO 408 #endif 408 END IF 409 409 410 410 -
LMDZ6/trunk/libf/dyn3dmem/call_dissip_mod.F90
r3435 r5258 87 87 USE vampir 88 88 USE write_field_loc 89 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 89 90 IMPLICIT NONE 90 91 INCLUDE 'comgeom.h' … … 144 145 CALL dissip_loc(vcov,ucov,teta,p,dvdis,dudis,dtetadis) 145 146 146 #ifdef DEBUG_IO 147 IF (CPPKEY_DEBUGIO) THEN 147 148 CALL WriteField_u('dudis',dudis) 148 149 CALL WriteField_v('dvdis',dvdis) 149 150 CALL WriteField_u('dtetadis',dtetadis) 150 #endif 151 END IF 151 152 152 153 ! CALL FTRACE_REGION_END("dissip") -
LMDZ6/trunk/libf/dyn3dmem/dissip_loc.F90
r5246 r5258 8 8 USE dissip_mod, ONLY: dissip_allocate 9 9 USE comconst_mod, ONLY: dtdiss 10 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 10 11 IMPLICIT NONE 11 12 … … 116 117 ENDIF 117 118 118 #ifdef DEBUG_IO 119 IF (CPPKEY_DEBUGIO) THEN 119 120 call WriteField_u('gdx',gdx) 120 121 call WriteField_v('gdy',gdy) 121 #endif 122 END IF 122 123 123 124 ijb=ij_begin … … 161 162 ENDIF 162 163 163 #ifdef DEBUG_IO 164 IF (CPPKEY_DEBUGIO) THEN 164 165 call WriteField_u('grx',grx) 165 166 call WriteField_v('gry',gry) 166 #endif 167 END IF 167 168 168 169 … … 215 216 ENDIF 216 217 217 #ifdef DEBUG_IO 218 IF (CPPKEY_DEBUGIO) THEN 218 219 call WriteField_u('gdx',gdx) 219 #endif 220 END IF 220 221 221 222 -
LMDZ6/trunk/libf/dyn3dmem/exner_hyb_loc_m.F90
r2600 r5258 37 37 USE comconst_mod, ONLY: cpp, kappa, r, jmp1 38 38 USE comvert_mod, ONLY: preff 39 39 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 40 40 IMPLICIT NONE 41 41 ! … … 185 185 jjb=jj_begin 186 186 jje=jj_end 187 #ifdef DEBUG_IO 187 IF (CPPKEY_DEBUGIO) THEN 188 188 call WriteField_u('pkf',pkf) 189 #endif 189 END IF 190 190 CALL filtreg_p ( pkf,jjb_u,jje_u,jjb,jje, jmp1, llm, & 191 191 2, 1, .TRUE., 1 ) 192 #ifdef DEBUG_IO 192 IF (CPPKEY_DEBUGIO) THEN 193 193 call WriteField_u('pkf',pkf) 194 #endif 194 END IF 195 195 end if 196 196 -
LMDZ6/trunk/libf/dyn3dmem/groupe_loc.F90
r5246 r5258 4 4 USE groupe_mod 5 5 USE comconst_mod, ONLY: ngroup 6 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 6 7 implicit none 7 8 … … 71 72 !$OMP END DO NOWAIT 72 73 73 #ifdef DEBUG_IO 74 IF (CPPKEY_DEBUGIO) THEN 74 75 CALL WriteField_v('pbarvm',reshape(pbarvm,(/ip1jm,llm/))) 75 #endif 76 END IF 76 77 call groupeun_loc(jjm,llm,jjb_v,jje_v,jjb,jje,pbarvm) 77 #ifdef DEBUG_IO 78 IF (CPPKEY_DEBUGIO) THEN 78 79 CALL WriteField_v('pbarvm',reshape(pbarvm,(/ip1jm,llm/))) 79 #endif 80 END IF 80 81 ! Champs 3D 81 82 -
LMDZ6/trunk/libf/dyn3dmem/integrd_loc.f90
r5246 r5258 15 15 USE comvert_mod, ONLY: ap, bp 16 16 USE temps_mod, ONLY: dt 17 17 USE strings_mod, ONLY: int2str 18 18 IMPLICIT NONE 19 19 -
LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.F90
r5252 r5258 1 !2 ! $Id$3 !4 !5 !6 #define DEBUG_IO7 #undef DEBUG_IO8 9 10 1 SUBROUTINE leapfrog_loc(ucov0,vcov0,teta0,ps0, & 11 2 masse0,phis0,q0,time_0) … … 46 37 xios_set_current_context, & 47 38 using_xios 48 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_REPROBUS 39 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_REPROBUS, CPPKEY_DEBUGIO 40 USE strings_mod, ONLY: int2str 49 41 50 42 IMPLICIT NONE … … 696 688 !$OMP END MASTER 697 689 698 #ifdef DEBUG_IO 690 IF (CPPKEY_DEBUGIO) THEN 699 691 call WriteField_u('du',du) 700 692 call WriteField_v('dv',dv) … … 707 699 call WriteField_u('masse',masse) 708 700 call WriteField_u('pk',pk) 709 #endif 701 END IF 710 702 !----------------------------------------------------------------------- 711 703 ! calcul des tendances advection des traceurs (dont l'humidite) … … 747 739 call VTb(VTintegre) 748 740 !$OMP END MASTER 749 #ifdef DEBUG_IO 741 IF (CPPKEY_DEBUGIO) THEN 750 742 if (true_itau>20) then 751 743 call WriteField_u('ucovm1',ucovm1) … … 758 750 call WriteField_u('ps_int',ps) 759 751 endif 760 #endif 752 END IF 761 753 !$OMP BARRIER 762 754 ! CALL FTRACE_REGION_BEGIN("integrd") … … 775 767 ! CALL FTRACE_REGION_END("integrd") 776 768 !$OMP BARRIER 777 #ifdef DEBUG_IO 769 IF (CPPKEY_DEBUGIO) THEN 778 770 call WriteField_u('ucovm1',ucovm1) 779 771 call WriteField_v('vcovm1',vcovm1) … … 784 776 call WriteField_u('teta_int',teta) 785 777 call WriteField_u('ps_int',ps) 786 #endif 778 END IF 787 779 788 780 call check_isotopes(q,ijb_u,ije_u,'leapfrog 775') -
LMDZ6/trunk/libf/dyn3dmem/vlspltgen_loc.F90
r5246 r5258 33 33 USE comconst_mod, ONLY: cpp 34 34 USE logic_mod, ONLY: adv_qsat_liq 35 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO 35 36 IMPLICIT NONE 36 37 … … 187 188 ENDDO 188 189 189 #ifdef DEBUG_IO190 CALL WriteField_u('mu',mu)191 CALL WriteField_v('mv',mv)192 CALL WriteField_u('mw',mw)193 CALL WriteField_u('qsat',qsat)194 #endif195 196 190 ! ! verif temporaire 197 191 ijb=ij_begin … … 204 198 IF(tracers(iq)%parent /= 'air') CYCLE 205 199 ! !write(*,*) 'vlspltgen 192: iq,iadv=',iq,tracers(iq)%iadv 206 #ifdef DEBUG_IO207 CALL WriteField_u('zq',zq(:,:,iq))208 CALL WriteField_u('zm',zm(:,:,iq))209 #endif210 200 SELECT CASE(tracers(iq)%iadv) 211 201 CASE(0); CYCLE … … 340 330 IF(tracers(iq)%parent /= 'air') CYCLE 341 331 ! !write(*,*) 'vlspltgen 321: iq=',iq 342 #ifdef DEBUG_IO343 CALL WriteField_u('zq',zq(:,:,iq))344 CALL WriteField_u('zm',zm(:,:,iq))345 #endif346 332 347 333 SELECT CASE(tracers(iq)%iadv) … … 360 346 IF(tracers(iq)%parent /= 'air') CYCLE 361 347 ! !write(*,*) 'vlspltgen 349: iq=',iq 362 #ifdef DEBUG_IO363 CALL WriteField_u('zq',zq(:,:,iq))364 CALL WriteField_u('zm',zm(:,:,iq))365 #endif366 348 SELECT CASE(tracers(iq)%iadv) 367 349 CASE(0); CYCLE … … 462 444 IF(tracers(iq)%parent /= 'air') CYCLE 463 445 ! !write(*,*) 'vlspltgen 449: iq=',iq 464 #ifdef DEBUG_IO465 CALL WriteField_u('zq',zq(:,:,iq))466 CALL WriteField_u('zm',zm(:,:,iq))467 #endif468 446 SELECT CASE(tracers(iq)%iadv) 469 447 CASE(0); CYCLE … … 481 459 IF(tracers(iq)%parent /= 'air') CYCLE 482 460 ! !write(*,*) 'vlspltgen 477: iq=',iq 483 #ifdef DEBUG_IO484 CALL WriteField_u('zq',zq(:,:,iq))485 CALL WriteField_u('zm',zm(:,:,iq))486 #endif487 461 SELECT CASE(tracers(iq)%iadv) 488 462 CASE(0); CYCLE … … 508 482 DO iq=1,nqtot 509 483 ! !write(*,*) 'vlspltgen_loc 561, iq=',iq 510 #ifdef DEBUG_IO511 CALL WriteField_u('zq',zq(:,:,iq))512 CALL WriteField_u('zm',zm(:,:,iq))513 #endif514 484 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 515 485 DO l=1,llm
Note: See TracChangeset
for help on using the changeset viewer.