Changeset 2088 for LMDZ5/trunk/libf/dyn3d_common
- Timestamp:
- Jul 11, 2014, 2:52:53 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d_common/grid_atob.F
r1944 r2088 52 52 REAL zzmin 53 53 #endif 54 include "iniprint.h" 54 55 c 55 56 IF (imar.GT.2200 .OR. jmar.GT.1100) THEN … … 118 119 sortie(i,j) = sortie(i,j) / number(i,j) 119 120 ELSE 120 PRINT*, 'probleme,i,j=', i,j121 if (prt_level >= 1) PRINT*, 'probleme,i,j=', i,j 121 122 ccc CALL ABORT_GCM("", "", 1) 122 123 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) … … 135 136 j_proche = (ij_proche-1)/imdep + 1 136 137 i_proche = ij_proche - (j_proche-1)*imdep 137 PRINT*, "solution:", ij_proche, i_proche, j_proche 138 if (prt_level >= 1) PRINT*, "solution:", ij_proche, i_proche, 139 $ j_proche 138 140 sortie(i,j) = entree(i_proche,j_proche) 139 141 ENDIF … … 449 451 REAL zzmin 450 452 #endif 453 include "iniprint.h" 451 454 c 452 455 IF (imar.GT.400 .OR. jmar.GT.400) THEN … … 512 515 sortie(i,j) = EXP(sortie(i,j)) 513 516 ELSE 514 PRINT*, 'probleme,i,j=', i,j517 if (prt_level >= 1) PRINT*, 'probleme,i,j=', i,j 515 518 ccc CALL ABORT_GCM("", "", 1) 516 519 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) … … 529 532 j_proche = (ij_proche-1)/imdep + 1 530 533 i_proche = ij_proche - (j_proche-1)*imdep 531 PRINT*, "solution:", ij_proche, i_proche, j_proche 534 if (prt_level >= 1) PRINT*, "solution:", ij_proche, i_proche, 535 $ j_proche 532 536 sortie(i,j) = entree(i_proche,j_proche) 533 537 ENDIF … … 574 578 REAL zzmin 575 579 #endif 580 include "iniprint.h" 576 581 c 577 582 IF (imar.GT.400 .OR. jmar.GT.400) THEN … … 641 646 ENDIF 642 647 ELSE 643 PRINT*, 'probleme,i,j=', i,j648 if (prt_level >= 1) PRINT*, 'probleme,i,j=', i,j 644 649 ccc CALL ABORT_GCM("", "", 1) 645 650 CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans) … … 658 663 j_proche = (ij_proche-1)/imdep + 1 659 664 i_proche = ij_proche - (j_proche-1)*imdep 660 PRINT*, "solution:", ij_proche, i_proche, j_proche 665 if (prt_level >= 1) PRINT*, "solution:", ij_proche, i_proche, 666 $ j_proche 661 667 IF (NINT(glace01(i_proche,j_proche)).EQ.1 ) THEN 662 668 frac_ice(i,j) = 1.0 … … 710 716 INTEGER i_proche, j_proche, ij_proche 711 717 c 718 include "iniprint.h" 719 712 720 IF (immod.GT.2200 .OR. jmmod.GT.1100) THEN 713 721 PRINT*, 'immod ou jmmod trop grand', immod, jmmod … … 874 882 rugs(i,j) = EXP(rugs(i,j)) 875 883 ELSE 876 PRINT*, 'probleme,i,j=', i,j884 if (prt_level >= 1) PRINT*, 'probleme,i,j=', i,j 877 885 ccc CALL ABORT_GCM("", "", 1) 878 886 CALL dist_sphe(xmod(i),ymod(j),xtmp,ytmp,imtmp,jmtmp,distans) … … 891 899 j_proche = (ij_proche-1)/imtmp + 1 892 900 i_proche = ij_proche - (j_proche-1)*imtmp 893 PRINT*, "solution:", ij_proche, i_proche, j_proche 901 if (prt_level >= 1) PRINT*, "solution:", ij_proche, i_proche, 902 $ j_proche 894 903 rugs(i,j) = LOG(MAX(0.001_8,cham2tmp(i_proche,j_proche))) 895 904 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.