- Timestamp:
- 12/02/14 19:21:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3d_common/inigeom.F
r222 r298 16 16 c 17 17 c 18 use mod_phys_lmdz_para, only : is_master 18 19 IMPLICIT NONE 19 20 c … … 160 161 c 161 162 c 162 WRITE(6,3)163 if (is_master) WRITE(6,3) 163 164 3 FORMAT( // 10x,' .... INIGEOM date du 01/06/98 ..... ', 164 165 * //5x,' Calcul des elongations cu et cv comme sommes des 4 ' / … … 183 184 ENDIF 184 185 185 WRITE(6,*) ' gamdi_gd ',gamdi_gdiv,gamdi_grot,gamdi_h,coefdis,186 * nitergdiv,nitergrot,niterh186 if (is_master) WRITE(6,*) ' gamdi_gd ',gamdi_gdiv,gamdi_grot, 187 * gamdi_h,coefdis,nitergdiv,nitergrot,niterh 187 188 c 188 189 pi = 2.* ASIN(1.) 189 190 c 190 WRITE(6,990)191 if (is_master) WRITE(6,990) 191 192 192 193 c ---------------------------------------------------------------- … … 197 198 IF( ysinus ) THEN 198 199 c 199 WRITE(6,*) ' *** Inigeom , Y = Sinus ( Latitude ) *** ' 200 if (is_master) WRITE(6,*) ' *** Inigeom , Y = Sinus ', 201 * '( Latitude ) *** ' 200 202 c 201 203 c .... utilisation de f(x,y ) avec y = sinus de la latitude ..... … … 207 209 ELSE 208 210 c 209 WRITE(6,*) '*** Inigeom , Y = Latitude , der. sinusoid . ***' 211 if (is_master) WRITE(6,*) '*** Inigeom , Y = Latitude ,', 212 * ' der. sinusoid . ***' 210 213 211 214 c .... utilisation de f(x,y) a tangente sinusoidale , y etant la latit. ... … … 262 265 c ..................................................................... 263 266 264 WRITE(6,*)'*** Inigeom , Y = Latitude , der.tg. hyperbolique ***' 267 if (is_master) WRITE(6,*)'*** Inigeom , Y = Latitude ,', 268 * ' der.tg. hyperbolique ***' 265 269 266 270 CALL fxyhyper( clat, grossismy, dzoomy, tauy , … … 656 660 c----------------------------------------------------------------------- 657 661 c 662 if (is_master) then 658 663 WRITE(6,*) ' *** Coordonnees de la grille *** ' 659 664 WRITE(6,995) … … 661 666 WRITE(6,*) ' LONGITUDES aux pts. V ( degres ) ' 662 667 WRITE(6,995) 668 end if 663 669 DO i=1,iip1 664 670 rlonvv(i) = rlonv(i)*180./pi 665 671 ENDDO 666 WRITE(6,400) rlonvv 667 c 668 WRITE(6,995) 669 WRITE(6,*) ' LATITUDES aux pts. V ( degres ) ' 670 WRITE(6,995) 672 if (is_master) WRITE(6,400) rlonvv 673 c 674 if (is_master) WRITE(6,995) 675 if (is_master) WRITE(6,*) ' LATITUDES aux pts. V ', 676 * '( degres ) ' 677 if (is_master) WRITE(6,995) 671 678 DO i=1,jjm 672 679 rlatuu(i)=rlatv(i)*180./pi 673 680 ENDDO 674 WRITE(6,400) (rlatuu(i),i=1,jjm)681 if (is_master) WRITE(6,400) (rlatuu(i),i=1,jjm) 675 682 c 676 683 DO i=1,iip1 677 684 rlonvv(i)=rlonu(i)*180./pi 678 685 ENDDO 686 if (is_master) then 679 687 WRITE(6,995) 680 688 WRITE(6,*) ' LONGITUDES aux pts. U ( degres ) ' … … 685 693 WRITE(6,*) ' LATITUDES aux pts. U ( degres ) ' 686 694 WRITE(6,995) 695 end if 696 if (is_master) WRITE(6,995) 687 697 DO i=1,jjp1 688 698 rlatuu(i)=rlatu(i)*180./pi 689 699 ENDDO 690 WRITE(6,400) (rlatuu(i),i=1,jjp1)691 WRITE(6,995)700 if (is_master) WRITE(6,400) (rlatuu(i),i=1,jjp1) 701 if (is_master) WRITE(6,995) 692 702 c 693 703 444 format(f10.3,f6.0)
Note: See TracChangeset
for help on using the changeset viewer.