Changeset 3125
- Timestamp:
- Dec 14, 2017, 8:55:21 AM (7 years ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90
r3002 r3125 45 45 USE ioipsl_getin_p_mod, ONLY: getin_p 46 46 USE slab_heat_transp_mod, ONLY: ini_slab_transp_geom 47 #ifdef REPROBUS 48 USE CHEM_REP, ONLY : Init_chem_rep_phys 49 #endif 47 50 IMPLICIT NONE 48 51 … … 177 180 rlonudyn,rlatudyn,rlonvdyn,rlatvdyn) 178 181 #endif 182 IF (type_trac == 'repr') THEN 183 #ifdef REPROBUS 184 CALL Init_chem_rep_phys(klon_omp,nbp_lev) 185 #endif 186 END IF 179 187 END IF 180 188 -
LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
r3124 r3125 1730 1730 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac(:) 1731 1731 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_trac_cum(:) 1732 #ifdef REPROBUS 1733 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_nas(:) 1734 #endif 1732 1735 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_vdf(:) 1733 1736 TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_the(:) -
LMDZ6/trunk/libf/phylmd/phys_output_mod.F90
r3082 r3125 35 35 USE iophy 36 36 USE dimphy 37 USE infotrac_phy, ONLY: nqtot, nqo, niadv, tname, ttext 37 USE infotrac_phy, ONLY: nqtot, nqo, niadv, tname, ttext, type_trac 38 38 USE ioipsl 39 39 USE phys_cal_mod, only : hour, calend … … 46 46 USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt 47 47 USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref 48 #ifdef REPROBUS 49 USE chem_rep, ONLY: nbnas, tnamenas, ttextnas 50 #endif 48 51 #ifdef CPP_XIOS 49 52 ! ug Pour les sorties XIOS … … 156 159 IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot)) 157 160 IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot)) 161 #ifdef REPROBUS 162 IF (.NOT. ALLOCATED(o_nas)) ALLOCATE(o_nas(nbnas)) 163 #endif 158 164 ALLOCATE(o_dtr_the(nqtot),o_dtr_con(nqtot),o_dtr_lessi_impa(nqtot)) 159 165 ALLOCATE(o_dtr_lessi_nucl(nqtot),o_dtr_insc(nqtot),o_dtr_bcscav(nqtot)) … … 572 578 ENDDO 573 579 ENDIF 580 IF (type_trac=='repr') THEN 581 #ifdef REPROBUS 582 DO iiq=1,nbnas 583 o_nas(iiq) = ctrl_out((/ 4, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), & 584 tnamenas(iiq),ttextnas(iiq), "-", & 585 (/ '', '', '', '', '', '', '', '', '', '' /)) 586 ENDDO 587 #endif 588 ENDIF 574 589 575 590 ENDDO ! iff -
LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90
r3124 r3125 168 168 o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, & 169 169 o_wxT, o_uxu, o_vxv, o_TxT, o_trac, & 170 #ifdef REPROBUS 171 o_nas, & 172 #endif 170 173 o_dtr_vdf, o_dtr_the, o_dtr_con, & 171 174 o_dtr_lessi_impa, o_dtr_lessi_nucl, & … … 316 319 surf_PM25_sulf, tau_strat_550, tausum_strat, & 317 320 vsed_aer, tau_strat_1020, f_r_wet 321 #endif 322 323 #ifdef REPROBUS 324 USE CHEM_REP, ONLY : nas, nbnas, tnamenas, ttextnas 318 325 #endif 319 326 … … 2051 2058 ENDDO 2052 2059 ENDIF 2060 #ifndef REPROBUS 2053 2061 CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d) 2062 #endif 2054 2063 ENDIF 2055 2064 ENDDO 2056 2065 ENDIF 2066 2067 IF (type_trac == 'repr') THEN 2068 #ifdef REPROBUS 2069 DO iq=1,nbnas 2070 CALL histwrite_phy(o_nas(iq), nas(:,:,iq)) 2071 ENDDO 2072 #endif 2073 ENDIF 2074 2057 2075 ENDIF !(iflag_phytrac==1) 2058 2076 -
LMDZ6/trunk/libf/phylmd/phytrac_mod.F90
r3109 r3125 508 508 flag_cvltr(it)=.false. 509 509 ENDIF 510 511 CASE('repr') 512 flag_cvltr(it)=.false. 510 513 511 514 CASE('inca') -
LMDZ6/trunk/libf/phylmd/tracreprobus_mod.F90
r2320 r3125 72 72 WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr 73 73 IF (it < nqtot) THEN 74 WRITE(*,*)'ici ',it,nqtot74 WRITE(*,*)'iciav',it,nqtot 75 75 #ifdef REPROBUS 76 CALL minmaxqfi_chimie( tr_seri(1,1,it),0.,1.e10,'avant chimie ')76 CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'avant chimie ') 77 77 #endif 78 WRITE(*,*) 'apres'78 WRITE(*,*)iter,'avpres' 79 79 ENDIF 80 80 ENDDO … … 94 94 ! pdel, pdt_rep, daynum : definit dans phytrac et utilise dans chemmain 95 95 ! et transporte par CHEM_REP 96 97 DO it=1, nbtr 98 ! WRITE(lunout,*)it,' ',minval(tr_seri(:,:,it)),maxval(tr_seri(:,:,it)) 99 ! seulement pour les especes chimiques (pas l'age de l'air) 100 ! verif valeurs extremes 101 ! correction: a 1.e-30 quand =0 ou negatif et 102 ! call abort si >ou= 1.e10 103 WRITE(*,*)it,'nqtot',nqtot,'nbtr',nbtr 104 IF (it < nqtot) THEN 105 WRITE(*,*)'iciap',it,nqtot 106 CALL minmaxqfi_chimie(it,tr_seri(1,1,it),0.,1.e10,'apres chemmain') 107 WRITE(*,*)iter,'appres' 108 ENDIF 109 ENDDO 110 96 111 #endif 97 112
Note: See TracChangeset
for help on using the changeset viewer.