Ignore:
Timestamp:
Aug 2, 2024, 2:12:03 PM (3 months ago)
Author:
abarral
Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver.F90

    r4853 r5158  
    311311  tstart = omp_get_wtime()
    312312#endif
    313   do jrepeat = 1,driver_config%nrepeat
     313  DO jrepeat = 1,driver_config%nrepeat
    314314   
    315315    if (driver_config%do_parallel) then
    … …  
    321321     
    322322      !$OMP PARALLEL DO PRIVATE(istartcol, iendcol) SCHEDULE(RUNTIME)
    323       do jblock = 1, nblock
     323      DO jblock = 1, nblock
    324324        ! Specify the range of columns to process.
    325325        istartcol = (jblock-1) * driver_config%nblocksize &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver_config.F90

    r4773 r5158  
    400400    ! wavelength bound, noting that the number of diagnostics is one
    401401    ! fewer than the number of valid bounds
    402     do jdiag = 0,NMaxSpectralDiag
     402    DO jdiag = 0,NMaxSpectralDiag
    403403      if (this%sw_diag_wavelength_bound(jdiag+1) < 0.0_jprb) then
    404404        this%n_sw_diag = max(0,jdiag-1)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver_read_input.F90

    r4773 r5158  
    563563
    564564    ! Loop through all radiatively important gases
    565     do jgas = 1,NMaxGases
     565    DO jgas = 1,NMaxGases
    566566      if (jgas == IH2O) then
    567567        if (file%exists('q')) then

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_ifs_driver.F90

    r4773 r5158  
    369369  tstart = omp_get_wtime()
    370370#endif
    371   do jrepeat = 1,driver_config%nrepeat
     371  DO jrepeat = 1,driver_config%nrepeat
    372372
    373373!    if (driver_config%do_parallel) then
    … …  
    379379
    380380      !$OMP PARALLEL DO PRIVATE(istartcol, iendcol) SCHEDULE(RUNTIME)
    381       do jblock = 1, nblock
     381      DO jblock = 1, nblock
    382382        ! Specify the range of columns to process.
    383383        istartcol = (jblock-1) * driver_config%nblocksize &

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_ifs_driver_blocked.F90

    r4773 r5158  
    371371  tstart = omp_get_wtime()
    372372#endif
    373   do jrepeat = 1,driver_config%nrepeat
     373  DO jrepeat = 1,driver_config%nrepeat
    374374
    375375!    if (driver_config%do_parallel) then
    … …  
    378378      !$OMP PARALLEL DO SCHEDULE(DYNAMIC,1)&
    379379      !$OMP&PRIVATE(JRL,IBEG,IEND,IL,IB)
    380       do jrl=1,ncol,nproma
     380      DO jrl=1,ncol,nproma
    381381        ibeg=jrl
    382382        iend=min(ibeg+nproma-1,ncol)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ifs_blocking.F90

    r4773 r5158  
    329329  !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
    330330  !$OMP&PRIVATE(IB,IFLD)
    331   do ib=1,ngpblks
    332     do ifld=1,ifs_config%ifldstot
     331  DO ib=1,ngpblks
     332    DO ifld=1,ifs_config%ifldstot
    333333      zrgp(:,ifld,ib) = 0._jprb
    334334    enddo
    … …  
    351351    !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
    352352    !$OMP&PRIVATE(JRL,IBEG,IEND,IL,IB,JAER,JOFF,JLEV,JALB)
    353     do jrl=1,ncol,nproma
     353    DO jrl=1,ncol,nproma
    354354
    355355      ibeg=jrl
    … …  
    363363      zrgp(1:il,ifs_config%iamu0,ib)  =  single_level%cos_sza(ibeg:iend)   ! cosine of solar zenith ang (mu0)
    364364
    365       do jemiss=1,yderad%nlwemiss
     365      DO jemiss=1,yderad%nlwemiss
    366366        zrgp(1:il,ifs_config%iemiss+jemiss-1,ib)  =  single_level%lw_emissivity(ibeg:iend,jemiss)
    367367      enddo
    … …  
    378378      ! zrgp(1:il,islon,ib)    = ???
    379379
    380       do jalb=1,yderad%nsw
     380      DO jalb=1,yderad%nsw
    381381        zrgp(1:il,ifs_config%iald+jalb-1,ib)  =  single_level%sw_albedo(ibeg:iend,jalb)
    382382      enddo
    383383
    384384      if (allocated(single_level%sw_albedo_direct)) then
    385         do jalb=1,yderad%nsw
     385        DO jalb=1,yderad%nsw
    386386          zrgp(1:il,ifs_config%ialp+jalb-1,ib)  =  single_level%sw_albedo_direct(ibeg:iend,jalb)
    387387        end do
    388388      else
    389         do jalb=1,yderad%nsw
     389        DO jalb=1,yderad%nsw
    390390          zrgp(1:il,ifs_config%ialp+jalb-1,ib)  =  single_level%sw_albedo(ibeg:iend,jalb)
    391391        end do
    392392      end if
    393393     
    394       do jlev=1,nlev
     394      DO jlev=1,nlev
    395395        zrgp(1:il,ifs_config%iti+jlev-1,ib)   = temperature_fl(ibeg:iend,jlev) ! full level temperature
    396396        zrgp(1:il,ifs_config%ipr+jlev-1,ib)   = pressure_fl(ibeg:iend,jlev) ! full level pressure
    … …  
    398398      enddo
    399399
    400       do jlev=1,nlev
     400      DO jlev=1,nlev
    401401        zrgp(1:il,ifs_config%iwv+jlev-1,ib)   = gas%mixing_ratio(ibeg:iend,jlev,IH2O) ! this is already in MassMixingRatio units
    402402        if (rad_config%do_clouds) then
    … …  
    421421      if (yderad%naermacc == 1) then
    422422        joff=ifs_config%iaero
    423         do jaer=1,rad_config%n_aerosol_types
    424           do jlev=1,nlev
     423        DO jaer=1,rad_config%n_aerosol_types
     424          DO jlev=1,nlev
    425425            zrgp(1:il,joff,ib) = aerosol%mixing_ratio(ibeg:iend,jlev,jaer)
    426426            joff=joff+1
    … …  
    429429      endif
    430430
    431       do jlev=1,nlev+1
     431      DO jlev=1,nlev+1
    432432        ! zrgp(1:il,ihpr+jlev-1,ib)  = ???
    433433        zrgp(1:il,ifs_config%iaprs+jlev-1,ib) = thermodynamics%pressure_hl(ibeg:iend,jlev)
    … …  
    453453      ! local workaround variables for standalone input files
    454454      if (rad_config%do_clouds) then
    455         do jlev=1,nlev
     455        DO jlev=1,nlev
    456456          ! missing full-level temperature and pressure as well as land-sea-mask
    457457          zrgp(1:il,ifs_config%ire_liq+jlev-1,ib) = cloud%re_liq(ibeg:iend,jlev)
    458458          zrgp(1:il,ifs_config%ire_ice+jlev-1,ib) = cloud%re_ice(ibeg:iend,jlev)
    459459        enddo
    460         do jlev=1,nlev-1
     460        DO jlev=1,nlev-1
    461461          ! for the love of it, I can't figure this one out. Probably to do with
    462462          ! my crude approach of setting PGEMU?
    … …  
    465465        if(present(iseed)) iseed(1:il,ib) = single_level%iseed(ibeg:iend)
    466466      else
    467         do jlev=1,nlev
     467        DO jlev=1,nlev
    468468          ! missing full-level temperature and pressure as well as land-sea-mask
    469469          zrgp(1:il,ifs_config%ire_liq+jlev-1,ib) = 0._jprb
    470470          zrgp(1:il,ifs_config%ire_ice+jlev-1,ib) = 0._jprb
    471471        enddo
    472         do jlev=1,nlev-1
     472        DO jlev=1,nlev-1
    473473          zrgp(1:il,ifs_config%ioverlap+jlev-1,ib) = 0._jprb
    474474        enddo
    … …  
    531531    !$OMP PARALLEL DO SCHEDULE(RUNTIME)&
    532532    !$OMP&PRIVATE(JRL,IBEG,IEND,IL,IB,JLEV,JG)
    533     do jrl=1,ncol,nproma
     533    DO jrl=1,ncol,nproma
    534534      ibeg=jrl
    535535      iend=min(ibeg+nproma-1,ncol)
    … …  
    537537      ib=(jrl-1)/nproma+1
    538538
    539       do jlev=1,nlev+1
     539      DO jlev=1,nlev+1
    540540        flux%sw_up(ibeg:iend,jlev) = zrgp(1:il,ifs_config%ifrso+jlev-1,ib)
    541541        flux%lw_up(ibeg:iend,jlev) = zrgp(1:il,ifs_config%ifrth+jlev-1,ib)
    … …  
    561561      emissivity_out(ibeg:iend) = zrgp(1:il,ifs_config%iemit,ib)
    562562      if (yradiation%yrerad%lapproxswupdate) then
    563         do jg=1,yradiation%yrerad%nsw
     563        DO jg=1,yradiation%yrerad%nsw
    564564          flux_diffuse_band(ibeg:iend,jg) = zrgp(1:il,ifs_config%iswdiffuseband+jg-1,ib)
    565565          flux_direct_band(ibeg:iend,jg) = zrgp(1:il,ifs_config%iswdirectband+jg-1,ib)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_cloud_generator.F90

    r4773 r5158  
    5757       &  od_scaling, total_cloud_cover)
    5858
    59   do jlev = 1,nlev
    60     do jcol = 1,ncol
     59  DO jlev = 1,nlev
     60    DO jcol = 1,ncol
    6161      write(6,'(f5.2,a)',advance='no') od_scaling(jcol,jlev), ' '
    6262    end do

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_random_number_generator.F90

    r4773 r5158  
    2929  call random_number_generator%initialize(itype=IRngType, iseed=212075152, &
    3030       &                                  nmaxstreams=streammax)
    31   do jl = 1,1
     31  DO jl = 1,1
    3232    print *, 'initial_state = [ ', int(random_number_generator%istate(1:streammax),jpib), ' ]'
    3333    call random_number_generator%uniform_distribution(vec)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_spartacus_math.F90

    r4773 r5158  
    4545
    4646  write(*,*) 'A ='
    47   do j = 1,m
     47  DO j = 1,m
    4848     write(*,*) A(1,j,:)
    4949  end do
    5050
    5151  write(*,*) 'B ='
    52   do j = 1,m
     52  DO j = 1,m
    5353     write(*,*) B(1,j,:)
    5454  end do
    … …  
    6060
    6161  write(*,*) 'C = A*B ='
    62   do j = 1,m
     62  DO j = 1,m
    6363     write(*,*) C(1,j,:)
    6464  end do
    … …  
    7070
    7171  write(*,*) 'C = A\B ='
    72   do j = 1,m
     72  DO j = 1,m
    7373     write(*,*) C(1,j,:)
    7474  end do
    … …  
    7878
    7979  write(*,*) 'expm(A) ='
    80   do j = 1,m
     80  DO j = 1,m
    8181     write(*,*) A(1,j,:)
    8282  end do
    … …  
    9393   
    9494    write(*,*) 'fast_expm(A) = '
    95     do j = 1,m
     95    DO j = 1,m
    9696      write(*,*) A(1,j,:)
    9797    end do
    … …  
    106106   
    107107    write(*,*) 'expm(A) = '
    108     do j = 1,m
     108    DO j = 1,m
    109109      write(*,*) A(1,j,:)
    110110    end do
    … …  
    115115   
    116116    write(*,*) 'expm(zeros) = '
    117     do j = 1,m
     117    DO j = 1,m
    118118      write(*,*) A(1,j,:)
    119119    end do
    … …  
    124124   
    125125    write(*,*) 'fast_expm(A) = '
    126     do j = 1,m
     126    DO j = 1,m
    127127      write(*,*) A(1,j,:)
    128128    end do
    … …  
    140140   
    141141    write(*,*) 'expm(A) = '
    142     do j = 1,m
     142    DO j = 1,m
    143143      write(*,*) A(1,j,:)
    144144    end do
    … …  
    149149   
    150150    write(*,*) 'expm(zeros) = '
    151     do j = 1,m
     151    DO j = 1,m
    152152      write(*,*) A(1,j,:)
    153153    end do
Note: See TracChangeset for help on using the changeset viewer.