Ignore:
Timestamp:
Jul 19, 2024, 7:54:50 PM (6 months ago)
Author:
abarral
Message:

convert labeled do (f77) to do .. end do

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd
Files:
39 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90

    r3834 r5086  
    140140                  rain = null_array, qcol = rnl + hf + hlb - dels)
    141141          end if
    142        end do
     142       END DO
    143143    else
    144144       tkt = 0.

  • LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90

    r5082 r5086  
    731731         JX=0
    732732       ELSE
    733          DO 10 I=1,N
     733         DO I=1,N
    734734           IF (X<XC(I)) GO TO 20
    735  10      CONTINUE
     735         END DO
    736736         IER=1
    737737 20      JX=I-1
    … …  
    756756         JX=0
    757757       ELSE
    758          DO 10 I=1,N
     758         DO I=1,N
    759759           IF (XT>X(I)) GO TO 20
    760  10      CONTINUE
     760         END DO
    761761 20      JX=I
    762762       ENDIF

  • LMDZ6/branches/Amaury_dev/libf/phylmd/concvl.F90

    r4613 r5086  
    340340! em_wght(k)=wght_th(i,k)
    341341! print*,'em_wght=',em_wght(k),wght_th(i,k)
    342 ! end do
     342! END DO
    343343! END DO
    344344

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/icarus.F

    r5082 r5086  
    364364      enddo
    365365
    366       do 12 ilev=1,nlev
     366      do ilev=1,nlev
    367367        do j=1,npoints
    368368         if (pfull(j,ilev) < 40000. .and.
    … …  
    375375           end if
    376376        enddo
    377 12    continue
    378 
    379       do 13 ilev=1,nlev
     377      END DO
     378
     379      do ilev=1,nlev
    380380        do j=1,npoints
    381381           if (at(j,ilev) > atmax(j) .and.
    382382     &              ilev  >= itrop(j)) atmax(j)=at(j,ilev)
    383383        enddo
    384 13    continue
     384      END DO
    385385
    386386      end if
    … …  
    391391              meantb(j) = 0.
    392392              meantbclr(j) = 0.
    393           end do
     393          END DO
    394394      else
    395395          do j=1,npoints
    396396              meantb(j) = output_missing_value
    397397              meantbclr(j) = output_missing_value
    398           end do
     398          END DO
    399399      end if
    400400     
    … …  
    466466 
    467467      !initialize tau and albedocld to zero
    468       do 15 ibox=1,ncol
     468      do ibox=1,ncol
    469469        do j=1,npoints
    470470            tau(j,ibox)=0.
    … …  
    474474          box_cloudy(j,ibox)=.false.
    475475        enddo
    476 15    continue
     476      END DO
    477477
    478478      !compute total cloud optical depth for each column     
    … …  
    541541        if (ncolprint /= 0)
    542542     &         write(6,*)  'ilev   pw (kg/m2)   tauwv(j)      dem_wv'
    543         do 125 ilev=1,nlev
     543        do ilev=1,nlev
    544544          do j=1,npoints
    545545               !press and dpress are dyne/cm2 = Pascals *10
    … …  
    568568               enddo
    569569             endif
    570 125     continue
     570      END DO
    571571
    572572        !initialize variables
    … …  
    740740            bb(j)=1/( exp(1307.27/skt(j)) - 1. )
    741741            !bb(j)=5.67e-8*skt(j)**4
    742           end do
     742          END DO
    743743
    744744        do ibox=1,ncol
    … …  
    751751     &         * trans_layers_above(j,ibox)
    752752           
    753           end do
    754         end do
     753          END DO
     754        END DO
    755755
    756756        !calculate mean infrared brightness temperature
    … …  
    758758          do j=1,npoints
    759759            meantb(j) = meantb(j)+1307.27/(log(1.+(1./fluxtop(j,ibox))))
    760           end do
    761         end do
     760          END DO
     761        END DO
    762762          do j=1, npoints
    763763            meantb(j) = meantb(j) / real(ncol)
    764           end do       
     764          END DO
    765765
    766766        if (ncolprint/=0) then
    … …  
    784784          write (6,'(8f7.2)') (meantb(j),ibox=1,ncolprint)
    785785     
    786           end do
     786          END DO
    787787      endif
    788788   
    … …  
    925925
    926926      !compute cloud top pressure
    927       do 30 ibox=1,ncol
     927      do ibox=1,ncol
    928928        !segregate according to optical thickness
    929929        if (top_height == 1 .or. top_height == 3) then
    … …  
    933933            nmatch(j)=0
    934934          enddo
    935           do 29 k1=1,nlev-1
     935          do k1=1,nlev-1
    936936            if (top_height_direction == 2) then
    937937              ilev = nlev - k1
    … …  
    951951             end if                         
    952952            enddo
    953 29        continue
     953      END DO
    954954
    955955          do j=1,npoints
    … …  
    992992              levmatch(j,ibox)=ilev
    993993              end if
    994             end do
    995           end do
     994            END DO
     995          END DO
    996996        end if                           
    997997         
    … …  
    10031003        enddo
    10041004
    1005 30    continue
     1005      END DO
    10061006             
    10071007!
    … …  
    10321032
    10331033      !reset frequencies
    1034       do 38 ilev=1,7
     1034      do ilev=1,7
    10351035      do 38 ilev2=1,7
    10361036        do j=1,npoints !
    … …  
    10421042        enddo
    1043104338    continue
     1044      END DO
    10441045
    10451046      !reset variables need for averaging cloud properties
    … …  
    10601061      boxarea = 1./real(ncol)
    10611062     
    1062       do 39 ibox=1,ncol
     1063      do ibox=1,ncol
    10631064        do j=1,npoints
    10641065
    … …  
    11661167                       
    11671168        enddo ! j
    1168 39    continue
     1169      END DO
    11691170     
    11701171      !compute mean cloud properties
    … …  
    12271228     &           (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev)
    12281229     &           ,(cchar(acc(ilev,ibox)+1),ilev=1,nlev)
    1229              end do
     1230             END DO
    12301231             close(9)
    12311232

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/math_lib.F90

    r5082 r5086  
    273273       exit
    274274    end if
    275   end do
     275  END DO
    276276 
    277277  if (lerror) then
    … …  
    316316    end if
    317317    ilo = ilo + 1
    318   end do
     318  END DO
    319319
    320320  ilo = max ( 2, ilo )
    … …  
    326326    end if
    327327    ihi = ihi - 1
    328   end do
     328  END DO
    329329 
    330330  ihi = min ( ihi, ntab - 1 )
    … …  
    374374    syl = x2
    375375 
    376   end do
     376  END DO
    377377 
    378378  result = sum1 &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp_modis_simulator.F90

    r3233 r5086  
    176176              opticalThickness(i, j, k) = 0.   
    177177            end if
    178           end do
    179         end do
    180       end do
     178          END DO
     179        END DO
     180      END DO
    181181
    182182      !
    … …  
    197197          do i = 1, nSunlit
    198198            if(subCols%frac_out(sunlit(i), j, k) == I_CVC) opticalThickness(i, j, k) = gridBox%dtau_c(sunlit(i), k)
    199           end do
    200         end do
    201       end do
     199          END DO
     200        END DO
     201      END DO
    202202
    203203      !
    … …  
    220220                                retrievedPhase(i, :), retrievedCloudTopPressure(i, :),      &
    221221                                retrievedTau(i, :), retrievedSize(i, :))
    222      end do
     222     END DO
    223223     
    224224      ! DJS2015: Call L3 modis simulator used by cospv2.0

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_modis_sim.F90

    r5082 r5086  
    331331        retrievedTau(i) = R_UNDEF
    332332      end if
    333     end do
     333    END DO
    334334    where((retrievedSize(:) < 0.).and.(retrievedSize(:) /= R_UNDEF)) retrievedSize(:) = 1.0e-06*re_fill
    335335
    … …  
    802802        tauMask(:, :, i) = .false.
    803803      end where
    804     end do
     804    END DO
    805805
    806806    do i = 1, numPressureHistogramBins
    … …  
    811811        pressureMask(:, :, i) = .false.
    812812      end where
    813     end do
     813    END DO
    814814   
    815815    do i = 1, numPressureHistogramBins
    … …  
    817817        Optical_Thickness_vs_Cloud_Top_Pressure(:, j, i) = &
    818818          real(count(tauMask(:, :, j) .and. pressureMask(:, :, i), dim = 2)) / real(nSubcols)
    819       end do
    820     end do
     819      END DO
     820    END DO
    821821   
    822822  end subroutine modis_L3_simulator
    … …  
    851851      end if
    852852      if(totalTau >= tauLimit) exit
    853     end do
     853    END DO
    854854    cloud_top_pressure = totalProduct/totalTau
    855855  end function cloud_top_pressure
    … …  
    877877      end if
    878878      if(totalTau >= tauLimit) exit
    879     end do
     879    END DO
    880880    weight_by_extinction = totalProduct/totalTau
    881881  end function weight_by_extinction
    … …  
    11141114    do i = 1, size(cloudIndicies)
    11151115      call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
    1116     end do
     1116    END DO
    11171117                   
    11181118    call adding_doubling(Refl(:), Trans(:), Refl_tot, Trans_tot) 
    … …  
    12921292          Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1 - Refl_cumulative(i-1) * Refl(i))
    12931293          Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1 - Refl_cumulative(i-1) * Refl(i))
    1294       end do
     1294      END DO
    12951295     
    12961296      Refl_tot = Refl_cumulative(size(Refl))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/optics_lib.F90

    r5081 r5086  
    519519
    520520!   // region from 0.045 microns to 167.0 microns - no temperature depend
    521     do i=2,nwl
    522       if(alam < wl(i)) continue
    523     enddo
    524521    x1=log(wl(i-1))
    525522    x2=log(wl(i))
    … …  
    539536    if(tk > temref(1)) tk=temref(1)
    540537    if(tk < temref(4)) tk=temref(4)
    541     do 11 i=2,4
     538    do i=2,4
    542539      if(tk>=temref(i)) go to 12
    543     11 continue
     540    END DO
    544541    12 lt1=i
    545542    lt2=i-1
    546     do 13 i=2,nwlt
     543    do i=2,nwlt
    547544      if(alam<=wlt(i)) go to 14
    548     13 continue
     545    END DO
    549546    14 x1=log(wlt(i-1))
    550547    x2=log(wlt(i))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/scops.F

    r5082 r5086  
    161161     
    162162      !loop over vertical levels
    163       DO 200 ilev = 1,nlev
     163      DO ilev = 1,nlev
    164164                                 
    165165!     Initialise threshold
    … …  
    331331          endif
    332332
    333 200   CONTINUE    !loop over nlev
     333      END DO    !loop over nlev
    334334
    335335

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/zeff.F90

    r5082 r5086  
    113113      call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
    114114        dg, xs1, xs2, dph, err)
    115     end do
     115    END DO
    116116   
    117117  else

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp.F90

    r5082 r5086  
    731731                                  modisRetrievedCloudTopPressure(i,:),                   &
    732732                                  modisRetrievedTau(i,:),modisRetrievedSize(i,:))
    733           end do
     733          END DO
    734734       endif
    735735    endif

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/icarus.F90

    r5082 r5086  
    155155                  (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev),&
    156156                  (cchar(acc(ilev,ibox)+1),ilev=1,nlev)
    157           end do
     157          END DO
    158158          close(9)
    159159
    … …  
    322322       do ibox=1,ncol
    323323          fluxtop(1:npoints,ibox) = fluxtop(1:npoints,ibox) + emsfc_lw*bb(1:npoints)*trans_layers_above(1:npoints,ibox)
    324        end do
     324       END DO
    325325
    326326       ! All Sky brightness temperature
    … …  
    445445                levmatch(1:npoints,ibox)=ilev
    446446             endwhere
    447           end do
     447          END DO
    448448       end if
    449449       where(tau(1:npoints,ibox) <= tauchk)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90

    r3358 r5086  
    209209          exit
    210210       end if
    211     end do
     211    END DO
    212212   
    213213    if (lerror) then
    … …  
    244244       end if
    245245       ilo = ilo + 1
    246     end do
     246    END DO
    247247   
    248248    ilo = max ( 2, ilo )
    … …  
    254254       end if
    255255       ihi = ihi - 1
    256     end do
     256    END DO
    257257   
    258258    ihi = min ( ihi, ntab - 1 )
    … …  
    305305       syl = x2
    306306       
    307     end do
     307    END DO
    308308
    309309    result = sum1 + ca * ( b**3 - syl**3 ) / 3._wp &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/modis_simulator.F90

    r3358 r5086  
    222222          retrievedTau(i)              = R_UNDEF
    223223       end if
    224     end do
     224    END DO
    225225    where((retrievedSize(1:nSubCols) < 0.).and.(retrievedSize(1:nSubCols) /= R_UNDEF)) &
    226226         retrievedSize(1:nSubCols) = 1.0e-06_wp*re_fill
    … …  
    455455      end if
    456456      if(totalTau >= tauLimit) exit
    457     end do
     457    END DO
    458458
    459459    if (totalTau > 0._wp) then
    … …  
    489489      end if
    490490      if(totalTau >= tauLimit) exit
    491     end do
     491    END DO
    492492
    493493    if (totalTau > 0._wp) then
    … …  
    715715    do i = 1, size(cloudIndicies)
    716716       call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
    717     end do
     717    END DO
    718718   
    719719    call adding_doubling(count(tau(1:nLevels) > 0),Refl(:), Trans(:), Refl_tot, Trans_tot) 
    … …  
    897897       Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i))
    898898       Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1._wp - Refl_cumulative(i-1) * Refl(i))
    899     end do
     899    END DO
    900900   
    901901    Refl_tot = Refl_cumulative(size(Refl))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90

    r5082 r5086  
    539539    if (alam < cutice) then
    540540       ! Region from 0.045 microns to 167.0 microns - no temperature depend
    541        do i=2,nwl
    542           if(alam < wl(i)) continue
    543        enddo
    544541       x1  = log(wl(i-1))
    545542       x2  = log(wl(i))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/quickbeam_optics.F90

    r5082 r5086  
    10091009          call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
    10101010               dg, xs1, xs2, dph, err)
    1011        end do
     1011       END DO
    10121012
    10131013    else

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp.F90

    r5082 r5086  
    886886                                  modisRetrievedCloudTopPressure(i,:),                   &
    887887                                  modisRetrievedTau(i,:),modisRetrievedSize(i,:))
    888           end do
     888          END DO
    889889       endif
    890890    endif

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/icarus.F90

    r5082 r5086  
    155155                  (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev),&
    156156                  (cchar(acc(ilev,ibox)+1),ilev=1,nlev)
    157           end do
     157          END DO
    158158          close(9)
    159159
    … …  
    322322       do ibox=1,ncol
    323323          fluxtop(1:npoints,ibox) = fluxtop(1:npoints,ibox) + emsfc_lw*bb(1:npoints)*trans_layers_above(1:npoints,ibox)
    324        end do
     324       END DO
    325325
    326326       ! All Sky brightness temperature
    … …  
    445445                levmatch(1:npoints,ibox)=ilev
    446446             endwhere
    447           end do
     447          END DO
    448448       end if
    449449       where(tau(1:npoints,ibox) <= tauchk)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/lmdz_cosp_subsample_and_optics_mod.F90

    r5082 r5086  
    365365             endif
    366366             cospIN%g_vol_cloudsat(i,:,j)=g_vol(i,j)
    367           end do
    368        end do
     367          END DO
     368       END DO
    369369       
    370370       ! Loop over all subcolumns

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/math_lib.F90

    r3491 r5086  
    209209          exit
    210210       end if
    211     end do
     211    END DO
    212212   
    213213    if (lerror) then
    … …  
    244244       end if
    245245       ilo = ilo + 1
    246     end do
     246    END DO
    247247   
    248248    ilo = max ( 2, ilo )
    … …  
    254254       end if
    255255       ihi = ihi - 1
    256     end do
     256    END DO
    257257   
    258258    ihi = min ( ihi, ntab - 1 )
    … …  
    305305       syl = x2
    306306       
    307     end do
     307    END DO
    308308
    309309    result = sum1 + ca * ( b**3 - syl**3 ) / 3._wp &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90

    r3491 r5086  
    222222          retrievedTau(i)              = R_UNDEF
    223223       end if
    224     end do
     224    END DO
    225225    where((retrievedSize(1:nSubCols) < 0.).and.(retrievedSize(1:nSubCols) /= R_UNDEF)) &
    226226         retrievedSize(1:nSubCols) = 1.0e-06_wp*re_fill
    … …  
    455455      end if
    456456      if(totalTau >= tauLimit) exit
    457     end do
     457    END DO
    458458
    459459    if (totalTau > 0._wp) then
    … …  
    489489      end if
    490490      if(totalTau >= tauLimit) exit
    491     end do
     491    END DO
    492492
    493493    if (totalTau > 0._wp) then
    … …  
    715715    do i = 1, size(cloudIndicies)
    716716       call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
    717     end do
     717    END DO
    718718   
    719719    call adding_doubling(count(tau(1:nLevels) > 0),Refl(:), Trans(:), Refl_tot, Trans_tot) 
    … …  
    897897       Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i))
    898898       Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1._wp - Refl_cumulative(i-1) * Refl(i))
    899     end do
     899    END DO
    900900   
    901901    Refl_tot = Refl_cumulative(size(Refl))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/optics_lib.F90

    r5081 r5086  
    539539    if (alam < cutice) then
    540540       ! Region from 0.045 microns to 167.0 microns - no temperature depend
    541        do i=2,nwl
    542           if(alam < wl(i)) continue
    543        enddo
    544541       x1  = log(wl(i-1))
    545542       x2  = log(wl(i))

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/quickbeam_optics.F90

    r5081 r5086  
    983983          call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
    984984               dg, xs1, xs2, dph, err)
    985        end do
     985       END DO
    986986
    987987    else

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90

    r5082 r5086  
    15481548            ! !!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
    15491549            ! !!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
    1550             ! !!!      end do
     1550            ! !!!      END DO
    15511551            elij(il, i, j) = altem
    15521552            elij(il, i, j) = amax1(0.0, elij(il,i,j))
    … …  
    21442144            ! :            +tra(il,i,j)*(mp(il,i)-mp(il,i+1))
    21452145            ! trap(il,i,j)=trap(il,i,j)/mp(il,i)
    2146             ! end do
     2146            ! END DO
    21472147
    21482148          ELSE
    … …  
    21612161              ! do j=1,ntra
    21622162              ! trap(il,i,j)=trap(il,i+1,j)
    2163               ! end do
     2163              ! END DO
    21642164
    21652165            END IF

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90

    r5082 r5086  
    24312431!!!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
    24322432!!!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
    2433 !!!!      end do
     2433!!!!      END DO
    24342434            elij(il, i, j) = altem
    24352435            elij(il, i, j) = max(0.0, elij(il,i,j))
    … …  
    34243424!AC!        endif ! (i.lt.inb(il) .and. lwork(il))
    34253425!AC!       enddo
    3426 !AC!      end do
     3426!AC!      END DO
    34273427
    34283428400 END DO

  • LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h

    r5075 r5086  
    20102010          read (ilesfile,*) height(k),ugprof(k),vgprof(k),wfls(k),         &
    20112011     &                      dqtdxls(k),dqtdyls(k),dqtdtls(k),thlpcar(k)
    2012         end do
     2012        END DO
    20132013        do k=1,kmax
    20142014          if (height(k) .ne. height1(k)) then
    … …  
    20172017            stop
    20182018          endif
    2019         end do
     2019        END DO
    20202020        close(ilesfile)
    20212021
    … …  
    20352035        do k=1,kmax
    20362036          read (ilesfile,*) height(k),(tracer(k,itrac),itrac=nt1,nt2)
    2037         end do
     2037        END DO
    20382038        close(ilesfile)
    20392039        endif

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90

    r5075 r5086  
    470470                    depth = depth - dzsnSV(ikl,isl) / 2.
    471471                   
    472                 end do
     472                END DO
    473473
    474474            END DO

  • LMDZ6/branches/Amaury_dev/libf/phylmd/isccp_cloud_types.F90

    r1992 r5086  
    10271027    ! write (6,'(a)') '100.*f:'
    10281028    ! write (6,'(8f7.2)') (100.*fluxtop(j,ibox),ibox=1,ncolprint)
    1029     ! end do
     1029    ! END DO
    10301030    ! endif
    10311031
    … …  
    15691569  ! &           (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev)
    15701570  ! &           ,(cchar(acc(ilev,ibox)+1),ilev=1,nlev)
    1571   ! end do
     1571  ! END DO
    15721572  ! close(9)
    15731573

  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake.F90

    r4908 r5086  
    15071507    ! c      do i=1,klon
    15081508    ! c         print*,alpha(i)
    1509     ! c      end do
     1509    ! c      END DO
    15101510    ! cc
    15111511    DO k = 1, klev

  • LMDZ6/branches/Amaury_dev/libf/phylmd/o3_chem_m.F90

    r4103 r5086  
    2424
    2525    integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
    26     real, intent(in):: gmtime ! heure de la journée en fraction de jour
     26    real, intent(in):: gmtime ! heure de la journ�e en fraction de jour
    2727    real, intent(in):: t_seri(:, :) ! (klon, nbp_lev) temperature, in K
    2828
    … …  
    5959    real earth_long
    6060    ! (longitude vraie de la Terre dans son orbite solaire, par
    61     ! rapport au point vernal (21 mars), en degrés)
     61    ! rapport au point vernal (21 mars), en degr�s)
    6262
    6363    real pmu0(klon) ! mean of cosine of solar zenith angle during "pdtphys"
    … …  
    163163    do k =  nbp_lev - 1, 1, -1
    164164       sigma_mass(:, k) = sigma_mass(:, k+1) + zmasse(:, k) * q(:, k)
    165     end do
     165    END DO
    166166
    167167    o3_prod = c + b * q + a6_mass * sigma_mass

  • LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90

    r5075 r5086  
    321321        IF (.NOT. found) THEN
    322322           PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
    323            PRINT*, "          Il prend donc la valeur de surface"
     323           PRINT*, "          Il prEND DOnc la valeur de surface"
    324324           tsoil(:, isoil, :)=ftsol(:, :)
    325325        ENDIF

  • LMDZ6/branches/Amaury_dev/libf/phylmd/physiq_mod.F90

    r5082 r5086  
    17051705998          CONTINUE
    17061706             CLOSE(98)
    1707              CONTINUE
    17081707             IF(nCFMIP>npCFMIP) THEN
    17091708                print*,'nCFMIP > npCFMIP : augmenter npCFMIP et recompiler'

  • LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90

    r5075 r5086  
    182182    do i_v = 1, n_o3_param
    183183       call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))
    184     end do
     184    END DO
    185185
    186186    ! Create the output file and get the variable IDs:
    … …  
    225225            o3_par_out(nbp_lat:1:-1, :, :))
    226226       ! (The order of "rlatu" is inverted in the output file)
    227     end do
     227    END DO
    228228
    229229    call nf95_close(ncid_out)
    … …  
    309309            & varid_out(i))
    310310       call handle_err_copy_att("standard_name")
    311     end do
     311    END DO
    312312
    313313    ! Global attributes:

  • LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_int_m.F90

    r4489 r5086  
    9898                         v3(i, nbp_lev:1:-1))
    9999       ! (invert order of indices because "pplay" is in descending order)
    100     end do
     100    END DO
    101101
    102102  end subroutine regr_pr_int

  • LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_o3_m.F90

    r5075 r5086  
    7777            p3d(1, j, nbp_lev + 1:1:-1), o3_mob_regr(1, j, nbp_lev:1:-1))
    7878       ! (invert order of indices because "p3d" is in descending order)
    79     end do
     79    END DO
    8080
    8181    ! Other latitudes:
    … …  
    8585               p3d(i, j, nbp_lev + 1:1:-1), o3_mob_regr(i, j, nbp_lev:1:-1))
    8686          ! (invert order of indices because "p3d" is in descending order)
    87        end do
    88     end do
     87       END DO
     88    END DO
    8989
    9090    ! Duplicate pole values on all longitudes:

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F

    r5082 r5086  
    2727      REAL zmin,zmax,zllu,zllm
    2828      write(F1000,'(''(4x,'',I3,''(1H-))'')')im+3
    29       DO 10001 i=1,200
     29      DO i=1,200
    3030      jline(1+(i-1)*5:5*i)='.    '
    31 10001 CONTINUE
     31      END DO
    323210002 zmin=z(1,1)
    3333      imin=1
    … …  
    3737      jmax=1
    3838      kzero=0
    39       DO 10003 j=1,jm
    40       DO 10005 i=1,im
     39      DO j=1,jm
     40      DO i=1,im
    4141      IF(.NOT.( z(i,j)>zmax))GOTO 10007
    4242      zmax=z(i,j)
    … …  
    5050      kzero=kzero+1
    515110011 CONTINUE
    52 10005 CONTINUE
     52      END DO
    535310006 CONTINUE
    54 10003 CONTINUE
     54      END DO
    555510004 zsign=(sign(1.,zmin)*sign(1.,zmax)>0.)
    5656      WRITE(*,*)'>>> dump2ds: ',trim(nom_z)
    … …  
    7373      zinf=.false.
    7474      znan=.false.
    75       DO 10017 j=1,jm
    76       DO 10019 i=1,im
     75      DO j=1,jm
     76      DO i=1,im
    7777      az=abs(z(i,j))
    7878      IF(.NOT.( az==0.))GOTO 10021
    … …  
    9797      kchar(i)=32-kchar(i)
    989810027 CONTINUE
    99 10019 CONTINUE
     99      END DO
    10010010020 WRITE(*,'(1x,i3,''|'',1000a)')j+1,(iform(kchar(i)),i=1,im),'|','|
    101101     *'
    102 10017 CONTINUE
     102      END DO
    10310310018 write(*,F1000)
    104104      WRITE(*,'(5x,1000i1)')(mod(i,10),i=1,im)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90

    r5082 r5086  
    353353  do k = 2, n
    354354    w = w * (x - k)
    355   end do
     355  END DO
    356356else
    357357  w = 1
    358358  do k = 0, -n - 1
    359359    y = y * (x + k)
    360   end do
     360  END DO
    361361end if
    362362gamma_res = w / y

  • LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/set99.F

    r5082 r5086  
    1515      NIL=0
    1616      NHL=(N/2)-1
    17       DO 10 K=NIL,NHL
     17      DO K=NIL,NHL
    1818      ANGLE=FLOAT(K)*DEL
    1919      TRIGS(2*K+1)=COS(ANGLE)
    2020      TRIGS(2*K+2)=SIN(ANGLE)
    21    10 CONTINUE
     21      END DO
    2222C
    2323C     FIND FACTORS OF N (8,6,5,4,3,2; ONLY ONE 8 ALLOWED)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/simu_airs.F90

    r5082 r5086  
    3030        do i=1,N-1
    3131        slope(i)=-(T(i+1)-T(i))/(alt(i+1)-alt(i))
    32         end do
     32        END DO
    3333        slope(N)=slope(N-1)
    3434
    … …  
    6060        i=i+i_dir
    6161        if (i<=1.or.i>=N) exit_flag=1
    62         end do
     62        END DO
    6363
    6464        if (first_point<=0) P_tropo=65.4321
Note: See TracChangeset for help on using the changeset viewer.