Ignore:
Timestamp:
Aug 3, 2024, 2:56:58 PM (6 months ago)
Author:
abarral
Message:

Put .h into modules

Location:
LMDZ6/branches/Amaury_dev/libf/dyn3d_common
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90

    r5159 r5160  
    4141    !----------------------------------------------------------------------
    4242
    43     print *, "Call sequence information: fxhyp"
     43    PRINT *, "Call sequence information: fxhyp"
    4444
    4545    test_iim: if (iim==1) THEN
    … …  
    126126        END DO
    127127
    128         print *, "ffdx = ", ffdx
     128        PRINT *, "ffdx = ", ffdx
    129129        beta = (grossismx * ffdx - pi_d) / (ffdx - pi_d)
    130         print *, "beta = ", beta
     130        PRINT *, "beta = ", beta
    131131
    132132        IF (2. * beta - grossismx <= 0.) THEN
    133           print *, 'Bad choice of grossismx, taux, dzoomx.'
    134           print *, 'Decrease dzoomx or grossismx.'
     133          PRINT *, 'Bad choice of grossismx, taux, dzoomx.'
     134          PRINT *, 'Decrease dzoomx or grossismx.'
    135135          STOP 1
    136136        END IF
    … …  
    192192
    193193          IF (rlonm025(is2) < - pi) THEN
    194             print *, 'Rlonm025 plus petit que - pi !'
     194            PRINT *, 'Rlonm025 plus petit que - pi !'
    195195            STOP 1
    196196          end if
    … …  
    203203
    204204          IF (rlonm025(is2) > pi) THEN
    205             print *, 'Rlonm025 plus grand que pi !'
     205            PRINT *, 'Rlonm025 plus grand que pi !'
    206206            STOP 1
    207207          end if
    … …  
    215215
    216216      forall (i = 1:iim) d_rlonv(i) = rlonv(i + 1) - rlonv(i)
    217       print *, "Minimum longitude step:", MINval(d_rlonv) * 180. / pi, &
     217      PRINT *, "Minimum longitude step:", MINval(d_rlonv) * 180. / pi, &
    218218              "degrees"
    219       print *, "Maximum longitude step:", MAXval(d_rlonv) * 180. / pi, &
     219      PRINT *, "Maximum longitude step:", MAXval(d_rlonv) * 180. / pi, &
    220220              "degrees"
    221221
    … …  
    223223      DO i = 1, iim + 1
    224224        IF (rlonp025(i) < rlonv(i)) THEN
    225           print *, 'rlonp025(', i, ') = ', rlonp025(i)
    226           print *, "< rlonv(", i, ") = ", rlonv(i)
     225          PRINT *, 'rlonp025(', i, ') = ', rlonp025(i)
     226          PRINT *, "< rlonv(", i, ") = ", rlonv(i)
    227227          STOP 1
    228228        END IF
    229229
    230230        IF (rlonv(i) < rlonm025(i)) THEN
    231           print *, 'rlonv(', i, ') = ', rlonv(i)
    232           print *, "< rlonm025(", i, ") = ", rlonm025(i)
     231          PRINT *, 'rlonv(', i, ') = ', rlonv(i)
     232          PRINT *, "< rlonm025(", i, ") = ", rlonm025(i)
    233233          STOP 1
    234234        END IF
    235235
    236236        IF (rlonp025(i) > rlonu(i)) THEN
    237           print *, 'rlonp025(', i, ') = ', rlonp025(i)
    238           print *, "> rlonu(", i, ") = ", rlonu(i)
     237          PRINT *, 'rlonp025(', i, ') = ', rlonp025(i)
     238          PRINT *, "> rlonu(", i, ") = ", rlonu(i)
    239239          STOP 1
    240240        END IF

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fyhyp_m.F90

    r5159 r5160  
    5959    !-------------------------------------------------------------------
    6060
    61     print *, "Call sequence information: fyhyp"
     61    PRINT *, "Call sequence information: fyhyp"
    6262
    6363    pi = 2.*asin(1.)
    … …  
    6767    y0 = clat*pi/180.
    6868    dzoom = dzoomy*pi
    69     print *, 'yzoom(rad), grossismy, tauy, dzoom (rad):'
    70     print *, y0, grossismy, tauy, dzoom
     69    PRINT *, 'yzoom(rad), grossismy, tauy, dzoom (rad):'
     70    PRINT *, y0, grossismy, tauy, dzoom
    7171
    7272    DO i = 0, nmax2
    … …  
    132132
    133133    IF (2. * beta - grossismy <= 0.) THEN
    134        print *, 'Attention ! La valeur beta calculee dans la routine fyhyp ' &
     134       PRINT *, 'Attention ! La valeur beta calculee dans la routine fyhyp ' &
    135135            // 'est mauvaise. Modifier les valeurs de grossismy, tauy ou ' &
    136136            // 'dzoomy et relancer.'
    … …  
    209209          END DO
    210210          IF (abs(yi-yo1) > epsilon) THEN
    211              print *, 'Pas de solution.', j, ylon2
     211             PRINT *, 'Pas de solution.', j, ylon2
    212212             STOP 1
    213213          end if
    … …  
    220220       DO j = 1, jlat - 1
    221221          IF (yvrai(j + 1)<yvrai(j)) THEN
    222              print *, 'Problème avec rlat(', j + 1, ') plus petit que rlat(', &
     222             PRINT *, 'Problème avec rlat(', j + 1, ') plus petit que rlat(', &
    223223                  j, ')'
    224224             STOP 1
    … …  
    226226       END DO
    227227
    228        print *, 'Reorganisation des latitudes pour avoir entre - pi/2 et pi/2'
     228       PRINT *, 'Reorganisation des latitudes pour avoir entre - pi/2 et pi/2'
    229229
    230230       IF (ik==1) THEN
    … …  
    299299    DO j = 1, jjm
    300300       IF (rlatu1(j) <= rlatu2(j)) THEN
    301           print *, 'Attention ! rlatu1 < rlatu2 ', rlatu1(j), rlatu2(j), j
     301          PRINT *, 'Attention ! rlatu1 < rlatu2 ', rlatu1(j), rlatu2(j), j
    302302          STOP 13
    303303       ENDIF
    304304
    305305       IF (rlatu2(j) <= rlatu(j+1)) THEN
    306           print *, 'Attention ! rlatu2 < rlatup1 ', rlatu2(j), rlatu(j+1), j
     306          PRINT *, 'Attention ! rlatu2 < rlatup1 ', rlatu2(j), rlatu(j+1), j
    307307          STOP 14
    308308       ENDIF
    309309
    310310       IF (rlatu(j) <= rlatu1(j)) THEN
    311           print *, ' Attention ! rlatu < rlatu1 ', rlatu(j), rlatu1(j), j
     311          PRINT *, ' Attention ! rlatu < rlatu1 ', rlatu(j), rlatu1(j), j
    312312          STOP 15
    313313       ENDIF
    314314
    315315       IF (rlatv(j) <= rlatu2(j)) THEN
    316           print *, ' Attention ! rlatv < rlatu2 ', rlatv(j), rlatu2(j), j
     316          PRINT *, ' Attention ! rlatv < rlatu2 ', rlatv(j), rlatu2(j), j
    317317          STOP 16
    318318       ENDIF
    319319
    320320       IF (rlatv(j) >= rlatu1(j)) THEN
    321           print *, ' Attention ! rlatv > rlatu1 ', rlatv(j), rlatu1(j), j
     321          PRINT *, ' Attention ! rlatv > rlatu1 ', rlatv(j), rlatu1(j), j
    322322          STOP 17
    323323       ENDIF
    324324
    325325       IF (rlatv(j) >= rlatu(j)) THEN
    326           print *, ' Attention ! rlatv > rlatu ', rlatv(j), rlatu(j), j
     326          PRINT *, ' Attention ! rlatv > rlatu ', rlatv(j), rlatu(j), j
    327327          STOP 18
    328328       ENDIF
    329329    ENDDO
    330330
    331     print *, 'Latitudes'
     331    PRINT *, 'Latitudes'
    332332    print 3, champmin, champmax
    333333

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90

    r5159 r5160  
    411411
    412412    IF (ABS(abs(xo(1)) - 90) < 0.001 .OR. ABS(abs(xo(nmax)) - 90) < 0.001) THEN
    413       print *, "ord_coord"
     413      PRINT *, "ord_coord"
    414414      PRINT *, '"xi" should contain the latitudes of the boundaries of ' &
    415415              // 'grid cells, not the centers of grid cells.'

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90

    r5159 r5160  
    6767       IF (ABS(xvrai(i) - xo1) > my_eps) THEN
    6868          ! iter == 300
    69           print *, 'Pas de solution.'
    70           print *, i, xfi
     69          PRINT *, 'Pas de solution.'
     70          PRINT *, i, xfi
    7171          STOP 1
    7272       end if
    … …  
    7777    DO i = 1, iim -1
    7878       IF (xvrai(i + 1) < xvrai(i)) THEN
    79           print *, 'xvrai(', i + 1, ') < xvrai(', i, ')'
     79          PRINT *, 'xvrai(', i + 1, ') < xvrai(', i, ')'
    8080          STOP 1
    8181       END IF

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.f90

    r5159 r5160  
    905905        ( wk2(i,k  )*(c1*(A1 - A2)+A2*flux(i,k-1)) - &
    906906        wk2(i,k-1)*A1*flux(i,k)  )
    907    ! print *,'AL1',i,k, AL(i,k)
     907   ! PRINT *,'AL1',i,k, AL(i,k)
    908908  END DO
    909909  END DO
    … …  
    911911  DO i=1,IMR*NLAYM1
    912912  AR(i,1) = AL(i,2)
    913    ! print *,'AR1',i,AR(i,1)
     913   ! PRINT *,'AR1',i,AR(i,1)
    914914  END DO
    915915
    916916  DO i=1,IMR*NLAY
    917917  A6(i,1) = 3.*(wk1(i,1)+wk1(i,1) - (AL(i,1)+AR(i,1)))
    918    ! print *,'A61',i,A6(i,1)
     918   ! PRINT *,'A61',i,A6(i,1)
    919919  END DO
    920920
    … …  
    937937    flux(i,2) = AR(i,1)+0.5*CM*(AL(i,1)-AR(i,1)+A6(i,1)*(1.-R23*CM))
    938938  else
    939      ! print *,'test2-0',i,j,wz2(i,1),wk2(i,2)
     939     ! PRINT *,'test2-0',i,j,wz2(i,1),wk2(i,2)
    940940    CP= wz2(i,1) / wk2(i,2)
    941      ! print *,'testCP',CP
     941     ! PRINT *,'testCP',CP
    942942    flux(i,2) = AL(i,2)+0.5*CP*(AL(i,2)-AR(i,2)-A6(i,2)*(1.+R23*CP))
    943      ! print *,'test2',i, AL(i,2),AR(i,2),A6(i,2),R23
     943     ! PRINT *,'test2',i, AL(i,2),AR(i,2),A6(i,2),R23
    944944  ENDIF
    945945  END DO

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ugeostr.F90

    r5159 r5160  
    6767  end DO
    6868
    69   print *, 301
     69  PRINT *, 301
    7070
    7171END SUBROUTINE  ugeostr
Note: See TracChangeset for help on using the changeset viewer.