- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phytracr_spl_mod.F90
r5144 r5158 1343 1343 1344 1344 !JE_dbg 1345 doi = 1, klon1345 DO i = 1, klon 1346 1346 tsol(i) = 0.0 1347 doj = 1, nbsrf1347 DO j = 1, nbsrf 1348 1348 tsol(i) = tsol(i) + ftsol(i, j) * pctsrf(i, j) 1349 1349 enddo … … 1382 1382 id_scdu = -1 1383 1383 itr = 0 1384 doiq = 1, nqtot1384 DO iq = 1, nqtot 1385 1385 IF(.NOT.tracers(iq)%isInPhysics) CYCLE 1386 1386 itr = itr + 1 … … 1941 1941 !$OMP END MASTER 1942 1942 !$OMP BARRIER 1943 doi = 1, klon1943 DO i = 1, klon 1944 1944 aux_var2(i) = pctsrf(i, is_oce) 1945 1945 enddo … … 1954 1954 !$OMP BARRIER 1955 1955 1956 doi = 1, klon1956 DO i = 1, klon 1957 1957 aux_var2(i) = pctsrf(i, is_sic) 1958 1958 enddo … … 1967 1967 !$OMP BARRIER 1968 1968 1969 doi = 1, klon1969 DO i = 1, klon 1970 1970 aux_var2(i) = pctsrf(i, is_ter) 1971 1971 enddo … … 1980 1980 !$OMP BARRIER 1981 1981 1982 doi = 1, klon1982 DO i = 1, klon 1983 1983 aux_var2(i) = pctsrf(i, is_lic) 1984 1984 enddo … … 2066 2066 ! . MAXVAL(tr_seri(:,:,3)) 2067 2067 #ifdef IOPHYS_DUST 2068 doitr=1,nbtr2068 DO itr=1,nbtr 2069 2069 WRITE(str2,'(i2.2)') itr 2070 2070 CALL iophys_ecrit('sav'//str2,1,'SOURCE','',source_tr(:,itr)) 2071 2071 CALL iophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,itr)) 2072 2072 enddo 2073 doitr=1,nbtr2073 DO itr=1,nbtr 2074 2074 WRITE(str2,'(i2.2)') itr 2075 2075 CALL iophys_ecrit('TRB'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2099 2099 2100 2100 #ifdef IOPHYS_DUST 2101 doitr=1,nbtr2101 DO itr=1,nbtr 2102 2102 WRITE(str2,'(i2.2)') itr 2103 2103 CALL iophys_ecrit('sap'//str2,1,'SOURCE','',source_tr(:,itr)) … … 2196 2196 #ifdef IOPHYS_DUST 2197 2197 2198 doitr=1,nbtr2198 DO itr=1,nbtr 2199 2199 WRITE(str2,'(i2.2)') itr 2200 2200 CALL iophys_ecrit('tpr'//str2,1,'SOURCE','',source_tr(:,itr)) … … 2240 2240 2241 2241 #ifdef IOPHYS_DUST 2242 doitr=1,nbtr2242 DO itr=1,nbtr 2243 2243 WRITE(str2,'(i2.2)') itr 2244 2244 CALL iophys_ecrit('t'//str2,1,'SOURCE','',source_tr(:,itr)) … … 2296 2296 2297 2297 #ifdef IOPHYS_DUST 2298 doitr=1,nbtr2298 DO itr=1,nbtr 2299 2299 WRITE(str2,'(i2.2)') itr 2300 2300 CALL iophys_ecrit('TRC'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2336 2336 2337 2337 #ifdef IOPHYS_DUST 2338 doitr=1,nbtr2338 DO itr=1,nbtr 2339 2339 WRITE(str2,'(i2.2)') itr 2340 2340 CALL iophys_ecrit('TRD'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2410 2410 CALL iophys_ecrit('yv1',1,'yv1','',yv1) 2411 2411 CALL iophys_ecrit('delp',klev,'delp','',delp) 2412 doitr=1,nbtr2412 DO itr=1,nbtr 2413 2413 WRITE(str2,'(i2.2)') itr 2414 2414 CALL iophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2602 2602 !======================================================================= 2603 2603 #ifdef IOPHYS_DUST 2604 doitr=1,nbtr2604 DO itr=1,nbtr 2605 2605 WRITE(str2,'(i2.2)') itr 2606 2606 CALL iophys_ecrit('TRF'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2689 2689 2690 2690 #ifdef IOPHYS_DUST 2691 doitr=1,nbtr2691 DO itr=1,nbtr 2692 2692 WRITE(str2,'(i2.2)') itr 2693 2693 CALL iophys_ecrit('TRG'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 2854 2854 2855 2855 #ifdef IOPHYS_DUST 2856 doitr=1,nbtr2856 DO itr=1,nbtr 2857 2857 WRITE(str2,'(i2.2)') itr 2858 2858 CALL iophys_ecrit('TRH'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 3053 3053 CALL iophys_ecrit('wdtrainM',klev,'wdtrainM','',wdtrainM) 3054 3054 3055 doitr=1,nbtr3055 DO itr=1,nbtr 3056 3056 WRITE(str2,'(i2.2)') itr 3057 3057 CALL iophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 3174 3174 !======================================================================= 3175 3175 #ifdef IOPHYS_DUST 3176 doitr=1,nbtr3176 DO itr=1,nbtr 3177 3177 WRITE(str2,'(i2.2)') itr 3178 3178 CALL iophys_ecrit('TRJ'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 3440 3440 !====================================================================== 3441 3441 #ifdef IOPHYS_DUST 3442 doitr=1,nbtr3442 DO itr=1,nbtr 3443 3443 WRITE(str2,'(i2.2)') itr 3444 3444 CALL iophys_ecrit('TRK'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) … … 4306 4306 IF (is_mpi_root .AND. is_omp_root) THEN 4307 4307 OPEN(1, file = trim(adjustl(filescaleparams)), form = 'unformatted') 4308 dok = 1, nbreg4308 DO k = 1, nbreg 4309 4309 read(1) scale_param(k) 4310 4310 enddo … … 4454 4454 ierr = nf90_open(trim(adjustl(filescaleparams)), nf90_nowrite, nid) 4455 4455 IF (ierr == nf90_noerr) THEN 4456 doi = 1, nbreg4456 DO i = 1, nbreg 4457 4457 WRITE(aux_2s, '(i2.2)') i 4458 4458 varname = trim(adjustl(paramname)) // aux_2s
Note: See TracChangeset
for help on using the changeset viewer.