Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.f90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.f90
r5113 r5116 86 86 INTEGER :: ierr 87 87 88 ! write(*,*) 'integrd 88: entree, nq=',nq88 !WRITE(*,*) 'integrd 88: entree, nq=',nq 89 89 !----------------------------------------------------------------------- 90 90 … … 151 151 ! & MPI_LOGICAL,MPI_LOR,COMM_LMDZ,ierr) 152 152 IF( .NOT. checksum ) THEN 153 write(lunout,*) "integrd: ps = ", ps(stop_it)154 write(lunout,*) " at node ij =", stop_it153 WRITE(lunout,*) "integrd: ps = ", ps(stop_it) 154 WRITE(lunout,*) " at node ij =", stop_it 155 155 ! since ij=j+(i-1)*jjp1 , we have 156 156 j=modulo(stop_it,jjp1) 157 157 i=1+(stop_it-j)/jjp1 158 write(lunout,*) " lon = ",rlonv(i)*180./pi, " deg", &158 WRITE(lunout,*) " lon = ",rlonv(i)*180./pi, " deg", & 159 159 " lat = ",rlatu(j)*180./pi, " deg" 160 160 CALL abort_gcm("integrd_loc", "negative surface pressure", 1) … … 163 163 !$OMP END MASTER 164 164 !$OMP BARRIER 165 ! write(*,*) 'integrd 170'165 !WRITE(*,*) 'integrd 170' 166 166 IF (.NOT. Checksum_all) THEN 167 167 CALL WriteField_v('int_vcov',vcov) … … 193 193 194 194 ! 195 ! ! write(*,*) 'integrd 200'195 ! !WRITE(*,*) 'integrd 200' 196 196 !$OMP MASTER 197 197 if (pole_nord) THEN … … 220 220 !$OMP END MASTER 221 221 !$OMP BARRIER 222 ! write(*,*) 'integrd 217'222 !WRITE(*,*) 'integrd 217' 223 223 ! 224 224 ! ... Calcul de la nouvelle masse d'air au dernier temps integre t+1 ... … … 284 284 ! 285 285 ! 286 ! ! write(*,*) 'integrd 291'286 ! !WRITE(*,*) 'integrd 291' 287 287 IF (pole_nord) THEN 288 288 … … 334 334 ije=ij_end 335 335 336 if (planet_type=="earth") then336 if (planet_type=="earth") THEN 337 337 ! Earth-specific treatment of first 2 tracers (water) 338 338 !$OMP BARRIER … … 349 349 CALL check_isotopes(q,ijb,ije,'integrd 342') 350 350 351 ! write(*,*) 'integrd 341'351 !WRITE(*,*) 'integrd 341' 352 352 CALL qminimum_loc( q, nq, deltap ) 353 ! write(*,*) 'integrd 343'353 !WRITE(*,*) 'integrd 343' 354 354 355 355 CALL check_isotopes(q,ijb,ije,'integrd 346') … … 421 421 ! ..... FIN de l'integration de q ....... 422 422 423 ! write(*,*) 'integrd 410'423 !WRITE(*,*) 'integrd 410' 424 424 425 425 !$OMP DO SCHEDULE(STATIC)
Note: See TracChangeset
for help on using the changeset viewer.