- Timestamp:
- Nov 5, 2018, 3:24:59 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
/LMDZ6/trunk removed
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_output_ctrlout_mod.F90
r3356 r3411 1 !2 ! $Id$3 !4 1 MODULE phys_output_ctrlout_mod 5 2 … … 19 16 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 20 17 ! CHARACTER(len=20), dimension(nfiles) :: TEF = type_ecri_files 21 22 !!! saving lon and lat as variables for CMIP6 DataRequest23 TYPE(ctrl_out), SAVE :: o_longitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &24 'io_lon', '', '', (/ ('once', i=1, 10) /))25 TYPE(ctrl_out), SAVE :: o_latitude = ctrl_out((/ 1, 1, 1, 1, 1, 1, 11, 11, 11, 11/), &26 'io_lat', '', '', (/ ('once', i=1, 10) /))27 18 28 19 !!! Comosantes de la coordonnee sigma-hybride … … 54 45 55 46 !!! 2D 56 TYPE(ctrl_out), SAVE :: o_sza = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), &57 'sza', 'Solar zenithal angle', 'degrees', (/ ('', i=1, 10) /))58 47 59 48 ! Marine … … 335 324 'ndayrain', 'Number of dayrain(liq+sol)', '-', & 336 325 (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /)) 337 TYPE(ctrl_out), SAVE :: o_rain_fall = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11, 11/), &338 'rain_fall', 'Precip Totale liq', 'kg/(s*m2)', (/ ('', i=1, 10) /))339 TYPE(ctrl_out), SAVE :: o_rain_con = ctrl_out((/ 7, 7, 7, 10, 7, 10, 11, 11, 11, 11/), &340 'rain_con', 'Precip liq conv.', 'kg/(s*m2)', (/ ('', i=1, 10) /))341 326 TYPE(ctrl_out), SAVE :: o_precip = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11, 11/), & 342 327 'precip', 'Precip Totale liq+sol', 'kg/(s*m2)', (/ ('', i=1, 10) /)) … … 396 381 TYPE(ctrl_out), SAVE :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11, 11/), & 397 382 'SWupTOAclr', 'SWup clear sky at TOA', 'W/m2', (/ ('', i=1, 10) /)) 398 TYPE(ctrl_out), SAVE :: o_SWupTOAcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &399 'SWupTOAcleanclr', 'SWup clear sky clean (no aerosol) at TOA', 'W/m2', (/ ('', i=1, 10) /))400 383 TYPE(ctrl_out), SAVE :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11, 11/), & 401 384 'SWdnTOA', 'SWdn at TOA', 'W/m2', (/ ('', i=1, 10) /)) … … 444 427 TYPE(ctrl_out), SAVE :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11, 11/), & 445 428 'SWupSFCclr', 'SWup clear sky at surface', 'W/m2', (/ ('', i=1, 10) /)) 446 TYPE(ctrl_out), SAVE :: o_SWupSFCcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &447 'SWupSFCcleanclr', 'SWup clear sky clean (no aerosol) at surface', 'W/m2', (/ ('', i=1, 10) /))448 429 TYPE(ctrl_out), SAVE :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), & 449 430 'SWdnSFC', 'SWdn at surface', 'W/m2', (/ ('', i=1, 10) /)) 450 431 TYPE(ctrl_out), SAVE :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11, 11/), & 451 432 'SWdnSFCclr', 'SWdn clear sky at surface', 'W/m2', (/ ('', i=1, 10) /)) 452 TYPE(ctrl_out), SAVE :: o_SWdnSFCcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &453 'SWdnSFCcleanclr', 'SWdn clear sky clean (no aerosol) at surface', 'W/m2', (/ ('', i=1, 10) /))454 433 TYPE(ctrl_out), SAVE :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11, 11/), & 455 434 'LWupSFC', 'Upwd. IR rad. at surface', 'W/m2', (/ ('', i=1, 10) /)) … … 460 439 TYPE(ctrl_out), SAVE :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 5, 10, 11, 11, 11, 11/), & 461 440 'LWdnSFCclr', 'Down. CS IR rad. at surface', 'W/m2', (/ ('', i=1, 10) /)) 462 TYPE(ctrl_out), SAVE :: o_LWupTOAcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &463 'LWupTOAcleanclr', 'Upward CS clean (no aerosol) IR rad. at TOA', 'W/m2', (/ ('', i=1, 10) /))464 TYPE(ctrl_out), SAVE :: o_LWdnSFCcleanclr = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &465 'LWdnSFCcleanclr', 'Downward CS clean (no aerosol) IR rad. at surface', 'W/m2', (/ ('', i=1, 10) /))466 441 TYPE(ctrl_out), SAVE :: o_bils = ctrl_out((/ 1, 2, 10, 5, 10, 10, 11, 11, 11, 11/), & 467 442 'bils', 'Surf. total heat flux', 'W/m2', (/ ('', i=1, 10) /)) … … 490 465 TYPE(ctrl_out), SAVE :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 491 466 'fqfonte', 'Land ice melt', 'kg/m2/s', (/ ('', i=1, 10) /)) 492 TYPE(ctrl_out), SAVE :: o_mrroli = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &493 'mrroli', 'Runoff flux over land ice', 'kg/m2/s', (/ ('', i=1, 10) /))494 467 TYPE(ctrl_out), SAVE :: o_runofflic = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 495 468 'runofflic', 'Land ice melt to ocean', 'kg/m2/s', (/ ('', i=1, 10) /)) … … 679 652 'iwp', 'Cloud ice water path', 'kg/m2', (/ ('', i=1, 10) /)) 680 653 TYPE(ctrl_out), SAVE :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 681 'ue', 'Zonal dry staticenergy transport', '-', (/ ('', i=1, 10) /))654 'ue', 'Zonal energy transport', '-', (/ ('', i=1, 10) /)) 682 655 TYPE(ctrl_out), SAVE :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 683 've', 'Merid dry staticenergy transport', '-', (/ ('', i=1, 10) /))656 've', 'Merid energy transport', '-', (/ ('', i=1, 10) /)) 684 657 TYPE(ctrl_out), SAVE :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 685 658 'uq', 'Zonal humidity transport', '-', (/ ('', i=1, 10) /)) 686 659 TYPE(ctrl_out), SAVE :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 687 660 'vq', 'Merid humidity transport', '-', (/ ('', i=1, 10) /)) 688 TYPE(ctrl_out), SAVE :: o_uwat = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &689 'uwat', 'Zonal total water transport', '-', (/ ('', i=1, 10) /))690 TYPE(ctrl_out), SAVE :: o_vwat = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &691 'vwat', 'Merid total water transport', '-', (/ ('', i=1, 10) /))692 661 TYPE(ctrl_out), SAVE :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 693 662 'cape', 'Conv avlbl pot ener', 'J/kg', (/ ('', i=1, 10) /)) … … 818 787 's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 10) /)) 819 788 820 !-------Déclenchement stochastiqu e789 !-------Déclenchement stochastiqu 821 790 TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), & 822 791 'proba_notrig', 'Probabilite de non-declenchement', ' ', (/ ('', i=1, 10) /)) … … 972 941 TYPE(ctrl_out), SAVE :: o_rh2m = ctrl_out((/ 5, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 973 942 'rh2m', 'Relative humidity at 2m', '%', (/ ('', i=1, 10) /)) 974 !TYPE(ctrl_out), SAVE :: o_rh2m_min = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &975 !'rh2m_min', 'Min Relative humidity at 2m', '%', &976 !(/ 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', &977 !'t_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)' /))978 !TYPE(ctrl_out), SAVE :: o_rh2m_max = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &979 !'rh2m_max', 'Max Relative humidity at 2m', '%', &980 !(/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', &981 !'t_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /))943 TYPE(ctrl_out), SAVE :: o_rh2m_min = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 944 'rh2m_min', 'Min Relative humidity at 2m', '%', & 945 (/ 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', & 946 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)', 't_min(X)' /)) 947 TYPE(ctrl_out), SAVE :: o_rh2m_max = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 948 'rh2m_max', 'Max Relative humidity at 2m', '%', & 949 (/ 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', & 950 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)', 't_max(X)' /)) 982 951 TYPE(ctrl_out), SAVE :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 983 952 'qsat2m', 'Saturant humidity at 2m', '%', (/ ('', i=1, 10) /)) … … 1071 1040 TYPE(ctrl_out), SAVE :: o_SWnetOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1072 1041 'SWnetOR', 'Sfce net SW radiation OR', 'W/m2', (/ ('', i=1, 10) /)) 1042 TYPE(ctrl_out), SAVE :: o_SWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1043 'SWdownOR', 'Sfce incident SW radiation OR', 'W/m2', (/ ('', i=1, 10) /)) 1073 1044 TYPE(ctrl_out), SAVE :: o_LWdownOR = ctrl_out((/ 10, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1074 1045 'LWdownOR', 'Sfce incident LW radiation OR', 'W/m2', (/ ('', i=1, 10) /)) … … 1316 1287 'vsed_aer', 'Strat. aerosol sedimentation velocity (mass-weighted)', 'm/s', (/ ('', i=1, 10) /)) 1317 1288 !--budget 2D 1318 TYPE(ctrl_out), SAVE :: o_budg_dep_dry_ocs = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &1319 'budg_dep_dry_ocs', 'OCS dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))1320 TYPE(ctrl_out), SAVE :: o_budg_dep_wet_ocs = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &1321 'budg_dep_wet_ocs', 'OCS wet deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))1322 TYPE(ctrl_out), SAVE :: o_budg_dep_dry_so2 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &1323 'budg_dep_dry_so2', 'SO2 dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))1324 TYPE(ctrl_out), SAVE :: o_budg_dep_wet_so2 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &1325 'budg_dep_wet_so2', 'SO2 wet deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))1326 1289 TYPE(ctrl_out), SAVE :: o_budg_dep_dry_h2so4 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), & 1327 1290 'budg_dep_dry_h2so4', 'H2SO4 dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /)) … … 1454 1417 TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1455 1418 'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /)) 1456 TYPE(ctrl_out), SAVE :: o_icc3dcon = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &1457 'icc3dcon', 'Convective cloud ice fraction', '1', (/ ('', i=1, 10) /))1458 TYPE(ctrl_out), SAVE :: o_icc3dstra = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &1459 'icc3dstra', 'Stratiform cloud ice fraction', '1', (/ ('', i=1, 10) /))1460 TYPE(ctrl_out), SAVE :: o_cldicemxrat = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &1461 'cldicemxrat', 'Cloud Ice Mixing Ratio', '1', (/ ('', i=1, 10) /))1462 TYPE(ctrl_out), SAVE :: o_cldwatmxrat = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &1463 'cldwatmxrat', 'Cloud Water Mixing Ratio', '1', (/ ('', i=1, 10) /))1464 TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &1465 'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /))1466 1419 TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 1467 1420 'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 10) /)) … … 1512 1465 TYPE(ctrl_out), SAVE :: o_mc = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1513 1466 'mc', 'Convective mass flux', 'kg/m2/s', (/ ('', i=1, 10) /)) 1514 TYPE(ctrl_out), SAVE :: o_ftime_deepcv = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &1515 'ftime_deepcv', 'Fraction of time deep convection Occurs', ' ', (/ ('', i=1, 10) /))1516 1467 TYPE(ctrl_out), SAVE :: o_ftime_con = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1517 'ftime_con', 'Fraction of time convection Occurs', ' ', (/ ('', i=1, 10) /)) 1518 !!jyg 'ftime_con', 'Fraction of time convection Occurs', ' ', & 1519 !!jyg (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', & 1520 !!jyg 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) 1468 'ftime_con', 'Fraction of time convection Occurs', ' ', & 1469 (/ 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', & 1470 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)', 'inst(X)' /)) 1521 1471 TYPE(ctrl_out), SAVE :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1522 1472 'dtdyn', 'Dynamics dT', 'K/s', (/ ('', i=1, 10) /)) … … 1553 1503 TYPE(ctrl_out), SAVE :: o_dqwak2d = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1554 1504 'dqwak2d', 'Wake dQ', '(kg/m2)/s', (/ ('', i=1, 10) /)) 1555 TYPE(ctrl_out), SAVE :: o_cv_gen = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &1556 'cv_gen', 'Cumulonimbus genesis', '1/(m2 s)', (/ ('', i=1, 10) /))1557 1505 TYPE(ctrl_out), SAVE :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1558 1506 'wake_h', 'wake_h', '-', (/ ('', i=1, 10) /)) 1559 TYPE(ctrl_out), SAVE :: o_wake_dens = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &1560 'wake_dens', 'number of wakes per m2', '1/m2', (/ ('', i=1, 10) /))1561 1507 TYPE(ctrl_out), SAVE :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1562 1508 'wake_s', 'wake_s', '-', (/ ('', i=1, 10) /)) … … 1566 1512 'wake_deltaq', 'wake_deltaq', ' ', (/ ('', i=1, 10) /)) 1567 1513 TYPE(ctrl_out), SAVE :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1568 'wake_omg', 'wake_omg', ' Pa/s', (/ ('', i=1, 10) /))1514 'wake_omg', 'wake_omg', ' ', (/ ('', i=1, 10) /)) 1569 1515 TYPE(ctrl_out), SAVE :: o_wdtrainA = ctrl_out((/ 4, 5, 10, 4, 10, 10, 11, 11, 11, 11 /), & 1570 1516 'wdtrainA', 'precipitation from AA', '-', (/ ('', i=1, 10) /)) … … 1744 1690 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac(:) 1745 1691 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac_cum(:) 1746 #ifdef REPROBUS1747 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_nas(:)1748 #endif1749 1692 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_vdf(:) 1750 1693 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_the(:) … … 1772 1715 TYPE(ctrl_out), SAVE :: o_rsucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1773 1716 'rsucs', 'SW CS upward radiation', 'W m-2', (/ ('', i=1, 10) /)) 1774 TYPE(ctrl_out), SAVE :: o_rsucsaf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &1775 'rsucsaf', 'SW CS clean (no aerosol) upward radiation', 'W m-2', (/ ('', i=1, 10) /))1776 1717 TYPE(ctrl_out), SAVE :: o_rsdcs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1777 1718 'rsdcs', 'SW CS downward radiation', 'W m-2', (/ ('', i=1, 10) /)) 1778 TYPE(ctrl_out), SAVE :: o_rsdcsaf = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &1779 'rsdcsaf', 'SW CS clean (no aerosol) downward radiation', 'W m-2', (/ ('', i=1, 10) /))1780 1719 TYPE(ctrl_out), SAVE :: o_rlucs = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), & 1781 1720 'rlucs', 'LW CS upward radiation', 'W m-2', (/ ('', i=1, 10) /))
Note: See TracChangeset
for help on using the changeset viewer.