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

Put .h into modules

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp
Files:
15 edited

Legend:

Unmodified
Added
Removed

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

    r5158 r5160  
    144144    iloc = infind(xsort,xxarr(i),dist=d)
    145145    if (d > tol) then
    146       print *, 'interpolation error'
     146      PRINT *, 'interpolation error'
    147147      stop
    148148    endif

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

    r5099 r5160  
    128128         
    129129  case default
    130     print *, 'Must enter a profile type'
     130    PRINT *, 'Must enter a profile type'
    131131    stop
    132132   

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

    r5099 r5160  
    9797        vu = (1/(0.2714 + 0.00057145*Np*rho_a*1E-6))**2 ! units of Nt = Np*rhoa = #/cm^3
    9898    else
    99         print *, 'Error: Must specify a value for Np in each volume', &
     99        PRINT *, 'Error: Must specify a value for Np in each volume', &
    100100             ' with Morrison/Martin Scheme.'
    101101            stop   
    … …  
    111111      ! vu isn't specified
    112112     
    113       print *, 'Error: Must specify a value for vu for Modified Gamma distribution'
     113      PRINT *, 'Error: Must specify a value for vu for Modified Gamma distribution'
    114114      stop   
    115115     
    … …  
    130130                local_Np = p1 ! total number concentration / pa --- units kg^-1
    131131            else
    132             print *, 'Error: Must specify Np or default value ', &
     132            PRINT *, 'Error: Must specify Np or default value ', &
    133133                 '(p1=Dm [um] or p2=Np [1/kg]) for ', &
    134134                 'Modified Gamma distribution'
    … …  
    169169    else
    170170   
    171     print *, 'Error: Must specify Np or default value ', &
     171    PRINT *, 'Error: Must specify Np or default value ', &
    172172         '(p1=No or p2=lambda) for Exponential distribution'
    173173        stop
    … …  
    191191    if(Np>0) then
    192192   
    193         print *, 'Variable Np not supported for ', &
     193        PRINT *, 'Variable Np not supported for ', &
    194194         'Power Law distribution'
    195195        stop
    … …  
    205205   
    206206        if(Np>0) then
    207         print *, 'Variable Np not supported for ', &
     207        PRINT *, 'Variable Np not supported for ', &
    208208         'Monodispersed distribution'
    209209        stop
    … …  
    225225            log_sigma_g = local_p3
    226226        else
    227             print *, 'Error: Must specify a value for sigma_g ', &
     227            PRINT *, 'Error: Must specify a value for sigma_g ', &
    228228             'when using a Log-Normal distribution'
    229229            stop
    … …  
    244244                local_Np=p1
    245245            else
    246                 print *, 'Error: Must specify Np or default value ', &
     246                PRINT *, 'Error: Must specify Np or default value ', &
    247247                 '(p2=Rg or p1=Np) for Log-Normal distribution'
    248248            endif

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

    r5158 r5160  
    649649   
    650650  if (i /= N_OUT_LIST) then
    651      print *, 'COSP_IO: wrong number of output diagnostics'
    652      print *, i,N_OUT_LIST
     651     PRINT *, 'COSP_IO: wrong number of output diagnostics'
     652     PRINT *, i,N_OUT_LIST
    653653     stop
    654654  endif

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

    r5158 r5160  
    127127        vu = (1/(0.2714 + 0.00057145*Np*rho_a*1E-6))**2.0 ! units of Nt = Np*rhoa = #/cm^3
    128128    else
    129         print *, 'Error: Must specify a value for Np in each volume', &
     129        PRINT *, 'Error: Must specify a value for Np in each volume', &
    130130             ' with Morrison/Martin Scheme.'
    131131            stop   
    … …  
    141141      ! vu isn't specified
    142142     
    143       print *, 'Error: Must specify a value for vu for Modified Gamma distribution'
     143      PRINT *, 'Error: Must specify a value for vu for Modified Gamma distribution'
    144144      stop   
    145145     
    … …  
    185185      else
    186186     
    187         print *, 'Error:  No default value for Dm or Np provided!  '
     187        PRINT *, 'Error:  No default value for Dm or Np provided!  '
    188188        stop
    189189       
    … …  
    255255 
    256256    if(Re>0) then
    257         print *, 'Variable Re not supported for ', &
     257        PRINT *, 'Variable Re not supported for ', &
    258258         'Power-Law distribution'
    259259    stop
    260260    elseif(Np>0) then
    261         print *, 'Variable Np not supported for ', &
     261        PRINT *, 'Variable Np not supported for ', &
    262262         'Power-Law distribution'
    263263    stop
    … …  
    319319      enddo
    320320
    321     ! print *,'test=',ahp,bhp,ahp/(rho_a*Q),D(100),N(100),bpm,dmin_mm,dmax_mm
     321    ! PRINT *,'test=',ahp,bhp,ahp/(rho_a*Q),D(100),N(100),bpm,dmin_mm,dmax_mm
    322322
    323323! ---------------------------------------------------------!
    … …  
    388388    else
    389389
    390       print *, 'Error: Must specify a value for sigma_g'
     390      PRINT *, 'Error: Must specify a value for sigma_g'
    391391      stop
    392392   

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

    r5158 r5160  
    248248
    249249      if ( npart .ne. 4 ) then
    250         print *,'Error in lidar_simulator, npart should be 4, not',npart
     250        PRINT *,'Error in lidar_simulator, npart should be 4, not',npart
    251251        stop
    252252      endif

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

    r5158 r5160  
    152152    path_integral = sumo
    153153  endif
    154   ! print *, sumo
     154  ! PRINT *, sumo
    155155
    156156  return

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

    r5158 r5160  
    141141  if ((.not. gbx%use_reff) .and. (reff_zero)) then ! No Reff in radar. Default in lidar
    142142        gbx%Reff = DEFAULT_LIDAR_REFF
    143         print *, '---------- COSP WARNING ------------'
    144         print *, ''
    145         print *, 'Using default Reff in lidar simulations'
    146         print *, ''
    147         print *, '----------------------------------'
     143        PRINT *, '---------- COSP WARNING ------------'
     144        PRINT *, ''
     145        PRINT *, 'Using default Reff in lidar simulations'
     146        PRINT *, ''
     147        PRINT *, '----------------------------------'
    148148  endif
    149149 
    … …  
    153153  if (Ncolumns == 1) then
    154154     if (gbx%use_precipitation_fluxes) then
    155         print *, '---------- COSP ERROR ------------'
    156         print *, ''
    157         print *, 'Use of precipitation fluxes not supported in CRM mode (Ncolumns=1)'
    158         print *, ''
    159         print *, '----------------------------------'
     155        PRINT *, '---------- COSP ERROR ------------'
     156        PRINT *, ''
     157        PRINT *, 'Use of precipitation fluxes not supported in CRM mode (Ncolumns=1)'
     158        PRINT *, ''
     159        PRINT *, '----------------------------------'
    160160        stop
    161161     endif
    162162     if ((maxval(gbx%dtau_c) > 0.0).or.(maxval(gbx%dem_c) > 0.0)) then
    163         print *, '---------- COSP ERROR ------------'
    164         print *, ''
    165         print *, ' dtau_c > 0.0 or dem_c > 0.0. In CRM mode (Ncolumns=1), '
    166         print *, ' the optical depth (emmisivity) of all clouds must be '
    167         print *, ' passed through dtau_s (dem_s)'
    168         print *, ''
    169         print *, '----------------------------------'
     163        PRINT *, '---------- COSP ERROR ------------'
     164        PRINT *, ''
     165        PRINT *, ' dtau_c > 0.0 or dem_c > 0.0. In CRM mode (Ncolumns=1), '
     166        PRINT *, ' the optical depth (emmisivity) of all clouds must be '
     167        PRINT *, ' passed through dtau_s (dem_s)'
     168        PRINT *, ''
     169        PRINT *, '----------------------------------'
    170170        stop
    171171     endif

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

    r5099 r5160  
    9090!   enddo
    9191!   enddo
    92 !  if (inconsistent)  print *, '%%%% COSP_SIMULATOR: inconsistency in mr_hydro and Reff'
     92!  if (inconsistent)  PRINT *, '%%%% COSP_SIMULATOR: inconsistency in mr_hydro and Reff'
    9393
    9494

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

    r5158 r5160  
    14001400    type(cosp_gridbox),intent(in) :: x
    14011401
    1402     print *, '%%%%----- Information on COSP_GRIDBOX ------'
     1402    PRINT *, '%%%%----- Information on COSP_GRIDBOX ------'
    14031403    ! Scalars and dimensions
    1404     print *,  x%Npoints
    1405     print *,  x%Nlevels
    1406     print *,  x%Ncolumns
    1407     print *,  x%Nhydro
    1408     print *,  x%Nprmts_max_hydro
    1409     print *,  x%Naero
    1410     print *,  x%Nprmts_max_aero
    1411     print *,  x%Npoints_it
     1404    PRINT *,  x%Npoints
     1405    PRINT *,  x%Nlevels
     1406    PRINT *,  x%Ncolumns
     1407    PRINT *,  x%Nhydro
     1408    PRINT *,  x%Nprmts_max_hydro
     1409    PRINT *,  x%Naero
     1410    PRINT *,  x%Nprmts_max_aero
     1411    PRINT *,  x%Npoints_it
    14121412   
    14131413    ! Time [days]
    1414     print *,  x%time
     1414    PRINT *,  x%time
    14151415   
    14161416    ! Radar ancillary info
    1417     print *,  x%radar_freq, &
     1417    PRINT *,  x%radar_freq, &
    14181418            x%k2
    1419     print *,  x%surface_radar, &
     1419    PRINT *,  x%surface_radar, &
    14201420              x%use_mie_tables, &
    14211421              x%use_gas_abs, &
    … …  
    14231423              x%melt_lay
    14241424
    1425 !               print *,  'shape(x%): ',shape(x%)
     1425!               PRINT *,  'shape(x%): ',shape(x%)
    14261426 
    14271427!     type(class_param) ::  hp  ! structure used by radar simulator to store Ze and N scaling constants and other information
    14281428!     type(mie)::  mt           ! structure used by radar simulator to store mie LUT information
    1429     print *,  x%nsizes
     1429    PRINT *,  x%nsizes
    14301430   
    14311431    ! Lidar
    1432     print *,  x%lidar_ice_type
     1432    PRINT *,  x%lidar_ice_type
    14331433   
    14341434    ! Radar
    1435     print *,  x%use_precipitation_fluxes
    1436     print *,  x%use_reff
     1435    PRINT *,  x%use_precipitation_fluxes
     1436    PRINT *,  x%use_reff
    14371437   
    14381438    ! Geolocation (Npoints)
    1439     print *,  'shape(x%longitude): ',shape(x%longitude)
    1440     print *,  'shape(x%latitude): ',shape(x%latitude)
     1439    PRINT *,  'shape(x%longitude): ',shape(x%longitude)
     1440    PRINT *,  'shape(x%latitude): ',shape(x%latitude)
    14411441    ! Gridbox information (Npoints,Nlevels)
    1442     print *,  'shape(x%zlev): ',shape(x%zlev)
    1443     print *,  'shape(x%zlev_half): ',shape(x%zlev_half)
    1444     print *,  'shape(x%dlev): ',shape(x%dlev)
    1445     print *,  'shape(x%p): ',shape(x%p)
    1446     print *,  'shape(x%ph): ',shape(x%ph)
    1447     print *,  'shape(x%T): ',shape(x%T)
    1448     print *,  'shape(x%q): ',shape(x%q)
    1449     print *,  'shape(x%sh): ',shape(x%sh)
    1450     print *,  'shape(x%dtau_s): ',shape(x%dtau_s)
    1451     print *,  'shape(x%dtau_c): ',shape(x%dtau_c)
    1452     print *,  'shape(x%dem_s): ',shape(x%dem_s)
    1453     print *,  'shape(x%dem_c): ',shape(x%dem_c)
    1454     print *,  'shape(x%mr_ozone): ',shape(x%mr_ozone)
     1442    PRINT *,  'shape(x%zlev): ',shape(x%zlev)
     1443    PRINT *,  'shape(x%zlev_half): ',shape(x%zlev_half)
     1444    PRINT *,  'shape(x%dlev): ',shape(x%dlev)
     1445    PRINT *,  'shape(x%p): ',shape(x%p)
     1446    PRINT *,  'shape(x%ph): ',shape(x%ph)
     1447    PRINT *,  'shape(x%T): ',shape(x%T)
     1448    PRINT *,  'shape(x%q): ',shape(x%q)
     1449    PRINT *,  'shape(x%sh): ',shape(x%sh)
     1450    PRINT *,  'shape(x%dtau_s): ',shape(x%dtau_s)
     1451    PRINT *,  'shape(x%dtau_c): ',shape(x%dtau_c)
     1452    PRINT *,  'shape(x%dem_s): ',shape(x%dem_s)
     1453    PRINT *,  'shape(x%dem_c): ',shape(x%dem_c)
     1454    PRINT *,  'shape(x%mr_ozone): ',shape(x%mr_ozone)
    14551455
    14561456    ! Point information (Npoints)
    1457     print *,  'shape(x%land): ',shape(x%land)
    1458     print *,  'shape(x%psfc): ',shape(x%psfc)
    1459     print *,  'shape(x%sunlit): ',shape(x%sunlit)
    1460     print *,  'shape(x%skt): ',shape(x%skt)
    1461     print *,  'shape(x%u_wind): ',shape(x%u_wind)
    1462     print *,  'shape(x%v_wind): ',shape(x%v_wind)
     1457    PRINT *,  'shape(x%land): ',shape(x%land)
     1458    PRINT *,  'shape(x%psfc): ',shape(x%psfc)
     1459    PRINT *,  'shape(x%sunlit): ',shape(x%sunlit)
     1460    PRINT *,  'shape(x%skt): ',shape(x%skt)
     1461    PRINT *,  'shape(x%u_wind): ',shape(x%u_wind)
     1462    PRINT *,  'shape(x%v_wind): ',shape(x%v_wind)
    14631463
    14641464    ! TOTAL and CONV cloud fraction for SCOPS
    1465     print *,  'shape(x%tca): ',shape(x%tca)
    1466     print *,  'shape(x%cca): ',shape(x%cca)
     1465    PRINT *,  'shape(x%tca): ',shape(x%tca)
     1466    PRINT *,  'shape(x%cca): ',shape(x%cca)
    14671467    ! Precipitation fluxes on model levels
    1468     print *,  'shape(x%rain_ls): ',shape(x%rain_ls)
    1469     print *,  'shape(x%rain_cv): ',shape(x%rain_cv)
    1470     print *,  'shape(x%snow_ls): ',shape(x%snow_ls)
    1471     print *,  'shape(x%snow_cv): ',shape(x%snow_cv)
    1472     print *,  'shape(x%grpl_ls): ',shape(x%grpl_ls)
     1468    PRINT *,  'shape(x%rain_ls): ',shape(x%rain_ls)
     1469    PRINT *,  'shape(x%rain_cv): ',shape(x%rain_cv)
     1470    PRINT *,  'shape(x%snow_ls): ',shape(x%snow_ls)
     1471    PRINT *,  'shape(x%snow_cv): ',shape(x%snow_cv)
     1472    PRINT *,  'shape(x%grpl_ls): ',shape(x%grpl_ls)
    14731473    ! Hydrometeors concentration and distribution parameters
    1474     print *,  'shape(x%mr_hydro): ',shape(x%mr_hydro)
    1475     print *,  'shape(x%dist_prmts_hydro): ',shape(x%dist_prmts_hydro)
     1474    PRINT *,  'shape(x%mr_hydro): ',shape(x%mr_hydro)
     1475    PRINT *,  'shape(x%dist_prmts_hydro): ',shape(x%dist_prmts_hydro)
    14761476    ! Effective radius [m]. (Npoints,Nlevels,Nhydro)
    1477     print *,  'shape(x%Reff): ',shape(x%Reff)
    1478     print *,  'shape(x%Np): ',shape(x%Np)       ! added by roj with Quickbeam V3
     1477    PRINT *,  'shape(x%Reff): ',shape(x%Reff)
     1478    PRINT *,  'shape(x%Np): ',shape(x%Np)       ! added by roj with Quickbeam V3
    14791479    ! Aerosols concentration and distribution parameters
    1480     print *,  'shape(x%conc_aero): ',shape(x%conc_aero)
    1481     print *,  'shape(x%dist_type_aero): ',shape(x%dist_type_aero)
    1482     print *,  'shape(x%dist_prmts_aero): ',shape(x%dist_prmts_aero)
     1480    PRINT *,  'shape(x%conc_aero): ',shape(x%conc_aero)
     1481    PRINT *,  'shape(x%dist_type_aero): ',shape(x%dist_type_aero)
     1482    PRINT *,  'shape(x%dist_prmts_aero): ',shape(x%dist_prmts_aero)
    14831483    ! ISCCP simulator inputs
    1484     print *, x%isccp_top_height
    1485     print *, x%isccp_top_height_direction
    1486     print *, x%isccp_overlap
    1487     print *, x%isccp_emsfc_lw
     1484    PRINT *, x%isccp_top_height
     1485    PRINT *, x%isccp_top_height_direction
     1486    PRINT *, x%isccp_overlap
     1487    PRINT *, x%isccp_emsfc_lw
    14881488 
    14891489    ! RTTOV inputs/options
    1490     print *, x%plat
    1491     print *, x%sat
    1492     print *, x%inst
    1493     print *, x%Nchan
    1494     print *,  'shape(x%Ichan): ',x%Ichan
    1495     print *,  'shape(x%Surfem): ',x%Surfem
    1496     print *, x%ZenAng
    1497     print *, x%co2,x%ch4,x%n2o,x%co
     1490    PRINT *, x%plat
     1491    PRINT *, x%sat
     1492    PRINT *, x%inst
     1493    PRINT *, x%Nchan
     1494    PRINT *,  'shape(x%Ichan): ',x%Ichan
     1495    PRINT *,  'shape(x%Surfem): ',x%Surfem
     1496    PRINT *, x%ZenAng
     1497    PRINT *, x%co2,x%ch4,x%n2o,x%co
    14981498               
    14991499END SUBROUTINE COSP_GRIDBOX_PRINT
    … …  
    15021502    type(cosp_misr),intent(in) :: x
    15031503
    1504     print *, '%%%%----- Information on COSP_MISR ------'
     1504    PRINT *, '%%%%----- Information on COSP_MISR ------'
    15051505               
    15061506     ! Dimensions
    1507     print *, x%Npoints
    1508     print *, x%Ntau
    1509     print *, x%Nlevels
     1507    PRINT *, x%Npoints
     1508    PRINT *, x%Ntau
     1509    PRINT *, x%Nlevels
    15101510
    15111511     ! --- (npoints,ntau,nlevels)
    15121512     !  the fraction of the model grid box covered by each of the MISR cloud types
    1513      print *,  'shape(x%fq_MISR): ',shape(x%fq_MISR)
     1513     PRINT *,  'shape(x%fq_MISR): ',shape(x%fq_MISR)
    15141514     
    15151515     ! --- (npoints)
    1516      print *,  'shape(x%MISR_meanztop): ',shape(x%MISR_meanztop)
    1517      print *,  'shape(x%MISR_cldarea): ',shape(x%MISR_cldarea)
     1516     PRINT *,  'shape(x%MISR_meanztop): ',shape(x%MISR_meanztop)
     1517     PRINT *,  'shape(x%MISR_cldarea): ',shape(x%MISR_cldarea)
    15181518     ! --- (npoints,nlevels)
    1519      print *,  'shape(x%MISR_dist_model_layertops): ',shape(x%MISR_dist_model_layertops)
     1519     PRINT *,  'shape(x%MISR_dist_model_layertops): ',shape(x%MISR_dist_model_layertops)
    15201520   
    15211521END SUBROUTINE COSP_MISR_PRINT
    … …  
    15241524    type(cosp_isccp),intent(in) :: x
    15251525           
    1526     print *, x%Npoints
    1527     print *, x%Ncolumns
    1528     print *, x%Nlevels
    1529 
    1530     print *, '%%%%----- Information on COSP_ISCCP ------'
    1531    
    1532      print *, 'shape(x%fq_isccp): ',shape(x%fq_isccp)
    1533      print *, 'shape(x%totalcldarea): ',shape(x%totalcldarea)
    1534      print *, 'shape(x%meantb): ',shape(x%meantb)
    1535      print *, 'shape(x%meantbclr): ',shape(x%meantbclr)
     1526    PRINT *, x%Npoints
     1527    PRINT *, x%Ncolumns
     1528    PRINT *, x%Nlevels
     1529
     1530    PRINT *, '%%%%----- Information on COSP_ISCCP ------'
     1531   
     1532     PRINT *, 'shape(x%fq_isccp): ',shape(x%fq_isccp)
     1533     PRINT *, 'shape(x%totalcldarea): ',shape(x%totalcldarea)
     1534     PRINT *, 'shape(x%meantb): ',shape(x%meantb)
     1535     PRINT *, 'shape(x%meantbclr): ',shape(x%meantbclr)
    15361536     
    1537      print *, 'shape(x%meanptop): ',shape(x%meanptop)
    1538      print *, 'shape(x%meantaucld): ',shape(x%meantaucld)
    1539      print *, 'shape(x%meanalbedocld): ',shape(x%meanalbedocld)
    1540      print *, 'shape(x%boxtau): ',shape(x%boxtau)
    1541      print *, 'shape(x%boxptop): ',shape(x%boxptop)
     1537     PRINT *, 'shape(x%meanptop): ',shape(x%meanptop)
     1538     PRINT *, 'shape(x%meantaucld): ',shape(x%meantaucld)
     1539     PRINT *, 'shape(x%meanalbedocld): ',shape(x%meanalbedocld)
     1540     PRINT *, 'shape(x%boxtau): ',shape(x%boxtau)
     1541     PRINT *, 'shape(x%boxptop): ',shape(x%boxptop)
    15421542END SUBROUTINE COSP_ISCCP_PRINT
    15431543
    … …  
    15451545    type(cosp_vgrid),intent(in) :: x
    15461546           
    1547     print *, '%%%%----- Information on COSP_VGRID ------'
    1548     print *, x%use_vgrid
    1549     print *, x%csat_vgrid
    1550     print *, x%Npoints
    1551     print *, x%Ncolumns
    1552     print *, x%Nlevels
    1553     print *, x%Nlvgrid
     1547    PRINT *, '%%%%----- Information on COSP_VGRID ------'
     1548    PRINT *, x%use_vgrid
     1549    PRINT *, x%csat_vgrid
     1550    PRINT *, x%Npoints
     1551    PRINT *, x%Ncolumns
     1552    PRINT *, x%Nlevels
     1553    PRINT *, x%Nlvgrid
    15541554    ! Array with dimensions (Nlvgrid)
    1555     print *, 'shape(x%z): ',shape(x%z)
    1556     print *, 'shape(x%zl): ',shape(x%zl)
    1557     print *, 'shape(x%zu): ',shape(x%zu)
     1555    PRINT *, 'shape(x%z): ',shape(x%z)
     1556    PRINT *, 'shape(x%zl): ',shape(x%zl)
     1557    PRINT *, 'shape(x%zu): ',shape(x%zu)
    15581558    ! Array with dimensions (Nlevels)
    1559     print *, 'shape(x%mz): ',shape(x%mz)
    1560     print *, 'shape(x%mzl): ',shape(x%mzl)
    1561     print *, 'shape(x%mzu): ',shape(x%mzu)
     1559    PRINT *, 'shape(x%mz): ',shape(x%mz)
     1560    PRINT *, 'shape(x%mzl): ',shape(x%mzl)
     1561    PRINT *, 'shape(x%mzu): ',shape(x%mzu)
    15621562END SUBROUTINE COSP_VGRID_PRINT
    15631563
    … …  
    15651565    type(cosp_sglidar),intent(in) :: x
    15661566           
    1567     print *, '%%%%----- Information on COSP_SGLIDAR ------'
     1567    PRINT *, '%%%%----- Information on COSP_SGLIDAR ------'
    15681568    ! Dimensions
    1569     print *, x%Npoints
    1570     print *, x%Ncolumns
    1571     print *, x%Nlevels
    1572     print *, x%Nhydro
    1573     print *, x%Nrefl
     1569    PRINT *, x%Npoints
     1570    PRINT *, x%Ncolumns
     1571    PRINT *, x%Nlevels
     1572    PRINT *, x%Nhydro
     1573    PRINT *, x%Nrefl
    15741574    ! Arrays with dimensions (Npoints,Nlevels)
    1575     print *, 'shape(x%beta_mol): ',shape(x%beta_mol)
     1575    PRINT *, 'shape(x%beta_mol): ',shape(x%beta_mol)
    15761576    ! Arrays with dimensions (Npoints,Ncolumns,Nlevels)
    1577     print *, 'shape(x%beta_tot): ',shape(x%beta_tot)
    1578     print *, 'shape(x%tau_tot): ',shape(x%tau_tot)
     1577    PRINT *, 'shape(x%beta_tot): ',shape(x%beta_tot)
     1578    PRINT *, 'shape(x%tau_tot): ',shape(x%tau_tot)
    15791579    ! Arrays with dimensions (Npoints,Ncolumns,Nrefl)
    1580     print *, 'shape(x%refl): ',shape(x%refl)
     1580    PRINT *, 'shape(x%refl): ',shape(x%refl)
    15811581END SUBROUTINE COSP_SGLIDAR_PRINT
    15821582
    … …  
    15841584    type(cosp_sgradar),intent(in) :: x
    15851585           
    1586     print *, '%%%%----- Information on COSP_SGRADAR ------'
    1587     print *, x%Npoints
    1588     print *, x%Ncolumns
    1589     print *, x%Nlevels
    1590     print *, x%Nhydro
     1586    PRINT *, '%%%%----- Information on COSP_SGRADAR ------'
     1587    PRINT *, x%Npoints
     1588    PRINT *, x%Ncolumns
     1589    PRINT *, x%Nlevels
     1590    PRINT *, x%Nhydro
    15911591    ! output vertical levels: spaceborne radar -> from TOA to SURFACE
    15921592    ! Arrays with dimensions (Npoints,Nlevels)
    1593     print *, 'shape(x%att_gas): ', shape(x%att_gas)
     1593    PRINT *, 'shape(x%att_gas): ', shape(x%att_gas)
    15941594    ! Arrays with dimensions (Npoints,Ncolumns,Nlevels)
    1595     print *, 'shape(x%Ze_tot): ', shape(x%Ze_tot)
     1595    PRINT *, 'shape(x%Ze_tot): ', shape(x%Ze_tot)
    15961596END SUBROUTINE COSP_SGRADAR_PRINT
    15971597
    … …  
    15991599    type(cosp_radarstats),intent(in) :: x
    16001600           
    1601     print *, '%%%%----- Information on COSP_SGRADAR ------'
    1602     print *, x%Npoints
    1603     print *, x%Ncolumns
    1604     print *, x%Nlevels
    1605     print *, x%Nhydro
    1606     print *, 'shape(x%cfad_ze): ',shape(x%cfad_ze)
    1607     print *, 'shape(x%radar_lidar_tcc): ',shape(x%radar_lidar_tcc)
    1608     print *, 'shape(x%lidar_only_freq_cloud): ',shape(x%lidar_only_freq_cloud)
     1601    PRINT *, '%%%%----- Information on COSP_SGRADAR ------'
     1602    PRINT *, x%Npoints
     1603    PRINT *, x%Ncolumns
     1604    PRINT *, x%Nlevels
     1605    PRINT *, x%Nhydro
     1606    PRINT *, 'shape(x%cfad_ze): ',shape(x%cfad_ze)
     1607    PRINT *, 'shape(x%radar_lidar_tcc): ',shape(x%radar_lidar_tcc)
     1608    PRINT *, 'shape(x%lidar_only_freq_cloud): ',shape(x%lidar_only_freq_cloud)
    16091609END SUBROUTINE COSP_RADARSTATS_PRINT
    16101610
    … …  
    16121612    type(cosp_lidarstats),intent(in) :: x
    16131613           
    1614     print *, '%%%%----- Information on COSP_SGLIDAR ------'
    1615     print *, x%Npoints
    1616     print *, x%Ncolumns
    1617     print *, x%Nlevels
    1618     print *, x%Nhydro
    1619     print *, x%Nrefl
     1614    PRINT *, '%%%%----- Information on COSP_SGLIDAR ------'
     1615    PRINT *, x%Npoints
     1616    PRINT *, x%Ncolumns
     1617    PRINT *, x%Nlevels
     1618    PRINT *, x%Nhydro
     1619    PRINT *, x%Nrefl
    16201620   
    16211621    ! Arrays with dimensions (SR_BINS)
    1622     print *, 'shape(x%srbval): ',shape(x%srbval)
     1622    PRINT *, 'shape(x%srbval): ',shape(x%srbval)
    16231623    ! Arrays with dimensions (Npoints,SR_BINS,Nlevels)
    1624     print *, 'shape(x%cfad_sr): ',shape(x%cfad_sr)
     1624    PRINT *, 'shape(x%cfad_sr): ',shape(x%cfad_sr)
    16251625!    ! Arrays with dimensions (Npoints,Ncolumns,Nlevels) !TIBO
    1626 !    print *, 'shape(x%profSR): ',shape(x%profSR)        !TIBO
     1626!    PRINT *, 'shape(x%profSR): ',shape(x%profSR)        !TIBO
    16271627    ! Arrays with dimensions (Npoints,Nlevels,Ncolumns) !TIBO2
    1628     print *, 'shape(x%profSR): ',shape(x%profSR)        !TIBO2
     1628    PRINT *, 'shape(x%profSR): ',shape(x%profSR)        !TIBO2
    16291629    ! Arrays with dimensions (Npoints,Nlevels)
    1630     print *, 'shape(x%lidarcld): ',shape(x%lidarcld)
    1631     print *, 'shape(x%proftemp): ',shape(x%proftemp)    !TIBO
     1630    PRINT *, 'shape(x%lidarcld): ',shape(x%lidarcld)
     1631    PRINT *, 'shape(x%proftemp): ',shape(x%proftemp)    !TIBO
    16321632    ! Arrays with dimensions (Npoints,LIDAR_NCAT)
    1633     print *, 'shape(x%cldlayer): ',shape(x%cldlayer)
     1633    PRINT *, 'shape(x%cldlayer): ',shape(x%cldlayer)
    16341634    ! Arrays with dimensions (Npoints,LIDAR_NTYPE)            !OPAQ
    1635     print *, 'shape(x%cldtype): ',shape(x%cldtype)            !OPAQ
     1635    PRINT *, 'shape(x%cldtype): ',shape(x%cldtype)            !OPAQ
    16361636    ! Arrays with dimensions (Npoints,PARASOL_NREFL)
    1637     print *, 'shape(x%parasolrefl): ',shape(x%parasolrefl)
     1637    PRINT *, 'shape(x%parasolrefl): ',shape(x%parasolrefl)
    16381638     ! Arrays with dimensions (Npoints,Nlevels,Nphase)
    1639     print *, 'shape(x%lidarcldphase): ',shape(x%lidarcldphase)
     1639    PRINT *, 'shape(x%lidarcldphase): ',shape(x%lidarcldphase)
    16401640     ! Arrays with dimensions (Npoints,Nlevels,LIDAR_NTYPE+1) !OPAQ
    1641     print *, 'shape(x%lidarcldtype): ',shape(x%lidarcldtype)  !OPAQ
     1641    PRINT *, 'shape(x%lidarcldtype): ',shape(x%lidarcldtype)  !OPAQ
    16421642     ! Arrays with dimensions (Npoints,LIDAR_NCAT,Nphase)
    1643     print *, 'shape(x%cldlayerphase): ',shape(x%cldlayerphase)
     1643    PRINT *, 'shape(x%cldlayerphase): ',shape(x%cldlayerphase)
    16441644     ! Arrays with dimensions (Npoints,Ntemps,Nphase)
    1645     print *, 'shape(x%lidarcldphase): ',shape(x%lidarcldtmp)
     1645    PRINT *, 'shape(x%lidarcldphase): ',shape(x%lidarcldtmp)
    16461646
    16471647END SUBROUTINE COSP_LIDARSTATS_PRINT
    … …  
    16501650    type(cosp_subgrid),intent(in) :: x
    16511651           
    1652     print *, '%%%%----- Information on COSP_SUBGRID ------'
    1653     print *, x%Npoints
    1654     print *, x%Ncolumns
    1655     print *, x%Nlevels
    1656     print *, x%Nhydro
    1657    
    1658     print *, 'shape(x%prec_frac): ',shape(x%prec_frac)
    1659     print *, 'shape(x%frac_out): ',shape(x%frac_out)
     1652    PRINT *, '%%%%----- Information on COSP_SUBGRID ------'
     1653    PRINT *, x%Npoints
     1654    PRINT *, x%Ncolumns
     1655    PRINT *, x%Nlevels
     1656    PRINT *, x%Nhydro
     1657   
     1658    PRINT *, 'shape(x%prec_frac): ',shape(x%prec_frac)
     1659    PRINT *, 'shape(x%frac_out): ',shape(x%frac_out)
    16601660END SUBROUTINE COSP_SUBGRID_PRINT
    16611661
    … …  
    16631663    type(cosp_sghydro),intent(in) :: x
    16641664           
    1665     print *, '%%%%----- Information on COSP_SGHYDRO ------'
    1666     print *, x%Npoints
    1667     print *, x%Ncolumns
    1668     print *, x%Nlevels
    1669     print *, x%Nhydro
    1670    
    1671     print *, 'shape(x%mr_hydro): ',shape(x%mr_hydro)
    1672     print *, 'shape(x%Reff): ',shape(x%Reff)
    1673     print *, 'shape(x%Np): ',shape(x%Np)         ! added by roj with Quickbeam V3
     1665    PRINT *, '%%%%----- Information on COSP_SGHYDRO ------'
     1666    PRINT *, x%Npoints
     1667    PRINT *, x%Ncolumns
     1668    PRINT *, x%Nlevels
     1669    PRINT *, x%Nhydro
     1670   
     1671    PRINT *, 'shape(x%mr_hydro): ',shape(x%mr_hydro)
     1672    PRINT *, 'shape(x%Reff): ',shape(x%Reff)
     1673    PRINT *, 'shape(x%Np): ',shape(x%Np)         ! added by roj with Quickbeam V3
    16741674END SUBROUTINE COSP_SGHYDRO_PRINT
    16751675

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

    r5158 r5160  
    273273    endif   
    274274   
    275     if (l_min) print *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
    276     if (l_max) print *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
     275    if (l_min) PRINT *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
     276    if (l_max) PRINT *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
    277277  END SUBROUTINE COSP_CHECK_INPUT_1D
    278278!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    … …  
    308308    endif   
    309309   
    310     if (l_min) print *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
    311     if (l_max) print *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
     310    if (l_min) PRINT *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
     311    if (l_max) PRINT *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
    312312  END SUBROUTINE COSP_CHECK_INPUT_2D
    313313!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    … …  
    343343    endif   
    344344   
    345     if (l_min) print *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
    346     if (l_max) print *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
     345    if (l_min) PRINT *,'----- WARNING: '//trim(pro_name)//': minimum value of '//trim(vname)//' set to: ',min_val
     346    if (l_max) PRINT *,'----- WARNING: '//trim(pro_name)//': maximum value of '//trim(vname)//' set to: ',max_val
    347347  END SUBROUTINE COSP_CHECK_INPUT_3D
    348348

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

    r5158 r5160  
    372372
    373373      if ( Ncat .ne. 4 ) then
    374          print *,'Error in lmd_ipsl_stats.cosp_cldfrac, Ncat must be 4, not',Ncat
     374         PRINT *,'Error in lmd_ipsl_stats.cosp_cldfrac, Ncat must be 4, not',Ncat
    375375         stop
    376376      endif
    … …  
    918918             enddo
    919919             checkcldlayerphase2=cldlayer(i,iz)-checkcldlayerphase
    920              if( (checkcldlayerphase2.gt.0.01).or.(checkcldlayerphase2.lt.-0.01) ) print *, checkcldlayerphase,cldlayer(i,iz)
     920             if( (checkcldlayerphase2.gt.0.01).or.(checkcldlayerphase2.lt.-0.01) ) PRINT *, checkcldlayerphase,cldlayer(i,iz)
    921921
    922922          endif
    … …  
    962962
    963963    if(checktemp.NE.lidarcldtemp(i,itemp,1))then
    964       print *, i,itemp
    965       print *, lidarcldtemp(i,itemp,1:4)
     964      PRINT *, i,itemp
     965      PRINT *, lidarcldtemp(i,itemp,1:4)
    966966    endif
    967967

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

    r5158 r5160  
    512512  alam=3E5/freq
    513513  if((alam < wlmin) .or. (alam > wlmax)) then
    514     print *, 'm_ice: wavelength out of bounds'
     514    PRINT *, 'm_ice: wavelength out of bounds'
    515515    stop
    516516  endif

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

    r5158 r5160  
    223223!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    224224! AI mars 2017
    225 !        print *, 'Allocating memory for gridbox type...'
     225!        PRINT *, 'Allocating memory for gridbox type...'
    226226
    227227! Surafce emissivity
    … …  
    240240!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    241241
    242 !        print *, 'Populating input structure...'
     242!        PRINT *, 'Populating input structure...'
    243243        gbx%longitude = lon
    244244        gbx%latitude = lat
    … …  
    324324        ! Define new vertical grid
    325325!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    326 !        print *, 'Defining new vertical grid...'
     326!        PRINT *, 'Defining new vertical grid...'
    327327        call construct_cosp_vgrid(gbx,Nlr,use_vgrid,csat_vgrid,vgrid)
    328328
    … …  
    330330       ! Allocate memory for other types
    331331!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    332 !        print *, 'Allocating memory for other types...'
     332!        PRINT *, 'Allocating memory for other types...'
    333333        call construct_cosp_subgrid(Npoints, Ncolumns, Nlevels, sgx)
    334334        call construct_cosp_sgradar(cfg,Npoints,Ncolumns,Nlevels,N_HYDRO,sgradar)
    … …  
    347347
    348348      !$OMP MASTER
    349 !        print *, ' Open outpts files and define axis'
     349!        PRINT *, ' Open outpts files and define axis'
    350350        call cosp_output_open(Nlevlmdz, Ncolumns, presnivs, dtime, freq_cosp, &
    351351                              ok_mensuelCOSP, ok_journeCOSP, ok_hfCOSP, ok_all_xml, &
    … …  
    358358        ! Call simulator
    359359!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    360 !        print *, 'Calling simulator...'
     360!        PRINT *, 'Calling simulator...'
    361361!! AI
    362362!        call cosp(overlap,Ncolumns,cfg,vgrid,gbx,sgx,sgradar,sglidar,isccp,misr,stradar,stlidar)
    … …  
    371371!!!!!!!!!!!!!!!!!! Ecreture des sorties Cosp !!!!!!!!!!!!!!r!!!!!!:!!!!!
    372372
    373 !       print *, 'Calling write output'
     373!       PRINT *, 'Calling write output'
    374374     if (.NOT. debut_cosp) call cosp_output_write(Nlevlmdz, Npoints, Ncolumns, itap, dtime, freq_COSP, missing_val, &
    375375                               cfg, gbx, vgrid, sglidar, sgradar, stlidar, stradar, &
    … …  
    379379        ! Deallocate memory in derived types
    380380!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    381 !        print *, 'Deallocating memory...'
     381!        PRINT *, 'Deallocating memory...'
    382382        call free_cosp_gridbox(gbx)
    383383        call free_cosp_subgrid(sgx)
    … …  
    397397  ! Time in s. Only for testing purposes
    398398!  call system_clock(t1,count_rate,count_max)
    399 !  print *,(t1-t0)*1.0/count_rate
     399!  PRINT *,(t1-t0)*1.0/count_rate
    400400    if (debut_cosp) then
    401401      debut_cosp=.false.

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

    r5158 r5160  
    213213          else
    214214            if (Np_matrix(tp,pr,k)>0) then
    215                print *, 'Warning: Re and Np set for the same ', &
     215               PRINT *, 'Warning: Re and Np set for the same ', &
    216216                        'volume & hydrometeor type.  Np is being ignored.'
    217217            endif
Note: See TracChangeset for help on using the changeset viewer.