Ignore:
Timestamp:
Jun 15, 2021, 1:18:14 PM (3 years ago)
Author:
crisi
Message:

replace files by symbloic liks from phylmdiso towards phylmd.
Many files at once

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmdiso/phys_output_ctrlout_mod.F90

    r3927 r3940  
    11!
    2 ! $Id: phys_output_ctrlout_mod.F90 3691 2020-05-30 15:37:19Z oboucher $
     2! $Id: phys_output_ctrlout_mod.F90 3888 2021-05-05 10:50:37Z jyg $
    33!
    44MODULE phys_output_ctrlout_mod
    … …  
    77  USE indice_sol_mod
    88  USE aero_mod
    9 
    109
    1110  IMPLICIT NONE
    … …  
    273272    't2m_sic', "Temp 2m "//clnsurf(4), "K", (/ ('', i=1, 10) /)) /)
    274273
     274  TYPE(ctrl_out), SAVE :: o_nt2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     275    'nt2mout', 'Nbt2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     276  TYPE(ctrl_out), SAVE :: o_nq2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     277    'nq2mout', 'Nbq2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     278  TYPE(ctrl_out), SAVE :: o_nu2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     279    'nu2mout', 'Nbu2m out of range complete computation', '-', (/ ('', i=1, 10) /))
     280
     281  TYPE(ctrl_out), SAVE :: o_nt2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     282    'nt2moutfg', 'Nbt2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
     283  TYPE(ctrl_out), SAVE :: o_nq2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     284    'nq2moutfg', 'Nbq2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
     285  TYPE(ctrl_out), SAVE :: o_nu2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), &
     286    'nu2moutfg', 'Nbu2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /))
     287
    275288  TYPE(ctrl_out), SAVE :: o_gusts = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11, 11/), &
    276289    'gusts', 'surface gustiness', 'm2/s2', (/ ('', i=1, 10) /))
    … …  
    281294    'wind100m', '100-m wind speed', 'm/s', (/ ('', i=1, 10) /))
    282295  TYPE(ctrl_out), SAVE :: o_loadfactor_wind_onshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    283     'load_factor_wind_onshore', 'Load factor for onshore windmill', '-', (/ ('', i=1, 10) /))
     296    'woncfr', 'Onshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
    284297  TYPE(ctrl_out), SAVE :: o_loadfactor_wind_offshore = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    285     'load_factor_wind_offshore', 'Load factor for offshore windmill', '-', (/ ('', i=1, 10) /))
     298    'wofcfr', 'Offshore Wind Capacity factor', 'kW/kW_installed', (/ ('', i=1, 10) /))
    286299  TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
    287300    'wind10max', '10m wind speed max', 'm/s', &
    … …  
    462475  TYPE(ctrl_out), SAVE :: o_SWupSFCcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
    463476    'SWupSFCcleanclr', 'SWup clear sky clean (no aerosol) at surface', 'W/m2', (/ ('', i=1, 10) /))
    464   TYPE(ctrl_out), SAVE :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), &
     477  TYPE(ctrl_out), SAVE :: o_fdiffSWdnSFC = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     478    'fdiffSWdnSFC', 'Fraction of diffuse SWdn at surface', 'W/m2', (/ ('', i=1, 10) /))
     479  TYPE(ctrl_out), SAVE :: o_SWdnSFC = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11, 11/), &
    465480    'SWdnSFC', 'SWdn at surface', 'W/m2', (/ ('', i=1, 10) /))
    466481  TYPE(ctrl_out), SAVE :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11, 11/), &
    … …  
    567582  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_evappot_srf  = (/ &
    568583      ctrl_out((/ 1, 6, 10, 10, 10, 10, 11, 11, 11, 11/),'evappot_ter',    &
    569       "Temperature"//clnsurf(1),"K", (/ ('', i=1, 10) /)),   &
     584      "Potential evaporation "//clnsurf(1),"kg/(m2*s)", (/ ('', i=1, 10) /)),   &
    570585      ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11, 11/),'evappot_lic',    &
    571       "Temperature"//clnsurf(2),"K", (/ ('', i=1, 10) /)),   &
     586      "Potential evaporation "//clnsurf(2),"kg/(m2*s)", (/ ('', i=1, 10) /)),   &
    572587      ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11, 11/),'evappot_oce',    &
    573       "Temperature"//clnsurf(3),"K", (/ ('', i=1, 10) /)),   &
     588      "Potential evaporation "//clnsurf(3),"kg/(m2*s)", (/ ('', i=1, 10) /)),   &
    574589      ctrl_out((/ 4, 6, 10, 10, 10, 10, 11, 11, 11, 11/),'evappot_sic',    &
    575       "Temperature"//clnsurf(4),"K", (/ ('', i=1, 10) /)) /)
     590      "Potential evaporation "//clnsurf(4),"kg/(m2*s)", (/ ('', i=1, 10) /)) /)
    576591
    577592  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_sens_srf     = (/          &
    … …  
    695710    'iwp', 'Cloud ice water path', 'kg/m2', (/ ('', i=1, 10) /))
    696711  TYPE(ctrl_out), SAVE :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    697     'ue', 'Zonal dry static energy transport', '-', (/ ('', i=1, 10) /))
     712    'ue', 'Zonal dry static energy transport', 'J/m/s', (/ ('', i=1, 10) /))
    698713  TYPE(ctrl_out), SAVE :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    699     've', 'Merid dry static energy transport', '-', (/ ('', i=1, 10) /))
     714    've', 'Merid dry static energy transport', 'J/m/s', (/ ('', i=1, 10) /))
    700715  TYPE(ctrl_out), SAVE :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    701     'uq', 'Zonal humidity transport', '-', (/ ('', i=1, 10) /))
     716    'uq', 'Zonal humidity transport', 'kg/m/s', (/ ('', i=1, 10) /))
    702717  TYPE(ctrl_out), SAVE :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    703     'vq', 'Merid humidity transport', '-', (/ ('', i=1, 10) /))
     718    'vq', 'Merid humidity transport', 'kg/m/s', (/ ('', i=1, 10) /))
    704719  TYPE(ctrl_out), SAVE :: o_uwat = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    705     'uwat', 'Zonal total water transport', '-', (/ ('', i=1, 10) /))
     720    'uwat', 'Zonal total water transport', 'kg/m/s', (/ ('', i=1, 10) /))
    706721  TYPE(ctrl_out), SAVE :: o_vwat = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    707     'vwat', 'Merid total water transport', '-', (/ ('', i=1, 10) /))
     722    'vwat', 'Merid total water transport', 'kg/m/s', (/ ('', i=1, 10) /))
    708723  TYPE(ctrl_out), SAVE :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    709724    'cape', 'Conv avlbl pot ener', 'J/kg', (/ ('', i=1, 10) /))
    … …  
    799814'flat_w', 'flat within_wake', 'W/m2', (/ ('', i=1, 10) /))
    800815!!
    801   type(ctrl_out),save :: o_delta_tsurf    = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    802 'delta_tsurf', 'Temperature difference (w-x)', 'K', (/ ('', i=1, 10) /))
    803816  type(ctrl_out),save :: o_cdragh_x       = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    804817'cdragh_x', 'cdragh off-wake', '', (/ ('', i=1, 10) /))
    … …  
    10041017  TYPE(ctrl_out), SAVE :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    10051018    'tke ', 'TKE', 'm2/s2', (/ ('', i=1, 10) /))
     1019  TYPE(ctrl_out), SAVE :: o_tke_dissip = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1020    'tke_dissip ', 'TKE DISSIPATION', 'm2/s3', (/ ('', i=1, 10) /))   
    10061021  TYPE(ctrl_out), SAVE :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    10071022    'tke_max', 'TKE max', 'm2/s2',                                  &
    10081023      (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &
    10091024         't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))
    1010 
    10111025  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_srf      = (/             &
    10121026      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'tke_ter',       &
    … …  
    10501064      "min PBL mixing length "//clnsurf(4),"m", (/ ('', i=1, 10) /)) /)
    10511065
     1066
    10521067  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_tke_max_srf  = (/                          &
    10531068      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'tke_max_ter',                &
    … …  
    10771092      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'dltpbltke_sic',       &
    10781093      "TKE difference (w - x) "//clnsurf(4),"-", (/ ('', i=1, 10) /)) /)
     1094
     1095  TYPE(ctrl_out), SAVE :: o_delta_tsurf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1096    'delta_tsurf ', 'T_surf difference (w - x)', 'K', (/ ('', i=1, 10) /))
     1097  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_delta_tsurf_srf      = (/             &
     1098      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'delta_tsurf_ter',       &
     1099      "T_surf difference (w - x) "//clnsurf(1),"-", (/ ('', i=1, 10) /)), &
     1100      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'delta_tsurf_lic',       &
     1101      "T_surf difference (w - x) "//clnsurf(2),"-", (/ ('', i=1, 10) /)), &
     1102      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'delta_tsurf_oce',       &
     1103      "T_surf difference (w - x) "//clnsurf(3),"-", (/ ('', i=1, 10) /)), &
     1104      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11, 11/),'delta_tsurf_sic',       &
     1105      "T_surf difference (w - x) "//clnsurf(4),"-", (/ ('', i=1, 10) /)) /)
    10791106
    10801107  TYPE(ctrl_out), SAVE :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    … …  
    13041331  TYPE(ctrl_out), SAVE :: o_flx_co2_land = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
    13051332    'flx_co2_land', 'CO2 flux from the land', '1', (/ ('', i=1, 10) /))
     1333  TYPE(ctrl_out), SAVE :: o_flx_co2_ocean_cor = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1334    'flx_co2_ocean_cor', 'correction of the CO2 flux from the ocean', 'kg CO2 m-2 s-1', (/ ('', i=1, 10) /))
     1335  TYPE(ctrl_out), SAVE :: o_flx_co2_land_cor = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     1336    'flx_co2_land_cor', 'correction of the CO2 flux from the land', 'kg CO2 m-2 s-1', (/ ('', i=1, 10) /))
    13061337
    13071338#ifdef CPP_StratAer
    … …  
    14331464  TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14341465    'rhum', 'Relative humidity', '-', (/ ('', i=1, 10) /))
     1466  TYPE(ctrl_out), SAVE :: o_rhl = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1467    'rhl', 'Relative humidity wrt liquid', '%', (/ ('', i=1, 10) /))
     1468  TYPE(ctrl_out), SAVE :: o_rhi = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     1469    'rhi', 'Relative humidity wrt ice', '%', (/ ('', i=1, 10) /))
    14351470  TYPE(ctrl_out), SAVE :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
    14361471    'ozone', 'Ozone mole fraction', '-', (/ ('', i=1, 10) /))
    … …  
    19702005#endif
    19712006
     2007   type(ctrl_out), save:: o_delta_sst &
     2008        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_SST', &
     2009        "ocean-air interface temperature minus bulk SST", "K", '')
     2010
     2011   type(ctrl_out), save:: o_delta_sal &
     2012        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'delta_sal', &
     2013        "ocean-air interface salinity minus bulk salinity", "ppt", '')
     2014
     2015   type(ctrl_out), save:: o_ds_ns &
     2016        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'dS_ns', &
     2017        "subskin salinity minus foundation salinity", "ppt", '')
     2018
     2019   type(ctrl_out), save:: o_dt_ns &
     2020        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'dT_ns', &
     2021        "subskin temperature minus foundation temperature", "K", '')
     2022
     2023   type(ctrl_out), save:: o_dter &
     2024        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'dTer', &
     2025        "ocean-air interface temperature minus sub-skin temperature", "K", '')
     2026
     2027   type(ctrl_out), save:: o_dser &
     2028        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'dSer', &
     2029        "ocean-air interface salinity minus sub-skin salinity", "ppt", '')
     2030
     2031   type(ctrl_out), save:: o_tkt &
     2032        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'tkt', &
     2033        "thickness of thermal microlayer", "m", '')
     2034
     2035   type(ctrl_out), save:: o_tks &
     2036        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'tks', &
     2037        "thickness of salinity microlayer", "m", '')
     2038
     2039   type(ctrl_out), save:: o_taur &
     2040        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'taur', &
     2041        "momentum flux due to rain", "Pa", '')
     2042
     2043   type(ctrl_out), save:: o_sss &
     2044        = ctrl_out([1, 10, 10, 1, 10, 10, 11, 11, 11, 11], 'SSS', &
     2045        "bulk sea-surface salinity", "ppt", '')
     2046
    19722047END MODULE phys_output_ctrlout_mod
Note: See TracChangeset for help on using the changeset viewer.