[1335] | 1 | ! |
---|
| 2 | ! phys_local_var_mod.F90 1327 2010-03-17 15:33:56Z idelkadi $ |
---|
| 3 | |
---|
[1938] | 4 | MODULE phys_output_var_mod |
---|
[1335] | 5 | |
---|
[1938] | 6 | use dimphy |
---|
| 7 | ! Variables outputs pour les ecritures des sorties |
---|
| 8 | !====================================================================== |
---|
| 9 | ! |
---|
| 10 | ! |
---|
| 11 | !====================================================================== |
---|
| 12 | ! Declaration des variables |
---|
[1335] | 13 | |
---|
[1938] | 14 | REAL, SAVE, ALLOCATABLE :: snow_o(:), zfra_o(:) |
---|
| 15 | !$OMP THREADPRIVATE(snow_o, zfra_o) |
---|
| 16 | INTEGER, SAVE, ALLOCATABLE :: itau_con(:) ! Nombre de pas ou rflag <= 1 |
---|
| 17 | !$OMP THREADPRIVATE(itau_con) |
---|
| 18 | REAL, ALLOCATABLE :: bils_ec(:) ! Contribution of energy conservation |
---|
[2042] | 19 | REAL, ALLOCATABLE :: bils_ech(:) ! Contribution of energy conservation |
---|
[1938] | 20 | REAL, ALLOCATABLE :: bils_tke(:) ! Contribution of energy conservation |
---|
| 21 | REAL, ALLOCATABLE :: bils_diss(:) ! Contribution of energy conservation |
---|
| 22 | REAL, ALLOCATABLE :: bils_kinetic(:) ! bilan de chaleur au sol, kinetic |
---|
| 23 | REAL, ALLOCATABLE :: bils_enthalp(:) ! bilan de chaleur au sol |
---|
| 24 | REAL, ALLOCATABLE :: bils_latent(:) ! bilan de chaleur au sol |
---|
[2042] | 25 | !$OMP THREADPRIVATE(bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent) |
---|
[1753] | 26 | |
---|
[2580] | 27 | ! Marine |
---|
| 28 | ! Variables de sortie du simulateur AIRS |
---|
[1865] | 29 | |
---|
[2580] | 30 | REAL, ALLOCATABLE :: map_prop_hc(:),map_prop_hist(:) |
---|
| 31 | REAL, ALLOCATABLE :: alt_tropo(:) |
---|
| 32 | REAL, ALLOCATABLE :: & |
---|
| 33 | & map_emis_hc(:),map_iwp_hc(:),map_deltaz_hc(:),& |
---|
| 34 | & map_pcld_hc(:),map_tcld_hc(:) |
---|
| 35 | REAL, ALLOCATABLE :: map_emis_hist(:),map_iwp_hist(:),map_deltaz_hist(:),& |
---|
| 36 | & map_rad_hist(:) |
---|
| 37 | REAL, ALLOCATABLE :: map_ntot(:),map_hc(:),map_hist(:) |
---|
| 38 | REAL, ALLOCATABLE :: map_Cb(:),map_ThCi(:),map_Anv(:) |
---|
| 39 | REAL, ALLOCATABLE :: map_emis_Cb(:),map_pcld_Cb(:),map_tcld_Cb(:) |
---|
| 40 | REAL, ALLOCATABLE :: map_emis_ThCi(:),map_pcld_ThCi(:),map_tcld_ThCi(:) |
---|
| 41 | REAL, ALLOCATABLE :: map_emis_Anv(:),map_pcld_Anv(:),map_tcld_Anv(:) |
---|
| 42 | |
---|
| 43 | |
---|
[1938] | 44 | ! ug Plein de variables venues de phys_output_mod |
---|
| 45 | INTEGER, PARAMETER :: nfiles = 9 |
---|
| 46 | LOGICAL, DIMENSION(nfiles), SAVE :: clef_files |
---|
| 47 | LOGICAL, DIMENSION(nfiles), SAVE :: clef_stations |
---|
| 48 | INTEGER, DIMENSION(nfiles), SAVE :: lev_files |
---|
| 49 | INTEGER, DIMENSION(nfiles), SAVE :: nid_files |
---|
| 50 | INTEGER, DIMENSION(nfiles), SAVE :: nnid_files |
---|
| 51 | !$OMP THREADPRIVATE(clef_files, clef_stations, lev_files,nid_files,nnid_files) |
---|
| 52 | INTEGER, DIMENSION(nfiles), SAVE :: nnhorim |
---|
[1761] | 53 | |
---|
[1938] | 54 | INTEGER, DIMENSION(nfiles), SAVE :: nhorim, nvertm |
---|
| 55 | INTEGER, DIMENSION(nfiles), SAVE :: nvertap, nvertbp, nvertAlt |
---|
| 56 | REAL, DIMENSION(nfiles), SAVE :: zoutm |
---|
| 57 | CHARACTER(LEN=20), DIMENSION(nfiles), SAVE :: type_ecri |
---|
| 58 | !$OMP THREADPRIVATE(nnhorim, nhorim, nvertm, zoutm,type_ecri) |
---|
| 59 | CHARACTER(LEN=20), DIMENSION(nfiles), SAVE :: type_ecri_files, phys_out_filetypes |
---|
| 60 | !$OMP THREADPRIVATE(type_ecri_files, phys_out_filetypes) |
---|
| 61 | CHARACTER(LEN=20), DIMENSION(nfiles), SAVE :: phys_out_filenames |
---|
| 62 | !$OMP THREADPRIVATE(phys_out_filenames) |
---|
[1791] | 63 | |
---|
[1938] | 64 | ! swaero_diag : flag indicates if it is necessary to do calculation for some aerosol diagnostics |
---|
[2529] | 65 | !--OB: this needs to be set to TRUE by default and changed back to FALSE after first radiation call |
---|
| 66 | !-- and corrected back to TRUE based on output requests |
---|
| 67 | LOGICAL, SAVE :: swaero_diag=.TRUE. |
---|
[1938] | 68 | !$OMP THREADPRIVATE(swaero_diag) |
---|
[1791] | 69 | |
---|
[1938] | 70 | INTEGER, SAVE:: levmin(nfiles) = 1 |
---|
| 71 | INTEGER, SAVE:: levmax(nfiles) |
---|
| 72 | !$OMP THREADPRIVATE(levmin, levmax) |
---|
[1791] | 73 | |
---|
[1938] | 74 | REAL, SAVE :: zdtime_moy |
---|
| 75 | !$OMP THREADPRIVATE(zdtime_moy) |
---|
[1797] | 76 | |
---|
[1938] | 77 | LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL |
---|
[1807] | 78 | |
---|
[2333] | 79 | REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon) |
---|
| 80 | REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon) |
---|
[1938] | 81 | REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon) |
---|
[2341] | 82 | !$OMP THREADPRIVATE(zustr_gwd_hines, zvstr_gwd_hines) |
---|
| 83 | !$OMP THREADPRIVATE(zustr_gwd_front, zvstr_gwd_front) |
---|
[1938] | 84 | !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando) |
---|
| 85 | |
---|
[1791] | 86 | TYPE ctrl_out |
---|
| 87 | INTEGER,DIMENSION(nfiles) :: flag |
---|
| 88 | CHARACTER(len=20) :: name |
---|
| 89 | CHARACTER(len=150) :: description |
---|
| 90 | CHARACTER(len=20) :: unit |
---|
| 91 | CHARACTER(len=20),DIMENSION(nfiles) :: type_ecrit |
---|
| 92 | END TYPE ctrl_out |
---|
[1938] | 93 | |
---|
[2538] | 94 | REAL, SAVE, ALLOCATABLE :: sens_prec_liq_o(:,:), sens_prec_sol_o(:,:) |
---|
| 95 | REAL, SAVE, ALLOCATABLE :: lat_prec_liq_o(:,:), lat_prec_sol_o(:,:) |
---|
| 96 | !$OMP THREADPRIVATE(sens_prec_liq_o, sens_prec_sol_o,lat_prec_liq_o,lat_prec_sol_o) |
---|
| 97 | |
---|
[1335] | 98 | CONTAINS |
---|
| 99 | |
---|
[1938] | 100 | !====================================================================== |
---|
| 101 | SUBROUTINE phys_output_var_init |
---|
| 102 | use dimphy |
---|
[1335] | 103 | |
---|
[1938] | 104 | IMPLICIT NONE |
---|
[1335] | 105 | |
---|
[1938] | 106 | include "clesphys.h" |
---|
[1335] | 107 | |
---|
[1938] | 108 | !------------------------------------------------ |
---|
[1335] | 109 | |
---|
[1938] | 110 | allocate(snow_o(klon), zfra_o(klon)) |
---|
| 111 | allocate(itau_con(klon)) |
---|
[2538] | 112 | allocate(sens_prec_liq_o(klon,2)) |
---|
| 113 | allocate(sens_prec_sol_o(klon,2)) |
---|
| 114 | allocate(lat_prec_liq_o(klon,2)) |
---|
| 115 | allocate(lat_prec_sol_o(klon,2)) |
---|
| 116 | sens_prec_liq_o = 0.0 ; sens_prec_sol_o = 0.0 |
---|
| 117 | lat_prec_liq_o = 0.0 ; lat_prec_sol_o = 0.0 |
---|
| 118 | |
---|
[2042] | 119 | allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon)) |
---|
[1335] | 120 | |
---|
[2580] | 121 | ! Marine |
---|
| 122 | ! Variables de sortie simulateur AIRS |
---|
| 123 | |
---|
[2582] | 124 | if (ok_airs) then |
---|
[2580] | 125 | allocate (map_prop_hc(klon),map_prop_hist(klon)) |
---|
| 126 | allocate (alt_tropo(klon)) |
---|
[2582] | 127 | allocate (map_emis_hc(klon),map_iwp_hc(klon),map_deltaz_hc(klon)) |
---|
[2580] | 128 | allocate (map_pcld_hc(klon),map_tcld_hc(klon)) |
---|
| 129 | allocate (map_emis_hist(klon),map_iwp_hist(klon),map_deltaz_hist(klon)) |
---|
| 130 | allocate (map_rad_hist(klon)) |
---|
| 131 | allocate (map_ntot(klon),map_hc(klon),map_hist(klon)) |
---|
| 132 | allocate (map_Cb(klon),map_ThCi(klon),map_Anv(klon)) |
---|
| 133 | allocate (map_emis_Cb(klon),map_pcld_Cb(klon),map_tcld_Cb(klon)) |
---|
| 134 | allocate (map_emis_ThCi(klon),map_pcld_ThCi(klon),map_tcld_ThCi(klon)) |
---|
| 135 | allocate (map_emis_Anv(klon),map_pcld_Anv(klon),map_tcld_Anv(klon)) |
---|
[2582] | 136 | endif |
---|
[2580] | 137 | |
---|
[2333] | 138 | IF (ok_hines) allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon)) |
---|
| 139 | IF (.not.ok_hines.and.ok_gwd_rando) & |
---|
| 140 | allocate(zustr_gwd_front(klon), zvstr_gwd_front(klon)) |
---|
[1938] | 141 | IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon)) |
---|
[1335] | 142 | |
---|
[1938] | 143 | END SUBROUTINE phys_output_var_init |
---|
[1335] | 144 | |
---|
[1938] | 145 | !====================================================================== |
---|
| 146 | SUBROUTINE phys_output_var_end |
---|
| 147 | use dimphy |
---|
| 148 | IMPLICIT NONE |
---|
| 149 | |
---|
[2582] | 150 | include "clesphys.h" |
---|
| 151 | |
---|
[1938] | 152 | deallocate(snow_o,zfra_o,itau_con) |
---|
[2042] | 153 | deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent) |
---|
[1938] | 154 | |
---|
[2580] | 155 | ! Marine |
---|
| 156 | ! Variables de sortie simulateur AIRS |
---|
| 157 | |
---|
[2582] | 158 | if (ok_airs) then |
---|
[2580] | 159 | deallocate (map_prop_hc,map_prop_hist) |
---|
| 160 | deallocate (alt_tropo) |
---|
[2582] | 161 | deallocate (map_emis_hc,map_iwp_hc,map_deltaz_hc) |
---|
[2580] | 162 | deallocate (map_pcld_hc,map_tcld_hc) |
---|
| 163 | deallocate (map_emis_hist,map_iwp_hist,map_deltaz_hist) |
---|
| 164 | deallocate (map_rad_hist) |
---|
| 165 | deallocate (map_ntot,map_hc,map_hist) |
---|
| 166 | deallocate (map_Cb,map_ThCi,map_Anv) |
---|
| 167 | deallocate (map_emis_Cb,map_pcld_Cb,map_tcld_Cb) |
---|
| 168 | deallocate (map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi) |
---|
| 169 | deallocate (map_emis_Anv,map_pcld_Anv,map_tcld_Anv) |
---|
[2582] | 170 | endif |
---|
[2580] | 171 | |
---|
[1938] | 172 | END SUBROUTINE phys_output_var_end |
---|
| 173 | |
---|
[1335] | 174 | END MODULE phys_output_var_mod |
---|