Changeset 5101 for LMDZ6/branches/Amaury_dev/libf/phylmd
- Timestamp:
- Jul 23, 2024, 8:22:55 AM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd
- Files:
-
- 192 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/coarsemission.F
r5099 r5101 28 28 USE indice_sol_mod 29 29 USE infotrac 30 USE dustemission_mod, ONLY : dustemission31 ! USE phytracr_spl_mod, ONLY 30 USE dustemission_mod, ONLY: dustemission 31 ! USE phytracr_spl_mod, ONLY: nbreg_dust, nbreg_ind, nbreg_bb 32 32 IMPLICIT NONE 33 33 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/condsurfc.f90
r5099 r5101 4 4 lmt_omnat) 5 5 USE dimphy 6 USE netcdf, ONLY 6 USE netcdf, ONLY: nf90_close, nf90_noerr, nf90_inq_varid, nf90_open, nf90_nowrite, nf90_get_var 7 7 IMPLICIT none 8 8 … … 38 38 write(6, *)' Pb d''ouverture du fichier limitbc.nc' 39 39 write(6, *)' ierr = ', ierr 40 callexit(1)40 CALL exit(1) 41 41 endif 42 42 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/condsurfc_new.f90
r5099 r5101 6 6 USE mod_phys_lmdz_para 7 7 USE dimphy 8 USE netcdf, ONLY 8 USE netcdf, ONLY: nf90_get_var, nf90_close, nf90_noerr, nf90_inq_varid, nf90_open, nf90_nowrite 9 9 IMPLICIT none 10 10 … … 61 61 write(6, *)' Pb d''ouverture du fichier limitbc.nc' 62 62 write(6, *)' ierr = ', ierr 63 callexit(1)63 CALL exit(1) 64 64 endif 65 65 … … 192 192 !$OMP END MASTER 193 193 !$OMP BARRIER 194 callscatter(lmt_bcff_glo, lmt_bcff)195 callscatter(lmt_bcnff_glo, lmt_bcnff)196 callscatter(lmt_bcbbl_glo, lmt_bcbbl)197 callscatter(lmt_bcbbh_glo, lmt_bcbbh)198 callscatter(lmt_bcba_glo, lmt_bcba)199 callscatter(lmt_omff_glo, lmt_omff)200 callscatter(lmt_omnff_glo, lmt_omnff)201 callscatter(lmt_ombbl_glo, lmt_ombbl)202 callscatter(lmt_ombbh_glo, lmt_ombbh)203 callscatter(lmt_omba_glo, lmt_omba)204 callscatter(lmt_terp_glo, lmt_terp)205 callscatter(lmt_omnat_glo, lmt_omnat)194 CALL scatter(lmt_bcff_glo, lmt_bcff) 195 CALL scatter(lmt_bcnff_glo, lmt_bcnff) 196 CALL scatter(lmt_bcbbl_glo, lmt_bcbbl) 197 CALL scatter(lmt_bcbbh_glo, lmt_bcbbh) 198 CALL scatter(lmt_bcba_glo, lmt_bcba) 199 CALL scatter(lmt_omff_glo, lmt_omff) 200 CALL scatter(lmt_omnff_glo, lmt_omnff) 201 CALL scatter(lmt_ombbl_glo, lmt_ombbl) 202 CALL scatter(lmt_ombbh_glo, lmt_ombbh) 203 CALL scatter(lmt_omba_glo, lmt_omba) 204 CALL scatter(lmt_terp_glo, lmt_terp) 205 CALL scatter(lmt_omnat_glo, lmt_omnat) 206 206 207 207 RETURN -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/condsurfs.f90
r5099 r5101 4 4 lmt_dmsbio, lmt_h2sbio, lmt_dms, lmt_dmsconc) 5 5 USE dimphy 6 USE netcdf, ONLY 6 USE netcdf, ONLY: nf90_close, nf90_noerr, nf90_inq_varid, nf90_open, & 7 7 nf90_nowrite, nf90_get_var 8 8 IMPLICIT none … … 39 39 write(6, *)' Pb d''ouverture du fichier limitsoufre.nc' 40 40 write(6, *)' ierr = ', ierr 41 callexit(1)41 CALL exit(1) 42 42 endif 43 43 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/condsurfs_new.f90
r5099 r5101 9 9 USE mod_phys_lmdz_para 10 10 USE dimphy 11 USE netcdf, ONLY 11 USE netcdf, ONLY: nf90_get_var, nf90_inq_varid, nf90_close, nf90_noerr, nf90_open, nf90_nowrite 12 12 IMPLICIT none 13 13 … … 64 64 write(6, *)' Pb d''ouverture du fichier sulphur_emissions_antro' 65 65 write(6, *)' ierr = ', ierr 66 callexit(1)66 CALL exit(1) 67 67 endif 68 68 … … 133 133 write(6, *)' Pb d''ouverture du fichier sulphur_emissions_nat' 134 134 write(6, *)' ierr = ', ierr 135 callexit(1)135 CALL exit(1) 136 136 endif 137 137 … … 188 188 write(6, *)' Pb d''ouverture du fichier sulphur_emissions_volc' 189 189 write(6, *)' ierr = ', ierr 190 callexit(1)190 CALL exit(1) 191 191 endif 192 192 … … 243 243 !$OMP END MASTER 244 244 !$OMP BARRIER 245 callscatter(lmt_so2b_glo, lmt_so2b)246 callscatter(lmt_so2h_glo, lmt_so2h)247 callscatter(lmt_so2bb_h_glo, lmt_so2bb_h)248 callscatter(lmt_so2bb_l_glo, lmt_so2bb_l)249 callscatter(lmt_so2ba_glo, lmt_so2ba)250 callscatter(lmt_so2nff_glo, lmt_so2nff)251 callscatter(lmt_dmsbio_glo, lmt_dmsbio)252 callscatter(lmt_h2sbio_glo, lmt_h2sbio)253 callscatter(lmt_dmsconc_glo, lmt_dmsconc)254 callscatter(lmt_dms_glo, lmt_dms)255 callscatter(lmt_so2volc_cont_glo, lmt_so2volc_cont)256 callscatter(lmt_altvolc_cont_glo, lmt_altvolc_cont)257 callscatter(lmt_so2volc_expl_glo, lmt_so2volc_expl)258 callscatter(lmt_altvolc_expl_glo, lmt_altvolc_expl)245 CALL scatter(lmt_so2b_glo, lmt_so2b) 246 CALL scatter(lmt_so2h_glo, lmt_so2h) 247 CALL scatter(lmt_so2bb_h_glo, lmt_so2bb_h) 248 CALL scatter(lmt_so2bb_l_glo, lmt_so2bb_l) 249 CALL scatter(lmt_so2ba_glo, lmt_so2ba) 250 CALL scatter(lmt_so2nff_glo, lmt_so2nff) 251 CALL scatter(lmt_dmsbio_glo, lmt_dmsbio) 252 CALL scatter(lmt_h2sbio_glo, lmt_h2sbio) 253 CALL scatter(lmt_dmsconc_glo, lmt_dmsconc) 254 CALL scatter(lmt_dms_glo, lmt_dms) 255 CALL scatter(lmt_so2volc_cont_glo, lmt_so2volc_cont) 256 CALL scatter(lmt_altvolc_cont_glo, lmt_altvolc_cont) 257 CALL scatter(lmt_so2volc_expl_glo, lmt_so2volc_expl) 258 CALL scatter(lmt_altvolc_expl_glo, lmt_altvolc_expl) 259 259 260 260 RETURN -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/dustemission_mod.F90
r5099 r5101 223 223 224 224 IF( debutphy ) THEN 225 ! callgather(maskdust,maskdust_glo)225 ! CALL gather(maskdust,maskdust_glo) 226 226 ! !$OMP MASTER 227 227 ! IF (is_mpi_root .AND. is_omp_root) THEN … … 349 349 350 350 if (sizedustmin>sizeacclow .or. sizedustmax<sizescohigh) then 351 callabort_gcm('adaptdustemission', 'Dust range problem',1)351 CALL abort_gcm('adaptdustemission', 'Dust range problem',1) 352 352 endif 353 353 print *,'FINE DUST BIN: tuning EMISSION factor= ',tuningfactorfine … … 393 393 ! allocate(itvmean_glo(klon_glo,nbins)) 394 394 ! ALLOCATE(emisbinlocalmean2_glo(klon_glo,nbins)) 395 ! callgather(itv2,itv2_glo)396 ! callgather(itvmean,itvmean_glo)395 ! CALL gather(itv2,itv2_glo) 396 ! CALL gather(itvmean,itvmean_glo) 397 397 !!$OMP MASTER 398 398 ! IF (is_mpi_root .AND. is_omp_root) THEN … … 486 486 enddo 487 487 counter1=0 488 ! callgather(emisbinlocalmean2,emisbinlocalmean2_glo)488 ! CALL gather(emisbinlocalmean2,emisbinlocalmean2_glo) 489 489 !!$OMP MASTER 490 490 ! IF (is_mpi_root .AND. is_omp_root) THEN … … 747 747 !$OMP BARRIER 748 748 !print *,'JEOK10',mpi_rank,omp_rank 749 callbcast(solspe)749 CALL bcast(solspe) 750 750 ! Calcul de la distribution en taille des particules de Dust 751 751 ! Elle depend du nombre de classe des particules nclass. … … 951 951 ! DO nb=1,nbinsHR 952 952 ! vdHR(nb)=vdout(nb) 953 !! WRITE(18,*) ,binsHR(nb),vdHR(nb)953 !! WRITE(18,*) binsHR(nb),vdHR(nb) 954 954 ! END DO 955 955 … … 1083 1083 DO k=1,ndistb 1084 1084 DO nb=1,nbins 1085 write(15001,*) ,k,nb,massfrac(k,nb)1085 write(15001,*) k,nb,massfrac(k,nb) 1086 1086 ENDDO 1087 1087 ENDDO … … 1384 1384 ! print*,'ustarsalt = ',ustarsalt,'dsmin=',dsmin,'dsmax=',dsmax 1385 1385 ! dichotomy 1386 calldef_dichotomy(sizeclass,nclass,1,ncl,dsmin,ideb)1386 CALL def_dichotomy(sizeclass,nclass,1,ncl,dsmin,ideb) 1387 1387 ! print*,'ideb = ',ideb 1388 calldef_dichotomy(sizeclass,nclass,ideb,ncl,dsmax,ifin)1388 CALL def_dichotomy(sizeclass,nclass,ideb,ncl,dsmax,ifin) 1389 1389 ! print*,'ifin = ',ifin 1390 1390 ! readaptation of large sizes particles -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/finemission.F
r5082 r5101 14 14 USE indice_sol_mod 15 15 USE infotrac 16 ! USE phytracr_spl_mod, ONLY 16 ! USE phytracr_spl_mod, ONLY: nbreg_dust, nbreg_ind, nbreg_bb 17 17 IMPLICIT NONE 18 18 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_orig.F90
r5099 r5101 9 9 USE mod_phys_lmdz_para 10 10 USE traclmdz_mod 11 USE infotrac,ONLY 11 USE infotrac,ONLY: nbtr 12 12 USE iophy 13 13 USE lmdz_yomcst … … 97 97 98 98 ! inscav_fisrt=.true. 99 ! callgetin('inscav_fisrt',inscav_fisrt)99 ! CALL getin('inscav_fisrt',inscav_fisrt) 100 100 ! if(inscav_fisrt) then 101 101 ! print*,'beta from fisrtilp.F90, beta = (z_cond - z_oliq)/z_cond, inscav_fisrt=',inscav_fisrt -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_spl.F90
r5099 r5101 11 11 USE mod_phys_lmdz_para 12 12 USE traclmdz_mod 13 USE infotrac,ONLY 13 USE infotrac,ONLY: nbtr 14 14 USE iophy 15 15 USE lmdz_yomcst … … 107 107 108 108 ! inscav_fisrt=.true. 109 ! callgetin('inscav_fisrt',inscav_fisrt)109 ! CALL getin('inscav_fisrt',inscav_fisrt) 110 110 ! if(inscav_fisrt) then 111 111 ! print*,'beta from fisrtilp.F90, beta = (z_cond - z_oliq)/z_cond, inscav_fisrt=',inscav_fisrt -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/neutral.F
r5082 r5101 52 52 f1 = (1. + phi_inv) / 2. 53 53 f2 = (1. + phi_inv_sq)/2. 54 c following to avoid numerical overruns. re calltan(90deg)=infinity54 c following to avoid numerical overruns. reCALL tan(90deg)=infinity 55 55 dum1 = min (1.e24, phi_inv) 56 56 f3 = atan(dum1) -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phys_output_write_spl_mod.F90
r5099 r5101 10 10 USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy 11 11 12 USE phytracr_spl_mod, ONLY 12 USE phytracr_spl_mod, ONLY: ok_chimeredust, id_prec, id_fine, id_coss, & 13 13 id_codu, id_scdu , & 14 14 d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, & … … 68 68 flux_sparam_sscoa,u10m_ss,v10m_ss 69 69 70 USE dustemission_mod, ONLY : m1dflux, m2dflux, m3dflux71 72 ! USE phytrac_mod, ONLY 70 USE dustemission_mod, ONLY: m1dflux, m2dflux, m3dflux 71 72 ! USE phytrac_mod, ONLY: d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, & 73 73 ! d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls, & 74 74 ! d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav … … 394 394 USE lmdz_xios, ONLY: xios_update_calendar, using_xios 395 395 USE wxios, ONLY: wxios_closedef, missing_val_xios => missing_val 396 USE phys_cal_mod, ONLY 396 USE phys_cal_mod, ONLY: mth_len 397 397 USE lmdz_yomcst 398 398 … … 466 466 IF (vars_defined) THEN 467 467 IF (prt_level >= 10) THEN 468 write(lunout,*)"phys_output_write: callxios_update_calendar, itau_w=",itau_w468 write(lunout,*)"phys_output_write: CALL xios_update_calendar, itau_w=",itau_w 469 469 ENDIF 470 470 ! CALL xios_update_calendar(itau_w) -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phytracr_spl_mod.F90
r5099 r5101 424 424 SUBROUTINE phytracr_spl_out_init() 425 425 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 426 !AS : This subroutine centralises the ALLOCATE needed for the 1st callof426 !AS : This subroutine centralises the ALLOCATE needed for the 1st CALL of 427 427 ! phys_output_write_spl in physiq 428 428 429 429 USE dimphy 430 USE infotrac_phy, ONLY 431 USE dustemission_mod, ONLY 430 USE infotrac_phy, ONLY: nbtr 431 USE dustemission_mod, ONLY: dustemis_out_init 432 432 433 433 ! pour les variables m[1-3]dflux … … 699 699 fileregionsdimsbb = 'regions_bb_meta' 700 700 fileregionsdimswstar = 'regions_pwstarwake_meta' 701 callreadregionsdims2_spl(nbreg_ind, fileregionsdimsind)702 callreadregionsdims2_spl(nbreg_dust, fileregionsdimsdust)703 callreadregionsdims2_spl(nbreg_bb, fileregionsdimsbb)704 callreadregionsdims2_spl(nbreg_wstardust, fileregionsdimswstar)701 CALL readregionsdims2_spl(nbreg_ind, fileregionsdimsind) 702 CALL readregionsdims2_spl(nbreg_dust, fileregionsdimsdust) 703 CALL readregionsdims2_spl(nbreg_bb, fileregionsdimsbb) 704 CALL readregionsdims2_spl(nbreg_wstardust, fileregionsdimswstar) 705 705 ENDIF ! ASSIM 706 706 ! fin debranchage … … 792 792 793 793 USE mod_phys_lmdz_transfert_para 794 USE lmdz_thermcell_dq, ONLY 795 USE phys_cal_mod, only: jD_1jan, year_len, mth_len, days_elapsed, jh_1jan, year_cur, &794 USE lmdz_thermcell_dq, ONLY: thermcell_dq 795 USE phys_cal_mod, ONLY: jD_1jan, year_len, mth_len, days_elapsed, jh_1jan, year_cur, & 796 796 mth_cur, phys_cal_update 797 797 … … 806 806 ! Remarques en vrac: 807 807 ! ------------------ 808 ! 1/ le callphytrac se fait avec nqmax-2 donc nous avons bien808 ! 1/ le CALL phytrac se fait avec nqmax-2 donc nous avons bien 809 809 ! les vrais traceurs (nbtr) dans phytrac (pas la vapeur ni eau liquide) 810 810 !! AS : nqmax-2 devrait etre nqmax-3 apres introducton de H2Oi ; … … 1231 1231 itr = itr+1 1232 1232 write(str2,'(i2.2)') itrr 1233 calliophys_ecrit('TRA'//str2,klev,'SOURCE','',tr_seri(:,:,itr))1233 CALL iophys_ecrit('TRA'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 1234 1234 enddo 1235 1235 #endif … … 1297 1297 !AS: La ligne suivante fait planter a l'execution : lmt_so2ff_l pas initialise 1298 1298 ! print *,'lmt_so2ff_l AVANT' , MINVAL(lmt_so2ff_l), MAXVAL(lmt_so2ff_l) 1299 callread_newemissions(ijulday, jH_cur, edgar, flag_dms, debutphy, & !I1299 CALL read_newemissions(ijulday, jH_cur, edgar, flag_dms, debutphy, & !I 1300 1300 pdtphys, lafin, nbjour, pctsrf, & !I 1301 1301 t_seri, rlat, rlon, & !I … … 1396 1396 if (ok_chimeredust) then 1397 1397 if (.not.(id_scdu>0 .and. id_codu>0 .and. id_fine>0)) then 1398 callabort_gcm('phytracr_mod', 'pb in ok_chimdust 0', 1)1398 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 0', 1) 1399 1399 endif 1400 1400 else 1401 1401 if (id_scdu>0) then 1402 callabort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU', 1)1402 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU', 1) 1403 1403 endif 1404 1404 if ((id_codu <= 0) .or. (id_fine<=0)) then 1405 callabort_gcm('phytracr_mod', 'pb in ok_chimdust 1', 1)1405 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 1', 1) 1406 1406 endif 1407 1407 endif … … 1660 1660 !Config 1661 1661 iflag_lscav_omp = 4 1662 callgetin('iflag_lscav', iflag_lscav_omp)1662 CALL getin('iflag_lscav', iflag_lscav_omp) 1663 1663 iflag_lscav = iflag_lscav_omp 1664 1664 ! initialiation for time computation … … 1826 1826 c_FullName1 = 'regions_dustacc' 1827 1827 !c_FullName1='regions_dust' 1828 callreadregions_spl(iregion_dust, c_FullName1)1828 CALL readregions_spl(iregion_dust, c_FullName1) 1829 1829 c_FullName1 = 'regions_ind' 1830 callreadregions_spl(iregion_ind, c_FullName1)1830 CALL readregions_spl(iregion_ind, c_FullName1) 1831 1831 c_FullName1 = 'regions_bb' 1832 callreadregions_spl(iregion_bb, c_FullName1)1832 CALL readregions_spl(iregion_bb, c_FullName1) 1833 1833 c_FullName1 = 'regions_pwstarwake' 1834 callreadregions_spl(iregion_wstardust, c_FullName1)1834 CALL readregions_spl(iregion_wstardust, c_FullName1) 1835 1835 1836 1836 !$OMP MASTER … … 1944 1944 aux_var2(i) = pctsrf(i, is_oce) 1945 1945 enddo 1946 callgather(aux_var2, auxklon_glo)1946 CALL gather(aux_var2, auxklon_glo) 1947 1947 !$OMP MASTER 1948 1948 IF (is_mpi_root .AND. is_omp_root) THEN … … 1957 1957 aux_var2(i) = pctsrf(i, is_sic) 1958 1958 enddo 1959 callgather(aux_var2, auxklon_glo)1959 CALL gather(aux_var2, auxklon_glo) 1960 1960 !$OMP MASTER 1961 1961 IF (is_mpi_root .AND. is_omp_root) THEN … … 1970 1970 aux_var2(i) = pctsrf(i, is_ter) 1971 1971 enddo 1972 callgather(aux_var2, auxklon_glo)1972 CALL gather(aux_var2, auxklon_glo) 1973 1973 !$OMP MASTER 1974 1974 IF (is_mpi_root .AND. is_omp_root) THEN … … 1983 1983 aux_var2(i) = pctsrf(i, is_lic) 1984 1984 enddo 1985 callgather(aux_var2, auxklon_glo)1985 CALL gather(aux_var2, auxklon_glo) 1986 1986 !$OMP MASTER 1987 1987 IF (is_mpi_root .AND. is_omp_root) THEN … … 2068 2068 do itr=1,nbtr 2069 2069 write(str2,'(i2.2)') itr 2070 calliophys_ecrit('sav'//str2,1,'SOURCE','',source_tr(:,itr))2071 calliophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,itr))2070 CALL iophys_ecrit('sav'//str2,1,'SOURCE','',source_tr(:,itr)) 2071 CALL iophys_ecrit('fav'//str2,1,'SOURCE','',source_tr(:,itr)) 2072 2072 enddo 2073 2073 do itr=1,nbtr 2074 2074 write(str2,'(i2.2)') itr 2075 calliophys_ecrit('TRB'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2075 CALL iophys_ecrit('TRB'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2076 2076 enddo 2077 2077 #endif … … 2101 2101 do itr=1,nbtr 2102 2102 write(str2,'(i2.2)') itr 2103 calliophys_ecrit('sap'//str2,1,'SOURCE','',source_tr(:,itr))2104 calliophys_ecrit('fap'//str2,1,'SOURCE','',source_tr(:,itr))2103 CALL iophys_ecrit('sap'//str2,1,'SOURCE','',source_tr(:,itr)) 2104 CALL iophys_ecrit('fap'//str2,1,'SOURCE','',source_tr(:,itr)) 2105 2105 enddo 2106 2106 #endif … … 2124 2124 #ifdef IOPHYS_DUST 2125 2125 print *,'INPUT TO PRECUREMISSION' 2126 calliophys_ecrit('ftsol',4,'ftsol','',ftsol)2127 calliophys_ecrit('u10m_ec',1,'u10m_ec','',u10m_ec)2128 calliophys_ecrit('v10m_ec',1,'v10m_ec','',v10m_ec)2129 calliophys_ecrit('pctsrf',4,'pctsrf','',pctsrf)2130 calliophys_ecrit('u_seri',klev,'u_seri','',u_seri)2131 calliophys_ecrit('v_seri',klev,'v_seri','',v_seri)2132 calliophys_ecrit('paprs',klev,'paprs','',paprs)2133 calliophys_ecrit('pplay',klev,'pplay','',pplay)2134 calliophys_ecrit('cdragh',1,'cdragh','',cdragh)2135 calliophys_ecrit('cdragm',1,'cdragm','',cdragm)2136 calliophys_ecrit('t_seri',klev,'t_seri','',t_seri)2137 calliophys_ecrit('q_seri',klev,'q_seri','',q_seri)2138 calliophys_ecrit('tsol',1,'tsol','',tsol)2126 CALL iophys_ecrit('ftsol',4,'ftsol','',ftsol) 2127 CALL iophys_ecrit('u10m_ec',1,'u10m_ec','',u10m_ec) 2128 CALL iophys_ecrit('v10m_ec',1,'v10m_ec','',v10m_ec) 2129 CALL iophys_ecrit('pctsrf',4,'pctsrf','',pctsrf) 2130 CALL iophys_ecrit('u_seri',klev,'u_seri','',u_seri) 2131 CALL iophys_ecrit('v_seri',klev,'v_seri','',v_seri) 2132 CALL iophys_ecrit('paprs',klev,'paprs','',paprs) 2133 CALL iophys_ecrit('pplay',klev,'pplay','',pplay) 2134 CALL iophys_ecrit('cdragh',1,'cdragh','',cdragh) 2135 CALL iophys_ecrit('cdragm',1,'cdragm','',cdragm) 2136 CALL iophys_ecrit('t_seri',klev,'t_seri','',t_seri) 2137 CALL iophys_ecrit('q_seri',klev,'q_seri','',q_seri) 2138 CALL iophys_ecrit('tsol',1,'tsol','',tsol) 2139 2139 print*,'fracso2emis,frach2sofso2,bateau',fracso2emis,frach2sofso2,bateau 2140 2140 print*,'kminbc,kmaxbc,pdtphys',kminbc,kmaxbc,pdtphys … … 2142 2142 print*,'iregion_ind,iregion_bb,nbreg_ind, nbreg_bb',iregion_ind,iregion_bb,nbreg_ind, nbreg_bb 2143 2143 print*,'id_prec,id_fine',id_prec,id_fine 2144 calliophys_ecrit('zdz',klev,'zdz','',zdz)2145 calliophys_ecrit('zalt',klev,'zalt','',zalt)2146 calliophys_ecrit('lmt_so2ff_l',1,'lmt_so2ff_l','',lmt_so2ff_l)2147 calliophys_ecrit('lmt_so2ff_h',1,'lmt_so2ff_h','',lmt_so2ff_h)2148 calliophys_ecrit('lmt_so2nff',1,'lmt_so2nff','',lmt_so2nff)2149 calliophys_ecrit('lmt_so2ba',1,'lmt_so2ba','',lmt_so2ba)2150 calliophys_ecrit('lmt_so2bb_l',1,'lmt_so2bb_l','',lmt_so2bb_l)2151 calliophys_ecrit('lmt_so2bb_h',1,'lmt_so2bb_h','',lmt_so2bb_h)2152 calliophys_ecrit('lmt_so2volc_cont',1,'lmt_so2volc_cont','',lmt_so2volc_cont)2153 calliophys_ecrit('lmt_altvolc_cont',1,'lmt_altvolc_cont','',lmt_altvolc_cont)2154 calliophys_ecrit('lmt_so2volc_expl',1,'lmt_so2volc_expl','',lmt_so2volc_expl)2155 calliophys_ecrit('lmt_altvolc_expl',1,'lmt_altvolc_expl','',lmt_altvolc_expl)2156 calliophys_ecrit('lmt_dmsbio',1,'lmt_dmsbio','',lmt_dmsbio)2157 calliophys_ecrit('lmt_h2sbio',1,'lmt_h2sbio','',lmt_h2sbio)2158 calliophys_ecrit('lmt_dmsconc',1,'lmt_dmsconc','',lmt_dmsconc)2159 calliophys_ecrit('lmt_dms',1,'lmt_dms','',lmt_dms)2160 calliophys_ecrit('flux_sparam_ind',1,'flux_sparam_ind','',flux_sparam_ind)2161 calliophys_ecrit('flux_sparam_bb',1,'flux_sparam_bb','',flux_sparam_bb)2144 CALL iophys_ecrit('zdz',klev,'zdz','',zdz) 2145 CALL iophys_ecrit('zalt',klev,'zalt','',zalt) 2146 CALL iophys_ecrit('lmt_so2ff_l',1,'lmt_so2ff_l','',lmt_so2ff_l) 2147 CALL iophys_ecrit('lmt_so2ff_h',1,'lmt_so2ff_h','',lmt_so2ff_h) 2148 CALL iophys_ecrit('lmt_so2nff',1,'lmt_so2nff','',lmt_so2nff) 2149 CALL iophys_ecrit('lmt_so2ba',1,'lmt_so2ba','',lmt_so2ba) 2150 CALL iophys_ecrit('lmt_so2bb_l',1,'lmt_so2bb_l','',lmt_so2bb_l) 2151 CALL iophys_ecrit('lmt_so2bb_h',1,'lmt_so2bb_h','',lmt_so2bb_h) 2152 CALL iophys_ecrit('lmt_so2volc_cont',1,'lmt_so2volc_cont','',lmt_so2volc_cont) 2153 CALL iophys_ecrit('lmt_altvolc_cont',1,'lmt_altvolc_cont','',lmt_altvolc_cont) 2154 CALL iophys_ecrit('lmt_so2volc_expl',1,'lmt_so2volc_expl','',lmt_so2volc_expl) 2155 CALL iophys_ecrit('lmt_altvolc_expl',1,'lmt_altvolc_expl','',lmt_altvolc_expl) 2156 CALL iophys_ecrit('lmt_dmsbio',1,'lmt_dmsbio','',lmt_dmsbio) 2157 CALL iophys_ecrit('lmt_h2sbio',1,'lmt_h2sbio','',lmt_h2sbio) 2158 CALL iophys_ecrit('lmt_dmsconc',1,'lmt_dmsconc','',lmt_dmsconc) 2159 CALL iophys_ecrit('lmt_dms',1,'lmt_dms','',lmt_dms) 2160 CALL iophys_ecrit('flux_sparam_ind',1,'flux_sparam_ind','',flux_sparam_ind) 2161 CALL iophys_ecrit('flux_sparam_bb',1,'flux_sparam_bb','',flux_sparam_bb) 2162 2162 #endif 2163 2163 … … 2198 2198 do itr=1,nbtr 2199 2199 write(str2,'(i2.2)') itr 2200 calliophys_ecrit('tpr'//str2,1,'SOURCE','',source_tr(:,itr))2201 calliophys_ecrit('fpr'//str2,1,'SOURCE','',flux_tr(:,itr))2200 CALL iophys_ecrit('tpr'//str2,1,'SOURCE','',source_tr(:,itr)) 2201 CALL iophys_ecrit('fpr'//str2,1,'SOURCE','',flux_tr(:,itr)) 2202 2202 enddo 2203 2203 #endif … … 2242 2242 do itr=1,nbtr 2243 2243 write(str2,'(i2.2)') itr 2244 calliophys_ecrit('t'//str2,1,'SOURCE','',source_tr(:,itr))2245 calliophys_ecrit('f'//str2,1,'SOURCE','',flux_tr(:,itr))2244 CALL iophys_ecrit('t'//str2,1,'SOURCE','',source_tr(:,itr)) 2245 CALL iophys_ecrit('f'//str2,1,'SOURCE','',flux_tr(:,itr)) 2246 2246 enddo 2247 2247 #endif … … 2298 2298 do itr=1,nbtr 2299 2299 write(str2,'(i2.2)') itr 2300 calliophys_ecrit('TRC'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2300 CALL iophys_ecrit('TRC'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2301 2301 enddo 2302 2302 #endif … … 2338 2338 do itr=1,nbtr 2339 2339 write(str2,'(i2.2)') itr 2340 calliophys_ecrit('TRD'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2340 CALL iophys_ecrit('TRD'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2341 2341 enddo 2342 2342 #endif … … 2406 2406 #ifdef IOPHYS_DUST 2407 2407 print*,'iflag_conv=',iflag_conv 2408 calliophys_ecrit('coefh',klev,'coefh','',coefh)2409 calliophys_ecrit('yu1',1,'yu1','',yu1)2410 calliophys_ecrit('yv1',1,'yv1','',yv1)2411 calliophys_ecrit('delp',klev,'delp','',delp)2408 CALL iophys_ecrit('coefh',klev,'coefh','',coefh) 2409 CALL iophys_ecrit('yu1',1,'yu1','',yu1) 2410 CALL iophys_ecrit('yv1',1,'yv1','',yv1) 2411 CALL iophys_ecrit('delp',klev,'delp','',delp) 2412 2412 do itr=1,nbtr 2413 2413 write(str2,'(i2.2)') itr 2414 calliophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2414 CALL iophys_ecrit('TRE'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2415 2415 enddo 2416 2416 #endif … … 2604 2604 do itr=1,nbtr 2605 2605 write(str2,'(i2.2)') itr 2606 calliophys_ecrit('TRF'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2606 CALL iophys_ecrit('TRF'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2607 2607 enddo 2608 2608 #endif … … 2691 2691 do itr=1,nbtr 2692 2692 write(str2,'(i2.2)') itr 2693 calliophys_ecrit('TRG'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2693 CALL iophys_ecrit('TRG'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2694 2694 enddo 2695 2695 #endif … … 2856 2856 do itr=1,nbtr 2857 2857 write(str2,'(i2.2)') itr 2858 calliophys_ecrit('TRH'//str2,klev,'SOURCE','',tr_seri(:,:,itr))2858 CALL iophys_ecrit('TRH'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 2859 2859 enddo 2860 2860 #endif … … 3033 3033 3034 3034 #ifdef IOPHYS_DUST 3035 calliophys_ecrit('da',klev,'da','',da)3036 calliophys_ecrit('phi',klev,'phi','',phi)3037 calliophys_ecrit('phi2',klev,'phi2','',phi2)3038 calliophys_ecrit('d1a',klev,'d1a','',d1a)3039 calliophys_ecrit('dam',klev,'dam','',dam)3040 calliophys_ecrit('mp',klev,'mp','',mp)3041 calliophys_ecrit('ep',klev,'ep','',ep)3042 calliophys_ecrit('sigd',klev,'sigd','',sigd)3043 calliophys_ecrit('sij',klev,'sij','',sij)3044 calliophys_ecrit('wght_cvfd',klev,'wght_cvfd','',wght_cvfd)3045 calliophys_ecrit('clw',klev,'clw','',clw)3046 calliophys_ecrit('elij',klev,'elij','',elij)3047 calliophys_ecrit('epmlmMm',klev,'epmlmMm','',epmlmMm)3048 calliophys_ecrit('eplaMm',klev,'eplaMm','',eplaMm)3049 calliophys_ecrit('pmflxr',klev,'pmflxr','',pmflxr)3050 calliophys_ecrit('pmflxs',klev,'pmflxs','',pmflxs)3051 calliophys_ecrit('evapls',klev,'evapls','',evapls)3052 calliophys_ecrit('wdtrainA',klev,'wdtrainA','',wdtrainA)3053 calliophys_ecrit('wdtrainM',klev,'wdtrainM','',wdtrainM)3035 CALL iophys_ecrit('da',klev,'da','',da) 3036 CALL iophys_ecrit('phi',klev,'phi','',phi) 3037 CALL iophys_ecrit('phi2',klev,'phi2','',phi2) 3038 CALL iophys_ecrit('d1a',klev,'d1a','',d1a) 3039 CALL iophys_ecrit('dam',klev,'dam','',dam) 3040 CALL iophys_ecrit('mp',klev,'mp','',mp) 3041 CALL iophys_ecrit('ep',klev,'ep','',ep) 3042 CALL iophys_ecrit('sigd',klev,'sigd','',sigd) 3043 CALL iophys_ecrit('sij',klev,'sij','',sij) 3044 CALL iophys_ecrit('wght_cvfd',klev,'wght_cvfd','',wght_cvfd) 3045 CALL iophys_ecrit('clw',klev,'clw','',clw) 3046 CALL iophys_ecrit('elij',klev,'elij','',elij) 3047 CALL iophys_ecrit('epmlmMm',klev,'epmlmMm','',epmlmMm) 3048 CALL iophys_ecrit('eplaMm',klev,'eplaMm','',eplaMm) 3049 CALL iophys_ecrit('pmflxr',klev,'pmflxr','',pmflxr) 3050 CALL iophys_ecrit('pmflxs',klev,'pmflxs','',pmflxs) 3051 CALL iophys_ecrit('evapls',klev,'evapls','',evapls) 3052 CALL iophys_ecrit('wdtrainA',klev,'wdtrainA','',wdtrainA) 3053 CALL iophys_ecrit('wdtrainM',klev,'wdtrainM','',wdtrainM) 3054 3054 3055 3055 do itr=1,nbtr 3056 3056 write(str2,'(i2.2)') itr 3057 calliophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,itr))3057 CALL iophys_ecrit('TRI'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 3058 3058 enddo 3059 3059 #endif … … 3176 3176 do itr=1,nbtr 3177 3177 write(str2,'(i2.2)') itr 3178 calliophys_ecrit('TRJ'//str2,klev,'SOURCE','',tr_seri(:,:,itr))3178 CALL iophys_ecrit('TRJ'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 3179 3179 enddo 3180 3180 #endif … … 3442 3442 do itr=1,nbtr 3443 3443 write(str2,'(i2.2)') itr 3444 calliophys_ecrit('TRK'//str2,klev,'SOURCE','',tr_seri(:,:,itr))3444 CALL iophys_ecrit('TRK'//str2,klev,'SOURCE','',tr_seri(:,:,itr)) 3445 3445 enddo 3446 3446 #endif … … 4384 4384 IF (test_sca == 0) THEN 4385 4385 ! READ file!! 4386 callread_scalenc(filescaleparams, paramname_ind, &4386 CALL read_scalenc(filescaleparams, paramname_ind, & 4387 4387 nbreg_ind, step_sca, & 4388 4388 scale_param_ind) 4389 callread_scalenc(filescaleparams, paramname_bb, &4389 CALL read_scalenc(filescaleparams, paramname_bb, & 4390 4390 nbreg_bb, step_sca, & 4391 4391 scale_param_bb) 4392 callread_scalenc(filescaleparams, paramname_ff, &4392 CALL read_scalenc(filescaleparams, paramname_ff, & 4393 4393 nbreg_ff, step_sca, & 4394 4394 scale_param_ff) 4395 callread_scalenc(filescaleparams, paramname_dustacc, &4395 CALL read_scalenc(filescaleparams, paramname_dustacc, & 4396 4396 nbreg_dustacc, step_sca, & 4397 4397 scale_param_dustacc) 4398 callread_scalenc(filescaleparams, paramname_dustcoa, &4398 CALL read_scalenc(filescaleparams, paramname_dustcoa, & 4399 4399 nbreg_dustcoa, step_sca, & 4400 4400 scale_param_dustcoa) 4401 callread_scalenc(filescaleparams, paramname_dustsco, &4401 CALL read_scalenc(filescaleparams, paramname_dustsco, & 4402 4402 nbreg_dustsco, step_sca, & 4403 4403 scale_param_dustsco) 4404 callread_scalenc(filescaleparams, paramname_wstarBL, &4404 CALL read_scalenc(filescaleparams, paramname_wstarBL, & 4405 4405 nbreg_wstardustBL, step_sca, & 4406 4406 param_wstarBLperregion) 4407 callread_scalenc(filescaleparams, paramname_wstarWAKE, &4407 CALL read_scalenc(filescaleparams, paramname_wstarWAKE, & 4408 4408 nbreg_wstardustWAKE, step_sca, & 4409 4409 param_wstarWAKEperregion) 4410 callread_scalenc(filescaleparams, paramname_ssacc, &4410 CALL read_scalenc(filescaleparams, paramname_ssacc, & 4411 4411 nbreg_ssacc, step_sca, & 4412 4412 scale_param_ssacc_tmp) 4413 callread_scalenc(filescaleparams, paramname_sscoa, &4413 CALL read_scalenc(filescaleparams, paramname_sscoa, & 4414 4414 nbreg_sscoa, step_sca, & 4415 4415 scale_param_sscoa_tmp) … … 4435 4435 USE mod_grid_phy_lmdz 4436 4436 USE mod_phys_lmdz_para 4437 USE netcdf, ONLY 4437 USE netcdf, ONLY: nf90_open, nf90_close, nf90_inq_varid, nf90_nowrite, nf90_noerr, nf90_get_var 4438 4438 IMPLICIT NONE 4439 4439 … … 4466 4466 print *, 'error ierr= ', ierr 4467 4467 CALL exit(1) 4468 callabort_gcm('read_scalenc', 'error reading variable', 1)4468 CALL abort_gcm('read_scalenc', 'error reading variable', 1) 4469 4469 ENDIF 4470 4470 … … 4483 4483 !$OMP BARRIER 4484 4484 ! CALL scatter(var local _glo,var local) o algo asi 4485 callbcast(scale_param)4485 CALL bcast(scale_param) 4486 4486 END SUBROUTINE read_scalenc 4487 4487 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/precuremission.F
r5082 r5101 19 19 USE indice_sol_mod 20 20 USE infotrac 21 ! USE phytracr_spl_mod, ONLY 21 ! USE phytracr_spl_mod, ONLY: nbreg_dust, nbreg_ind, nbreg_bb 22 22 IMPLICIT NONE 23 23 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_dust.f90
r5099 r5101 3 3 USE mod_grid_phy_lmdz 4 4 USE mod_phys_lmdz_para 5 USE netcdf, ONLY 5 USE netcdf, ONLY: nf90_get_var, nf90_nowrite, nf90_open, nf90_inq_varid 6 6 IMPLICIT NONE 7 7 … … 44 44 status = nf90_get_var(ncid1, varid1, dust_nc_glo, start, count) 45 45 46 ! callcorrectbid(iim,jjp1,dust_nc)47 callcorrectbid(nbp_lon, nbp_lat, dust_nc_glo)46 ! CALL correctbid(iim,jjp1,dust_nc) 47 CALL correctbid(nbp_lon, nbp_lat, dust_nc_glo) 48 48 49 49 !--upside down + physical grid -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_newemissions.F
r5099 r5101 181 181 !$OMP END MASTER 182 182 !$OMP BARRIER 183 callscatter(wth_glo,wth)184 callscatter(cly_glo,cly)185 callscatter(zprecipinsoil_glo,zprecipinsoil)183 CALL scatter(wth_glo,wth) 184 CALL scatter(cly_glo,cly) 185 CALL scatter(zprecipinsoil_glo,zprecipinsoil) 186 186 187 187 !JE20140908<< GOTO 1000 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_surface.F90
r5099 r5101 106 106 status=nf90_get_var(ncid,varid,tmp_dyn_glo,start,count) 107 107 108 ! calldump2d(iip1,jjp1,tmp_dyn,'tmp_dyn ')108 ! CALL dump2d(iip1,jjp1,tmp_dyn,'tmp_dyn ') 109 109 DO j=1, nbp_lat 110 110 DO ig=1, nbp_lon+1 … … 114 114 115 115 116 !JE20140522! callgr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi)116 !JE20140522! CALL gr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) 117 117 118 118 !JE20140526<< 119 ! callgr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi)119 ! CALL gr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) 120 120 if (isinversed) then 121 callgr_dyn_fi(1, nbp_lon+1, nbp_lat, klon_glo, &121 CALL gr_dyn_fi(1, nbp_lon+1, nbp_lat, klon_glo, & 122 122 tmp_dyn_invers_glo, tmp_fi_glo) 123 ! callgr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi)124 ! callgr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi)123 ! CALL gr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) 124 ! CALL gr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) 125 125 else 126 callgr_dyn_fi(1, nbp_lon+1, nbp_lat, klon_glo, &126 CALL gr_dyn_fi(1, nbp_lon+1, nbp_lat, klon_glo, & 127 127 tmp_dyn_glo, tmp_fi_glo) 128 ! callgr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn, tmp_fi)129 ! callgr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn, tmp_fi)128 ! CALL gr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn, tmp_fi) 129 ! CALL gr_dyn_fi_p(1, iip1, jjp1, klon, tmp_dyn, tmp_fi) 130 130 endif 131 131 !JE20140526>> 132 ! calldump2d(iim,jjm-1,tmp_fi(2),'tmp_fi ')132 ! CALL dump2d(iim,jjm-1,tmp_fi(2),'tmp_fi ') 133 133 134 134 DO j=1,klon_glo … … 145 145 !$OMP END MASTER 146 146 !$OMP BARRIER 147 callscatter(surfa_glo,surfa)147 CALL scatter(surfa_glo,surfa) 148 148 149 149 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_vent.f90
r5099 r5101 3 3 USE mod_grid_phy_lmdz 4 4 USE mod_phys_lmdz_para 5 USE netcdf, ONLY 5 USE netcdf, ONLY: nf90_get_var, nf90_open, nf90_inq_varid, nf90_nowrite 6 6 IMPLICIT NONE 7 7 INCLUDE "dimensions.h" … … 57 57 58 58 ! print *,status 59 ! callcorrectbid(iim,jjp1,u10m_nc)60 ! callcorrectbid(iim,jjp1,v10m_nc)61 callcorrectbid(nbp_lon, nbp_lat, u10m_nc_glo)62 callcorrectbid(nbp_lon, nbp_lat, v10m_nc_glo)59 ! CALL correctbid(iim,jjp1,u10m_nc) 60 ! CALL correctbid(iim,jjp1,v10m_nc) 61 CALL correctbid(nbp_lon, nbp_lat, u10m_nc_glo) 62 CALL correctbid(nbp_lon, nbp_lat, v10m_nc_glo) 63 63 64 64 ! print *,'afterbidcor u10m_nc', u10m_nc(1,jjp1) -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/splaeropt_5wv_rrtm.f90
r5099 r5101 8 8 USE DIMPHY 9 9 USE aero_mod 10 USE infotrac_phy, ONLY 11 USE phys_local_var_mod, ONLY 10 USE infotrac_phy, ONLY: nqtot, nbtr, tracers 11 USE phys_local_var_mod, ONLY: od550aer, od865aer, ec550aer, od550lt1aer 12 12 13 13 ! Olivier Boucher Jan 2017 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/splaeropt_6bands_rrtm.f90
r5099 r5101 8 8 USE dimphy 9 9 USE aero_mod 10 USE infotrac_phy, ONLY 11 USE phys_local_var_mod, ONLY 10 USE infotrac_phy, ONLY: nqtot, nbtr, tracers 11 USE phys_local_var_mod, ONLY: abs550aer 12 12 13 13 ! Olivier Boucher Jan 2017 -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/splaeropt_lw_rrtm.f90
r5099 r5101 10 10 USE dimphy 11 11 USE aero_mod 12 USE infotrac_phy, ONLY 13 USE phys_state_var_mod, ONLY 14 USE lmdz_yoerad, ONLY 12 USE infotrac_phy, ONLY: nqtot, nbtr, tracers 13 USE phys_state_var_mod, ONLY: tau_aero_lw_rrtm 14 USE lmdz_yoerad, ONLY: NLW 15 15 16 16 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/splaerosol_optic_rrtm.f90
r5099 r5101 12 12 USE dimphy 13 13 USE aero_mod 14 USE infotrac_phy, ONLY 15 USE lmdz_yomcst, ONLY 14 USE infotrac_phy, ONLY: nbtr, nqtot, tracers 15 USE lmdz_yomcst, ONLY: RD, RG 16 16 17 17 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90
r5086 r5101 108 108 if (Jwarm) then 109 109 if (rain_effect) then 110 callNear_Surface(al, t_subskin, s_subskin, ds_ns, dt_ns, &110 CALL Near_Surface(al, t_subskin, s_subskin, ds_ns, dt_ns, & 111 111 tau_with_min, taur, hlb, rhoa, xlv, dtime, t_ocean_1, s1, rain, & 112 112 q_pwp = fxp * rns - (hf + hlb + rnl + rf)) 113 113 else 114 callNear_Surface(al, t_subskin, s_subskin, ds_ns, dt_ns, &114 CALL Near_Surface(al, t_subskin, s_subskin, ds_ns, dt_ns, & 115 115 tau_with_min, taur, hlb, rhoa, xlv, dtime, t_ocean_1, s1, & 116 116 rain = null_array, q_pwp = fxp * rns - (hf + hlb + rnl)) … … 132 132 * (1. - exp(- tkt / 8e-4))) ! equation 16 Ohlmann 133 133 if (rain_effect) then 134 callMicrolayer(dter, dser, tkt, tks, hlb, tau_with_min, &134 CALL Microlayer(dter, dser, tkt, tks, hlb, tau_with_min, & 135 135 s_subskin, al, xlv, taur, rf, rain, & 136 136 qcol = rnl + hf + hlb - dels) 137 137 else 138 callMicrolayer(dter, dser, tkt, tks, hlb, tau_with_min, &138 CALL Microlayer(dter, dser, tkt, tks, hlb, tau_with_min, & 139 139 s_subskin, al, xlv, taur, rf = null_array, & 140 140 rain = null_array, qcol = rnl + hf + hlb - dels) -
LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/config_ocean_skin_m.F90
r4176 r5101 42 42 43 43 #ifdef IN_LMDZ 44 callgetin_p("activate_ocean_skin", activate_ocean_skin)45 callassert(activate_ocean_skin >= 0 .and. activate_ocean_skin <= 2, &44 CALL getin_p("activate_ocean_skin", activate_ocean_skin) 45 CALL assert(activate_ocean_skin >= 0 .and. activate_ocean_skin <= 2, & 46 46 "config_ocean_skin bad value of activate_ocean_skin") 47 47 if (activate_ocean_skin >= 1) then 48 callgetin_p("flag_ocean_skin", flag_ocean_skin)49 callgetin_p("depth_1", depth_1)48 CALL getin_p("flag_ocean_skin", flag_ocean_skin) 49 CALL getin_p("depth_1", depth_1) 50 50 end if 51 51 #else -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/aer_sedimnt.F90
r5099 r5101 19 19 USE phys_local_var_mod, ONLY: mdw, budg_sed_part, DENSO4, DENSO4B, f_r_wet, f_r_wetB, vsed_aer 20 20 USE strataer_local_var_mod, ONLY: flag_new_strat_compo 21 USE dimphy, ONLY 21 USE dimphy, ONLY: klon,klev 22 22 USE infotrac_phy 23 23 USE aerophys -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/coagulate.F90
r5099 r5101 23 23 ! ----------------------------------------------------------------------- 24 24 25 USE dimphy, ONLY 25 USE dimphy, ONLY: klon,klev 26 26 USE aerophys 27 27 USE infotrac_phy -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
r5099 r5101 28 28 USE aerophys 29 29 USE infotrac_phy 30 USE lmdz_yomcst, ONLY 31 USE sulfate_aer_mod, ONLY 32 USE strataer_local_var_mod, ONLY 30 USE lmdz_yomcst, ONLY: RPI 31 USE sulfate_aer_mod, ONLY: wph2so4, surftension, solh2so4, rpmvh2so4 32 USE strataer_local_var_mod, ONLY: ALPH2SO4, RRSI 33 33 34 34 IMPLICIT NONE … … 177 177 USE aerophys 178 178 USE infotrac_phy 179 USE lmdz_yomcst, ONLY 180 USE strataer_local_var_mod, ONLY 179 USE lmdz_yomcst, ONLY: RPI 180 USE strataer_local_var_mod, ONLY: ALPH2SO4, RRSI 181 181 182 182 IMPLICIT NONE … … 290 290 USE aerophys 291 291 USE infotrac_phy 292 USE lmdz_yomcst, ONLY 293 USE strataer_local_var_mod, ONLY 292 USE lmdz_yomcst, ONLY: RPI 293 USE strataer_local_var_mod, ONLY: RRSI,Vbin 294 294 295 295 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/interp_sulf_input.F90
r5099 r5101 9 9 10 10 USE mod_grid_phy_lmdz 11 USE mod_phys_lmdz_mpi_data, ONLY 12 USE mod_phys_lmdz_omp_data, ONLY 13 USE phys_local_var_mod, ONLY 14 USE phys_local_var_mod, ONLY 11 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 12 USE mod_phys_lmdz_omp_data, ONLY: is_omp_root 13 USE phys_local_var_mod, ONLY: budg_3D_backgr_ocs, budg_3D_backgr_so2 14 USE phys_local_var_mod, ONLY: OCS_lifetime, SO2_lifetime, H2SO4_lifetime, O3_clim 15 15 USE mod_phys_lmdz_para 16 16 USE dimphy … … 19 19 USE aerophys 20 20 USE lmdz_yomcst 21 USE strataer_local_var_mod, ONLY 21 USE strataer_local_var_mod, ONLY: flag_newclim_file,flag_verbose_strataer 22 22 23 23 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/micphy_tstep.F90
r5099 r5101 4 4 SUBROUTINE micphy_tstep(pdtphys,tr_seri,t_seri,pplay,paprs,rh,is_strato) 5 5 6 USE geometry_mod, ONLY 7 USE dimphy, ONLY 6 USE geometry_mod, ONLY: latitude_deg !NL- latitude corr. to local domain 7 USE dimphy, ONLY: klon,klev 8 8 USE aerophys 9 USE infotrac_phy, ONLY 9 USE infotrac_phy, ONLY: nbtr_bin, nbtr_sulgas, nbtr, id_H2SO4_strat 10 10 USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, budg_h2so4_to_part, R2SO4, DENSO4, & 11 11 f_r_wet, R2SO4B, DENSO4B, f_r_wetB 12 12 USE nucleation_tstep_mod 13 13 USE cond_evap_tstep_mod 14 USE sulfate_aer_mod, ONLY 15 USE lmdz_yomcst, ONLY 14 USE sulfate_aer_mod, ONLY: STRAACT 15 USE lmdz_yomcst, ONLY: RPI, RD, RG 16 16 USE print_control_mod, ONLY: lunout 17 17 USE strataer_local_var_mod ! contains also RRSI and Vbin -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/miecalc_aer.F90
r5099 r5101 18 18 USE aerophys, ONLY: dens_aer_dry, dens_aer_ref, V_rat 19 19 USE aero_mod 20 USE infotrac_phy, ONLY 20 USE infotrac_phy, ONLY: nbtr, nbtr_bin, nbtr_sulgas, id_SO2_strat 21 21 USE dimphy 22 USE lmdz_yomcst , ONLY 22 USE lmdz_yomcst , ONLY: RG, RPI 23 23 USE mod_phys_lmdz_para, only: gather, scatter, bcast 24 USE mod_grid_phy_lmdz, ONLY 25 USE mod_phys_lmdz_mpi_data, ONLY 24 USE mod_grid_phy_lmdz, ONLY: klon_glo 25 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 26 26 USE print_control_mod, ONLY: prt_level, lunout 27 27 -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/nucleation_tstep_mod.F90
r5099 r5101 10 10 USE aerophys 11 11 USE infotrac_phy 12 USE strataer_local_var_mod, ONLY 12 USE strataer_local_var_mod, ONLY: flag_new_nucl 13 13 USE lmdz_yomcst 14 14 … … 39 39 REAL ipr ! Ion pair production rate (cm-3 s-1) NOT IN USE 40 40 41 ! callnucleation routine41 ! CALL nucleation routine 42 42 IF (.NOT.flag_new_nucl) THEN 43 43 ! Use older routine from Hanna Vehkamäki (FMI) … … 73 73 USE aerophys 74 74 USE infotrac_phy 75 USE strataer_local_var_mod, ONLY 75 USE strataer_local_var_mod, ONLY: Vbin 76 76 77 77 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/ocs_to_so2.F90
r5099 r5101 4 4 SUBROUTINE ocs_to_so2(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) 5 5 6 USE dimphy, ONLY 6 USE dimphy, ONLY: klon,klev 7 7 USE aerophys 8 8 USE infotrac_phy 9 USE lmdz_yomcst, ONLY 10 USE phys_local_var_mod, ONLY 11 USE strataer_local_var_mod, ONLY 9 USE lmdz_yomcst, ONLY: RG 10 USE phys_local_var_mod, ONLY: OCS_lifetime, budg_3D_ocs_to_so2, budg_ocs_to_so2 11 USE strataer_local_var_mod, ONLY: flag_min_rreduce 12 12 13 13 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/so2_to_h2so4.F90
r5099 r5101 4 4 SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) 5 5 6 USE dimphy, ONLY 6 USE dimphy, ONLY: klon,klev 7 7 USE aerophys 8 8 USE infotrac_phy 9 USE lmdz_yomcst, ONLY 9 USE lmdz_yomcst, ONLY: RG, RD 10 10 ! lifetime (sec) et O3_clim (VMR) 11 USE phys_local_var_mod, ONLY 12 USE strataer_local_var_mod, ONLY 11 USE phys_local_var_mod, ONLY: SO2_lifetime, H2SO4_lifetime, O3_clim, budg_3D_so2_to_h2so4, budg_so2_to_h2so4 12 USE strataer_local_var_mod, ONLY: flag_OH_reduced, flag_H2SO4_photolysis, flag_min_rreduce 13 13 14 14 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratH2O_methox.F90
r5099 r5101 11 11 12 12 USE mod_grid_phy_lmdz 13 USE mod_phys_lmdz_mpi_data, ONLY 14 USE mod_phys_lmdz_omp_data, ONLY 13 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 14 USE mod_phys_lmdz_omp_data, ONLY: is_omp_root 15 15 16 16 USE mod_phys_lmdz_para 17 17 USE dimphy 18 USE phys_cal_mod, ONLY 18 USE phys_cal_mod, ONLY: mth_cur 19 19 USE infotrac_phy 20 20 USE aerophys 21 21 USE lmdz_yomcst 22 USE strataer_local_var_mod, ONLY 22 USE strataer_local_var_mod, ONLY: flag_newclim_file 23 23 24 24 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/strataer_emiss_mod.F90
r5098 r5101 7 7 8 8 USE strataer_local_var_mod 9 USE ioipsl_getin_p_mod, ONLY 10 USE print_control_mod, ONLY 11 USE mod_phys_lmdz_para, ONLY 9 USE ioipsl_getin_p_mod, ONLY: getin_p 10 USE print_control_mod, ONLY: lunout 11 USE mod_phys_lmdz_para, ONLY: is_master 12 12 13 13 ! Local variables … … 217 217 USE dimphy, ONLY: klon 218 218 USE mod_grid_phy_lmdz, ONLY: nbp_lat, nbp_lon 219 USE print_control_mod, ONLY 219 USE print_control_mod, ONLY: lunout 220 220 USE strataer_local_var_mod 221 221 USE lmdz_yomcst, ONLY: RPI -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/strataer_nuc_mod.F90
r5099 r5101 10 10 SUBROUTINE strataer_nuc_init() 11 11 12 USE ioipsl_getin_p_mod, ONLY 13 USE print_control_mod, ONLY 14 USE mod_phys_lmdz_para, ONLY 12 USE ioipsl_getin_p_mod, ONLY: getin_p 13 USE print_control_mod, ONLY: lunout 14 USE mod_phys_lmdz_para, ONLY: is_master 15 15 USE strataer_local_var_mod, ONLY: ALPH2SO4,flag_nuc_rate_box,nuclat_min,nuclat_max, & 16 16 nucpres_min,nucpres_max … … 50 50 51 51 USE infotrac_phy, ONLY: nbtr, nbtr_sulgas, id_H2SO4_strat 52 USE ioipsl, ONLY 53 USE print_control_mod, ONLY 52 USE ioipsl, ONLY: getin 53 USE print_control_mod, ONLY: lunout 54 54 55 55 ! Output variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratdistrib.F90
r5099 r5101 4 4 SUBROUTINE STRATDISTRIB(altLMDz,altemiss,sigma_alt,f_lay_emiss) 5 5 6 USE dimphy, ONLY 6 USE dimphy, ONLY: klon,klev 7 7 USE strataer_local_var_mod 8 USE lmdz_yomcst, only: RPI8 USE lmdz_yomcst, ONLY: RPI 9 9 10 10 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/stratemit.F90
r5099 r5101 6 6 stretchlong,ispeci,id_species_total) 7 7 8 USE dimphy, ONLY 8 USE dimphy, ONLY: klon,klev 9 9 USE strataer_local_var_mod 10 10 USE phys_cal_mod 11 11 USE phys_local_var_mod, ONLY: d_q_emiss 12 USE lmdz_yomcst, only: RD, RPI, RG13 USE geometry_mod, ONLY 12 USE lmdz_yomcst, ONLY: RD, RPI, RG 13 USE geometry_mod, ONLY: cell_area, boundslat 14 14 USE aerophys 15 15 USE infotrac_phy 16 USE mod_phys_lmdz_mpi_data, ONLY 16 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 17 17 USE mod_grid_phy_lmdz, ONLY: nbp_lon 18 18 -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90
r5099 r5101 25 25 ! assuming 'curved surface' composition (depends on aerosol size) 26 26 27 USE dimphy, ONLY 28 USE infotrac_phy, ONLY 27 USE dimphy, ONLY: klon,klev ! nb of longitude and altitude bands 28 USE infotrac_phy, ONLY: nbtr_bin 29 29 USE aerophys 30 30 USE phys_local_var_mod, ONLY: R2SO4, R2SO4B, DENSO4, DENSO4B, f_r_wet, f_r_wetB … … 164 164 ! R2SO4: aerosol H2SO4 weight fraction (percent) 165 165 166 USE dimphy, ONLY 166 USE dimphy, ONLY: klon,klev 167 167 USE aerophys 168 168 USE phys_local_var_mod, ONLY: R2SO4 … … 550 550 ! ACTSO4: H2SO4 activity (percent) 551 551 552 USE dimphy, ONLY 552 USE dimphy, ONLY: klon,klev 553 553 USE phys_local_var_mod, ONLY: R2SO4 554 554 … … 635 635 ! DENSO4: aerosol mass density (gr/cm3 = aerosol mass/aerosol volume) 636 636 637 USE dimphy, ONLY 637 USE dimphy, ONLY: klon,klev 638 638 USE phys_local_var_mod, ONLY: R2SO4, DENSO4 639 639 … … 986 986 ! mfh2so4 = Mh2so4*sulfmolal / (1000.+Mh2so4*sulfmolal) 987 987 ! wph2so4 (% mass fraction)= 100.*Mh2so4*sulfmolal / (1000.+Mh2so4*sulfmolal) 988 ! re call activity of i = a_i = P_i/P_pure_i and988 ! reCALL activity of i = a_i = P_i/P_pure_i and 989 989 ! activity coefficient of i = gamma_i = a_i/X_i (X_i: mole fraction of i) 990 990 ! so P_i = gamma_i*X_i*P_pure_i -
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/traccoag_mod.F90
r5099 r5101 15 15 16 16 USE dimphy 17 USE infotrac_phy, ONLY 17 USE infotrac_phy, ONLY: nbtr_bin, nbtr_sulgas, nbtr, id_SO2_strat 18 18 USE aerophys 19 USE geometry_mod, ONLY 19 USE geometry_mod, ONLY: cell_area, boundslat 20 20 USE mod_grid_phy_lmdz 21 USE mod_phys_lmdz_mpi_data, ONLY 21 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 22 22 USE mod_phys_lmdz_para, only: gather, scatter 23 USE phys_cal_mod, ONLY 23 USE phys_cal_mod, ONLY: year_len, year_cur, mth_cur, day_cur, hour 24 24 USE sulfate_aer_mod 25 25 USE phys_local_var_mod, ONLY: stratomask … … 325 325 CALL micphy_tstep(pdtphys,tr_seri,t_seri,pplay,paprs,rh,is_strato) 326 326 327 !-- call coagulation routine327 !--CALL coagulation routine 328 328 CALL coagulate(pdtphys,mdw,tr_seri,t_seri,pplay,dens_aer,is_strato) 329 329 330 !-- callsedimentation routine330 !--CALL sedimentation routine 331 331 CALL aer_sedimnt(pdtphys, t_seri, pplay, paprs, tr_seri, dens_aer) 332 332 -
LMDZ6/branches/Amaury_dev/libf/phylmd/acama_gwd_rando_m.F90
r5099 r5101 23 23 use dimphy, only: klon, klev 24 24 use assert_m, only: assert 25 USE ioipsl_getin_p_mod, ONLY 26 USE vertical_layers_mod, ONLY 25 USE ioipsl_getin_p_mod, ONLY: getin_p 26 USE vertical_layers_mod, ONLY: presnivs 27 27 28 28 include "YOMCST.h" … … 120 120 REAL, DIMENSION(klev+1) ::HREF 121 121 LOGICAL, SAVE :: gwd_reproductibilite_mpiomp=.true. 122 LOGICAL, SAVE :: first call= .TRUE.122 LOGICAL, SAVE :: firstCALL = .TRUE. 123 123 !$OMP THREADPRIVATE(firstcall,gwd_reproductibilite_mpiomp) 124 124 … … 199 199 200 200 ! ONLINE 201 callassert(klon == (/size(pp, 1), size(tt, 1), size(uu, 1), &201 CALL assert(klon == (/size(pp, 1), size(tt, 1), size(uu, 1), & 202 202 size(vv, 1), size(rot,1), size(zustr), size(zvstr), size(d_u, 1), & 203 203 size(d_v, 1), & 204 204 size(east_gwstress,1), size(west_gwstress,1) /), & 205 205 "ACAMA_GWD_RANDO klon") 206 callassert(klev == (/size(pp, 2), size(tt, 2), size(uu, 2), &206 CALL assert(klev == (/size(pp, 2), size(tt, 2), size(uu, 2), & 207 207 size(vv, 2), size(d_u, 2), size(d_v, 2), & 208 208 size(east_gwstress,2), size(west_gwstress,2) /), & -
LMDZ6/branches/Amaury_dev/libf/phylmd/add_phys_tend_mod.F90
r5099 r5101 100 100 101 101 USE dimphy, ONLY: klon, klev 102 USE phys_state_var_mod, ONLY 102 USE phys_state_var_mod, ONLY: phys_tstep 103 103 USE phys_local_var_mod, ONLY: u_seri, v_seri, ql_seri, qs_seri, qbs_seri, q_seri, t_seri 104 104 USE phys_state_var_mod, ONLY: ftsol … … 106 106 USE print_control_mod, ONLY: prt_level 107 107 USE cmp_seri_mod 108 USE phys_output_var_mod, ONLY 108 USE phys_output_var_mod, ONLY: d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 109 109 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 110 110 IMPLICIT none … … 192 192 IF (flag_inhib_tend > 0) THEN 193 193 ! print some diagnostics if xxx_seri have changed 194 callcmp_seri(flag_inhib_tend,text)194 CALL cmp_seri(flag_inhib_tend,text) 195 195 ENDIF 196 196 RETURN ! on n ajoute pas les tendance … … 292 292 write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k) 293 293 ENDDO 294 callprint_debug_phys(i,debug_level,text)294 CALL print_debug_phys(i,debug_level,text) 295 295 ENDIF 296 296 ENDDO … … 394 394 write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k) 395 395 ENDDO 396 callprint_debug_phys(i,debug_level,text)396 CALL print_debug_phys(i,debug_level,text) 397 397 ENDIF 398 398 ENDDO … … 412 412 write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k) 413 413 ENDDO 414 callprint_debug_phys(i,debug_level,text)414 CALL print_debug_phys(i,debug_level,text) 415 415 ENDIF 416 416 ENDDO … … 501 501 !====================================================================== 502 502 503 USE phys_state_var_mod, ONLY 503 USE phys_state_var_mod, ONLY: phys_tstep, ftsol 504 504 USE geometry_mod, ONLY: longitude_deg, latitude_deg 505 505 USE print_control_mod, ONLY: prt_level 506 506 USE cmp_seri_mod 507 USE phys_output_var_mod, ONLY 507 USE phys_output_var_mod, ONLY: d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 508 508 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 509 509 IMPLICIT none … … 725 725 726 726 USE dimphy, ONLY: klon, klev 727 USE phys_state_var_mod, ONLY 728 USE phys_state_var_mod, ONLY 727 USE phys_state_var_mod, ONLY: phys_tstep 728 USE phys_state_var_mod, ONLY: topsw, toplw, solsw, sollw, rain_con, snow_con, bs_fall 729 729 USE geometry_mod, ONLY: longitude_deg, latitude_deg 730 730 USE print_control_mod, ONLY: prt_level 731 731 USE cmp_seri_mod 732 USE phys_output_var_mod, ONLY 732 USE phys_output_var_mod, ONLY: d_qw_col, d_ql_col, d_qs_col, d_qbs_col, d_qt_col, d_ek_col, d_h_dair_col & 733 733 , d_h_qw_col, d_h_ql_col, d_h_qs_col, d_h_qbs_col, d_h_col 734 734 USE phys_local_var_mod, ONLY: evap, sens 735 735 USE phys_local_var_mod, ONLY: u_seri, v_seri, ql_seri, qs_seri, qbs_seri, q_seri, t_seri & 736 736 , rain_lsc, snow_lsc 737 USE climb_hq_mod, ONLY 737 USE climb_hq_mod, ONLY: d_h_col_vdf, f_h_bnd 738 738 IMPLICIT none 739 739 INCLUDE "YOMCST.h" -
LMDZ6/branches/Amaury_dev/libf/phylmd/add_wake_tend.F90
r5093 r5101 32 32 IF (prt_level >= 5) then 33 33 write (*,*) "In add_wake_tend, after ",text 34 callflush34 CALL flush 35 35 end if 36 36 -
LMDZ6/branches/Amaury_dev/libf/phylmd/alpale.F90
r5099 r5101 16 16 17 17 USE dimphy 18 USE ioipsl_getin_p_mod, ONLY 18 USE ioipsl_getin_p_mod, ONLY: getin_p 19 19 USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level 20 20 USE phys_local_var_mod, ONLY: zw2 ! Variables internes non sauvegardees de la physique -
LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_th.F90
r5099 r5101 18 18 19 19 USE dimphy 20 USE ioipsl_getin_p_mod, ONLY 20 USE ioipsl_getin_p_mod, ONLY: getin_p 21 21 USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level 22 22 … … 101 101 !! Control of the multiplication of no-trigger probabilities between calls 102 102 !! to the convection scheme. If multiply_proba_notrig is .false., then 103 !! proba_notrig is set to 1 at each call to alpale_th, so that only the last call103 !! proba_notrig is set to 1 at each CALL to alpale_th, so that only the last CALL 104 104 !! plays a role in the triggering of convection. If it is .true., then propa_notrig 105 105 !! is saved between calls to convection and is reset to 1 only after calling the -
LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_wk.F90
r5099 r5101 12 12 13 13 USE dimphy, ONLY: klon 14 USE ioipsl_getin_p_mod, ONLY 14 USE ioipsl_getin_p_mod, ONLY: getin_p 15 15 USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level 16 16 -
LMDZ6/branches/Amaury_dev/libf/phylmd/calbeta_clim.F90
r5099 r5101 9 9 !====================================================================== 10 10 11 !USE phys_local_var_mod, ONLY 11 !USE phys_local_var_mod, ONLY: ideal_beta !pour faire la variable dans le 12 12 ! physiq.f pour des sorties directes de beta 13 13 -
LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_fluxs_mod.F90
r5099 r5101 16 16 17 17 18 USE dimphy, ONLY 18 USE dimphy, ONLY: klon 19 19 USE indice_sol_mod 20 20 use sens_heat_rain_m, only: sens_heat_rain -
LMDZ6/branches/Amaury_dev/libf/phylmd/calltherm.F90
r5099 r5101 238 238 #ifdef ISOVERIF 239 239 if (iso_eau.gt.0) then 240 calliso_verif_egalite_vect2D( &240 CALL iso_verif_egalite_vect2D( & 241 241 xt_seri,q_seri, & 242 242 'calltherm 174',ntiso,klon,klev) … … 408 408 ! write(*,*) 'd_xt_the(iso_hdo,i,k),d_q_the(i,k)=', & 409 409 ! & d_xt_the(iso_hdo,i,k),d_q_the(i,k) 410 calliso_verif_aberrant_enc_vect2D( &410 CALL iso_verif_aberrant_enc_vect2D( & 411 411 xt_seri,q_seri, & 412 412 'calltherm 353, apres ajout d_xt_the',ntiso,klon,klev) … … 450 450 #ifdef ISOVERIF 451 451 if (iso_HDO.gt.0) then 452 calliso_verif_aberrant_enc_vect2D( &452 CALL iso_verif_aberrant_enc_vect2D( & 453 453 xt_seri,q_seri, & 454 454 'calltherm 393, apres bidouille q<0',ntiso,klon,klev) -
LMDZ6/branches/Amaury_dev/libf/phylmd/calwake.F90
r5099 r5101 29 29 USE indice_sol_mod, ONLY: is_oce 30 30 USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level 31 USE lmdz_wake, ONLY 31 USE lmdz_wake, ONLY: wake 32 32 IMPLICIT NONE 33 33 ! ====================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmd/carbon_cycle_mod.F90
r5099 r5101 219 219 SUBROUTINE carbon_cycle_init() 220 220 ! This subroutine is called from tracco2i_init, which is called from phytrac_init only at first timestep. 221 ! - Allocate variables. These variables must be allocated before first callto phys_output_write in physiq.221 ! - Allocate variables. These variables must be allocated before first CALL to phys_output_write in physiq. 222 222 223 223 USE dimphy … … 275 275 276 276 ! USE control_mod, ONLY: planet_type 277 USE phys_cal_mod, ONLY 277 USE phys_cal_mod, ONLY: mth_cur 278 278 USE mod_synchro_omp 279 279 USE mod_phys_lmdz_para, ONLY: is_mpi_root, is_omp_root -
LMDZ6/branches/Amaury_dev/libf/phylmd/cdrag_mod.F90
r5099 r5101 27 27 USE indice_sol_mod 28 28 USE print_control_mod, ONLY: lunout, prt_level 29 USE ioipsl_getin_p_mod, ONLY 30 USE lmdz_atke_turbulence_ini, ONLY 29 USE ioipsl_getin_p_mod, ONLY: getin_p 30 USE lmdz_atke_turbulence_ini, ONLY: smmin, ric, cinf, cepsilon, pr_slope, pr_asym, pr_neut 31 31 32 32 IMPLICIT NONE … … 173 173 174 174 175 LOGICAL, SAVE :: first call= .TRUE.175 LOGICAL, SAVE :: firstCALL = .TRUE. 176 176 !$OMP THREADPRIVATE(firstcall) 177 177 INTEGER, SAVE :: iflag_corr_sta … … 263 263 CALL getin_p('ok_cdrag_iter',ok_cdrag_iter) 264 264 265 first call= .FALSE.265 firstCALL = .FALSE. 266 266 ENDIF 267 267 … … 305 305 mixte = .true. 306 306 ENDIF 307 callclc_core_cp ( sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),t1(i),q1(i),&307 CALL clc_core_cp ( sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),t1(i),q1(i),& 308 308 zgeop1(i)/RG, zgeop1(i)/RG, zgeop1(i)/RG,& 309 309 psol(i),nit_bulk,mixte,& … … 323 323 PSSS = 0.0 324 324 ENDIF 325 callini_csts326 callecumev6_flux( z_0m,t1(i),tsurf(i),&325 CALL ini_csts 326 CALL ecumev6_flux( z_0m,t1(i),tsurf(i),& 327 327 q1(i),qsurf(i),sqrt(zdu2),zgeop1(i)/RG,PSSS,zgeop1(i)/RG,& 328 328 psol(i),pat1(i), OPRECIP, OPWEBB,& … … 342 342 LPRECIP = .false. 343 343 LPWG = .false. 344 callini_csts344 CALL ini_csts 345 345 block 346 346 real, dimension(1) :: z0m_1d, z_0h_1d, sqrt_zdu2_1d, zgeop1_rg_1d ! convert scalar to 1D for call … … 349 349 sqrt_zdu2_1d = sqrt(zdu2) 350 350 zgeop1_rg_1d=zgeop1(i)/RG 351 callcoare30_flux_cnrm(z0m_1d,t1(i),tsurf(i), q1(i), &351 CALL coare30_flux_cnrm(z0m_1d,t1(i),tsurf(i), q1(i), & 352 352 sqrt_zdu2_1d,zgeop1_rg_1d,zgeop1_rg_1d,psol(i),qsurf(i),PQSAT, & 353 353 PSFTH,PFSTQ,PUSTAR,PCD,PCDN,PCH,PCE,PRI, & … … 367 367 ENDIF 368 368 write(*,*) "debug size",size(coeffs) 369 callcoare_cp(sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),&369 CALL coare_cp(sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),& 370 370 t1(i),q1(i),& 371 371 zgeop1(i)/RG,zgeop1(i)/RG,zgeop1(i)/RG,& -
LMDZ6/branches/Amaury_dev/libf/phylmd/change_srf_frac_mod.F90
r5099 r5101 23 23 24 24 USE dimphy 25 USE surface_data, ONLY 25 USE surface_data, ONLY: type_ocean,version_ocean 26 26 USE limit_read_mod 27 USE pbl_surface_mod, ONLY 28 USE cpl_mod, ONLY 29 USE ocean_slab_mod, ONLY 27 USE pbl_surface_mod, ONLY: pbl_surface_newfrac 28 USE cpl_mod, ONLY: cpl_receive_frac 29 USE ocean_slab_mod, ONLY: fsic, ocean_slab_frac 30 30 USE indice_sol_mod 31 31 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/cltracrn.F90
r5099 r5101 8 8 9 9 USE dimphy 10 USE traclmdz_mod, ONLY 10 USE traclmdz_mod, ONLY: id_rn, id_pb 11 11 USE indice_sol_mod 12 12 -
LMDZ6/branches/Amaury_dev/libf/phylmd/coare30_flux_cnrm_mod.F90
r5099 r5101 70 70 71 71 72 !USE MODD_SEAFLUX_n, ONLY 72 !USE MODD_SEAFLUX_n, ONLY: SEAFLUX_t 73 73 74 74 !----------Rajout Olive --------- … … 79 79 !-------------------------------- 80 80 81 USE MODD_CSTS, ONLY 81 USE MODD_CSTS, ONLY: XKARMAN, XG, XSTEFAN, XRD, XRV, XPI, & 82 82 XLVTT, XCL, XCPD, XCPV, XRHOLW, XTT, & 83 83 XP00 84 84 85 !USE MODD_SURF_ATM, ONLY 86 87 !USE MODD_SURF_PAR, ONLY 85 !USE MODD_SURF_ATM, ONLY: XVZ0CM 86 87 !USE MODD_SURF_PAR, ONLY: XUNDEF, XSURF_EPSILON 88 88 !USE MODD_WATER_PAR 89 89 … … 98 98 99 99 100 !USE YOMHOOK ,ONLY 101 !USE PARKIND1 ,ONLY 100 !USE YOMHOOK ,ONLY: LHOOK, DR_HOOK 101 !USE PARKIND1 ,ONLY: JPRB 102 102 103 103 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/coare_cp_mod.F90
r5099 r5101 61 61 !taylor and yelland zo=1200 h*(L/h)**4.5 jwave=2 62 62 63 USE MODD_CSTS, ONLY 63 USE MODD_CSTS, ONLY: XKARMAN, XG, XSTEFAN, XRD, XRV, XPI, & 64 64 XLVTT, XCL, XCPD, XCPV, XRHOLW, XTT, & 65 65 XP00 -
LMDZ6/branches/Amaury_dev/libf/phylmd/concvl.F90
r5099 r5101 529 529 ! c e , t1,q1,ql,qs,u,v,paprs,pplay 530 530 ! c s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec) 531 ! c calldiagphy(paire,ztit,ip_ebil531 ! c CALL diagphy(paire,ztit,ip_ebil 532 532 ! c e , zero_v, zero_v, zero_v, zero_v, zero_v 533 533 ! c e , zero_v, rain, zero_v, ztsol -
LMDZ6/branches/Amaury_dev/libf/phylmd/conemav.F90
r2346 r5101 119 119 rflag(i) = iflag(i) 120 120 END DO 121 ! calldump2d(iim,jjm-1,rflag(2:klon-1),'FLAG CONVECTION ')121 ! CALL dump2d(iim,jjm-1,rflag(2:klon-1),'FLAG CONVECTION ') 122 122 ! if (klon.eq.1) then 123 123 ! print*,'IFLAG ',iflag -
LMDZ6/branches/Amaury_dev/libf/phylmd/conf_phys_m.F90
r5099 r5101 26 26 USE surface_data 27 27 USE phys_cal_mod 28 USE carbon_cycle_mod, ONLY 29 USE carbon_cycle_mod, ONLY 30 USE carbon_cycle_mod, ONLY 31 USE chemistry_cycle_mod, ONLY 32 USE mod_grid_phy_lmdz, ONLY 33 USE print_control_mod, ONLY 34 use config_ocean_skin_m, only: config_ocean_skin35 USE phys_state_var_mod, ONLY 36 USE infotrac_phy, ONLY 28 USE carbon_cycle_mod, ONLY: carbon_cycle_tr, carbon_cycle_cpl, carbon_cycle_rad, level_coupling_esm 29 USE carbon_cycle_mod, ONLY: read_fco2_ocean_cor, var_fco2_ocean_cor 30 USE carbon_cycle_mod, ONLY: read_fco2_land_cor, var_fco2_land_cor 31 USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl 32 USE mod_grid_phy_lmdz, ONLY: klon_glo 33 USE print_control_mod, ONLY: lunout 34 use config_ocean_skin_m, ONLY: config_ocean_skin 35 USE phys_state_var_mod, ONLY: phys_tstep 36 USE infotrac_phy, ONLY: type_trac 37 37 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER 38 38 … … 1464 1464 1465 1465 iflag_pbl_split_omp = 0 1466 callgetin('iflag_pbl_split', iflag_pbl_split_omp)1466 CALL getin('iflag_pbl_split', iflag_pbl_split_omp) 1467 1467 1468 1468 !Config Key = iflag_order2_sollw … … 2835 2835 2836 2836 !$OMP END MASTER 2837 callconfig_ocean_skin2837 CALL config_ocean_skin 2838 2838 2839 2839 END SUBROUTINE conf_phys … … 2846 2846 2847 2847 USE IOIPSL 2848 USE print_control_mod, ONLY 2848 USE print_control_mod, ONLY: lunout 2849 2849 IMPLICIT NONE 2850 2850 ! Configuration de l'interace atm/surf -
LMDZ6/branches/Amaury_dev/libf/phylmd/conflx.F90
r2346 r5101 539 539 END DO 540 540 541 i call= 0541 iCALL = 0 542 542 CALL flxadjtq(paph(1,k), ptenh(1,k), pqsenh(1,k), llflag, icall) 543 543 … … 660 660 END DO 661 661 662 i call= 1662 iCALL = 1 663 663 CALL flxadjtq(paph(1,k), ptu(1,k), pqu(1,k), llflag, icall) 664 664 … … 903 903 ! ---------------------------------------------------------------------- 904 904 905 i call= 1905 iCALL = 1 906 906 CALL flxadjtq(paph(1,k), ptu(1,k), pqu(1,k), llflag, icall) 907 907 … … 1359 1359 IF (is==0) GO TO 290 1360 1360 1361 i call= 21361 iCALL = 2 1362 1362 CALL flxadjtq(paph(1,k), ztenwb(1,k), zqenwb(1,k), llo2, icall) 1363 1363 … … 1485 1485 END DO 1486 1486 1487 i call= 21487 iCALL = 2 1488 1488 CALL flxadjtq(paph(1,k), ptd(1,k), pqd(1,k), llo2, icall) 1489 1489 … … 1512 1512 ! Objet: ajustement entre T et Q 1513 1513 ! ====================================================================== 1514 ! NOTE: INPUT PARAMETER k callDEFINES CALCULATION AS1514 ! NOTE: INPUT PARAMETER kCALL DEFINES CALCULATION AS 1515 1515 ! kcall=0 ENV. T AND QS IN*CUINI* 1516 1516 ! kcall=1 CONDENSATION IN UPDRAFTS (E.G. CUBASE, CUASC) -
LMDZ6/branches/Amaury_dev/libf/phylmd/convect1.F90
r2346 r5101 33 33 ! LIBRARIES OF RESIDENCE: /a/ops/lib/libfcst159.a 34 34 35 ! USAGE: callconvect1(len,nd,noff,minorig,35 ! USAGE: CALL convect1(len,nd,noff,minorig, 36 36 ! & t,q,qs,u,v, 37 37 ! & p,ph,iflag,ft, -
LMDZ6/branches/Amaury_dev/libf/phylmd/convect2.F90
r2346 r5101 36 36 ! LIBRARIES OF RESIDENCE: /a/ops/lib/libfcst159.a 37 37 38 ! USAGE: callconvect2(ncum,idcum,len,nd,nl,minorig,38 ! USAGE: CALL convect2(ncum,idcum,len,nd,nl,minorig, 39 39 ! & nk1,icb1, 40 40 ! & t1,q1,qs1,u1,v1,gz1,tv1,tp1,tvp1,clw1,h1, … … 656 656 ! K Emanuel fix 657 657 658 ! callzilch(byp,ncum)658 ! CALL zilch(byp,ncum) 659 659 ! do 530 k=minorig+1,nl-1 660 660 ! do 520 i=1,ncum -
LMDZ6/branches/Amaury_dev/libf/phylmd/convect3.F90
r2346 r5101 1343 1343 1344 1344 ! @$$cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 1345 ! @$$ callwriteg1d(1,klev,ma,'ma ','ma ')1346 ! @$$ callwriteg1d(1,klev,upwd,'upwd ','upwd ')1347 ! @$$ callwriteg1d(1,klev,dnwd,'dnwd ','dnwd ')1348 ! @$$ callwriteg1d(1,klev,dnwd0,'dnwd0 ','dnwd0 ')1349 ! @$$ callwriteg1d(1,klev,tvp,'tvp ','tvp ')1350 ! @$$ callwriteg1d(1,klev,tra(1:klev,3),'tra3 ','tra3 ')1351 ! @$$ callwriteg1d(1,klev,tra(1:klev,4),'tra4 ','tra4 ')1352 ! @$$ callwriteg1d(1,klev,tra(1:klev,5),'tra5 ','tra5 ')1353 ! @$$ callwriteg1d(1,klev,tra(1:klev,6),'tra6 ','tra6 ')1354 ! @$$ callwriteg1d(1,klev,tra(1:klev,7),'tra7 ','tra7 ')1355 ! @$$ callwriteg1d(1,klev,tra(1:klev,8),'tra8 ','tra8 ')1356 ! @$$ callwriteg1d(1,klev,tra(1:klev,9),'tra9 ','tra9 ')1357 ! @$$ callwriteg1d(1,klev,tra(1:klev,10),'tra10','tra10')1358 ! @$$ callwriteg1d(1,klev,tra(1:klev,11),'tra11','tra11')1359 ! @$$ callwriteg1d(1,klev,tra(1:klev,12),'tra12','tra12')1360 ! @$$ callwriteg1d(1,klev,tra(1:klev,13),'tra13','tra13')1361 ! @$$ callwriteg1d(1,klev,tra(1:klev,14),'tra14','tra14')1362 ! @$$ callwriteg1d(1,klev,tra(1:klev,15),'tra15','tra15')1363 ! @$$ callwriteg1d(1,klev,tra(1:klev,16),'tra16','tra16')1364 ! @$$ callwriteg1d(1,klev,tra(1:klev,17),'tra17','tra17')1365 ! @$$ callwriteg1d(1,klev,tra(1:klev,18),'tra18','tra18')1366 ! @$$ callwriteg1d(1,klev,tra(1:klev,19),'tra19','tra19')1367 ! @$$ callwriteg1d(1,klev,tra(1:klev,20),'tra20','tra20 ')1368 ! @$$ callwriteg1d(1,klev,trap(1:klev,1),'trp1','trp1')1369 ! @$$ callwriteg1d(1,klev,trap(1:klev,2),'trp2','trp2')1370 ! @$$ callwriteg1d(1,klev,trap(1:klev,3),'trp3','trp3')1371 ! @$$ callwriteg1d(1,klev,trap(1:klev,4),'trp4','trp4')1372 ! @$$ callwriteg1d(1,klev,trap(1:klev,5),'trp5','trp5')1373 ! @$$ callwriteg1d(1,klev,trap(1:klev,10),'trp10','trp10')1374 ! @$$ callwriteg1d(1,klev,trap(1:klev,12),'trp12','trp12')1375 ! @$$ callwriteg1d(1,klev,trap(1:klev,15),'trp15','trp15')1376 ! @$$ callwriteg1d(1,klev,trap(1:klev,20),'trp20','trp20')1377 ! @$$ callwriteg1d(1,klev,ftra(1:klev,1),'ftr1 ','ftr1 ')1378 ! @$$ callwriteg1d(1,klev,ftra(1:klev,2),'ftr2 ','ftr2 ')1379 ! @$$ callwriteg1d(1,klev,ftra(1:klev,3),'ftr3 ','ftr3 ')1380 ! @$$ callwriteg1d(1,klev,ftra(1:klev,4),'ftr4 ','ftr4 ')1381 ! @$$ callwriteg1d(1,klev,ftra(1:klev,5),'ftr5 ','ftr5 ')1382 ! @$$ callwriteg1d(1,klev,ftra(1:klev,6),'ftr6 ','ftr6 ')1383 ! @$$ callwriteg1d(1,klev,ftra(1:klev,7),'ftr7 ','ftr7 ')1384 ! @$$ callwriteg1d(1,klev,ftra(1:klev,8),'ftr8 ','ftr8 ')1385 ! @$$ callwriteg1d(1,klev,ftra(1:klev,9),'ftr9 ','ftr9 ')1386 ! @$$ callwriteg1d(1,klev,ftra(1:klev,10),'ftr10','ftr10')1387 ! @$$ callwriteg1d(1,klev,ftra(1:klev,11),'ftr11','ftr11')1388 ! @$$ callwriteg1d(1,klev,ftra(1:klev,12),'ftr12','ftr12')1389 ! @$$ callwriteg1d(1,klev,ftra(1:klev,13),'ftr13','ftr13')1390 ! @$$ callwriteg1d(1,klev,ftra(1:klev,14),'ftr14','ftr14')1391 ! @$$ callwriteg1d(1,klev,ftra(1:klev,15),'ftr15','ftr15')1392 ! @$$ callwriteg1d(1,klev,ftra(1:klev,16),'ftr16','ftr16')1393 ! @$$ callwriteg1d(1,klev,ftra(1:klev,17),'ftr17','ftr17')1394 ! @$$ callwriteg1d(1,klev,ftra(1:klev,18),'ftr18','ftr18')1395 ! @$$ callwriteg1d(1,klev,ftra(1:klev,19),'ftr19','ftr19')1396 ! @$$ callwriteg1d(1,klev,ftra(1:klev,20),'ftr20','ftr20 ')1397 ! @$$ callwriteg1d(1,klev,mp,'mp ','mp ')1398 ! @$$ callwriteg1d(1,klev,Mke,'Mke ','Mke ')1345 ! @$$ CALL writeg1d(1,klev,ma,'ma ','ma ') 1346 ! @$$ CALL writeg1d(1,klev,upwd,'upwd ','upwd ') 1347 ! @$$ CALL writeg1d(1,klev,dnwd,'dnwd ','dnwd ') 1348 ! @$$ CALL writeg1d(1,klev,dnwd0,'dnwd0 ','dnwd0 ') 1349 ! @$$ CALL writeg1d(1,klev,tvp,'tvp ','tvp ') 1350 ! @$$ CALL writeg1d(1,klev,tra(1:klev,3),'tra3 ','tra3 ') 1351 ! @$$ CALL writeg1d(1,klev,tra(1:klev,4),'tra4 ','tra4 ') 1352 ! @$$ CALL writeg1d(1,klev,tra(1:klev,5),'tra5 ','tra5 ') 1353 ! @$$ CALL writeg1d(1,klev,tra(1:klev,6),'tra6 ','tra6 ') 1354 ! @$$ CALL writeg1d(1,klev,tra(1:klev,7),'tra7 ','tra7 ') 1355 ! @$$ CALL writeg1d(1,klev,tra(1:klev,8),'tra8 ','tra8 ') 1356 ! @$$ CALL writeg1d(1,klev,tra(1:klev,9),'tra9 ','tra9 ') 1357 ! @$$ CALL writeg1d(1,klev,tra(1:klev,10),'tra10','tra10') 1358 ! @$$ CALL writeg1d(1,klev,tra(1:klev,11),'tra11','tra11') 1359 ! @$$ CALL writeg1d(1,klev,tra(1:klev,12),'tra12','tra12') 1360 ! @$$ CALL writeg1d(1,klev,tra(1:klev,13),'tra13','tra13') 1361 ! @$$ CALL writeg1d(1,klev,tra(1:klev,14),'tra14','tra14') 1362 ! @$$ CALL writeg1d(1,klev,tra(1:klev,15),'tra15','tra15') 1363 ! @$$ CALL writeg1d(1,klev,tra(1:klev,16),'tra16','tra16') 1364 ! @$$ CALL writeg1d(1,klev,tra(1:klev,17),'tra17','tra17') 1365 ! @$$ CALL writeg1d(1,klev,tra(1:klev,18),'tra18','tra18') 1366 ! @$$ CALL writeg1d(1,klev,tra(1:klev,19),'tra19','tra19') 1367 ! @$$ CALL writeg1d(1,klev,tra(1:klev,20),'tra20','tra20 ') 1368 ! @$$ CALL writeg1d(1,klev,trap(1:klev,1),'trp1','trp1') 1369 ! @$$ CALL writeg1d(1,klev,trap(1:klev,2),'trp2','trp2') 1370 ! @$$ CALL writeg1d(1,klev,trap(1:klev,3),'trp3','trp3') 1371 ! @$$ CALL writeg1d(1,klev,trap(1:klev,4),'trp4','trp4') 1372 ! @$$ CALL writeg1d(1,klev,trap(1:klev,5),'trp5','trp5') 1373 ! @$$ CALL writeg1d(1,klev,trap(1:klev,10),'trp10','trp10') 1374 ! @$$ CALL writeg1d(1,klev,trap(1:klev,12),'trp12','trp12') 1375 ! @$$ CALL writeg1d(1,klev,trap(1:klev,15),'trp15','trp15') 1376 ! @$$ CALL writeg1d(1,klev,trap(1:klev,20),'trp20','trp20') 1377 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,1),'ftr1 ','ftr1 ') 1378 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,2),'ftr2 ','ftr2 ') 1379 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,3),'ftr3 ','ftr3 ') 1380 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,4),'ftr4 ','ftr4 ') 1381 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,5),'ftr5 ','ftr5 ') 1382 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,6),'ftr6 ','ftr6 ') 1383 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,7),'ftr7 ','ftr7 ') 1384 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,8),'ftr8 ','ftr8 ') 1385 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,9),'ftr9 ','ftr9 ') 1386 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,10),'ftr10','ftr10') 1387 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,11),'ftr11','ftr11') 1388 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,12),'ftr12','ftr12') 1389 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,13),'ftr13','ftr13') 1390 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,14),'ftr14','ftr14') 1391 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,15),'ftr15','ftr15') 1392 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,16),'ftr16','ftr16') 1393 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,17),'ftr17','ftr17') 1394 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,18),'ftr18','ftr18') 1395 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,19),'ftr19','ftr19') 1396 ! @$$ CALL writeg1d(1,klev,ftra(1:klev,20),'ftr20','ftr20 ') 1397 ! @$$ CALL writeg1d(1,klev,mp,'mp ','mp ') 1398 ! @$$ CALL writeg1d(1,klev,Mke,'Mke ','Mke ') 1399 1399 1400 1400 -
LMDZ6/branches/Amaury_dev/libf/phylmd/cpl_mod.F90
r5099 r5101 15 15 ! Use statements 16 16 !************************************************************************************* 17 USE dimphy, ONLY 17 USE dimphy, ONLY: klon 18 18 USE mod_phys_lmdz_para 19 19 USE ioipsl … … 133 133 134 134 SUBROUTINE cpl_init(dtime, rlon, rlat) 135 USE carbon_cycle_mod, ONLY 135 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, fco2_ocn_day 136 136 USE surface_data 137 137 USE indice_sol_mod 138 USE mod_grid_phy_lmdz, ONLY 138 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid1dTo2d_glo, klon_glo, grid_type, unstructured, regular_lonlat 139 139 USE time_phylmdz_mod, ONLY: annee_ref, day_ini, itau_phy, itaufin_phy 140 140 USE print_control_mod, ONLY: lunout 141 USE geometry_mod, ONLY 141 USE geometry_mod, ONLY: longitude_deg, latitude_deg, ind_cell_glo, cell_area 142 142 USE ioipsl_getin_p_mod, ONLY: getin_p 143 143 use config_ocean_skin_m, only: activate_ocean_skin … … 407 407 ! are stored in this module. 408 408 USE surface_data 409 USE geometry_mod, ONLY 410 USE carbon_cycle_mod, ONLY 409 USE geometry_mod, ONLY: longitude_deg, latitude_deg 410 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 411 411 USE indice_sol_mod 412 412 USE time_phylmdz_mod, ONLY: start_time, itau_phy 413 USE mod_grid_phy_lmdz, ONLY 413 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 414 414 use config_ocean_skin_m, only: activate_ocean_skin 415 415 … … 537 537 ! The temperature is transformed into 1D array with valid points from index 1 to knon. 538 538 539 USE carbon_cycle_mod, ONLY 539 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, fco2_ocn_day 540 540 USE indice_sol_mod 541 541 use config_ocean_skin_m, only: activate_ocean_skin … … 659 659 ! fractions do not change between coupling time-steps. 660 660 661 USE carbon_cycle_mod, ONLY 661 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, co2_send 662 662 USE indice_sol_mod 663 USE mod_grid_phy_lmdz, ONLY 663 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 664 664 use config_ocean_skin_m, only: activate_ocean_skin 665 665 … … 919 919 ! the coupler. 920 920 921 USE carbon_cycle_mod, ONLY 921 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 922 922 USE indice_sol_mod 923 USE mod_grid_phy_lmdz, ONLY 923 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 924 924 925 925 ! Input arguments … … 1117 1117 ! (it is done in cpl_send_seaice_fields). 1118 1118 1119 USE mod_grid_phy_lmdz, ONLY 1119 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 1120 1120 1121 1121 ! Input arguments … … 1173 1173 ! will be done in cpl_send_seaice_fields. 1174 1174 1175 USE mod_grid_phy_lmdz, ONLY 1175 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 1176 1176 1177 1177 ! Input varibales … … 1227 1227 1228 1228 USE surface_data 1229 USE carbon_cycle_mod, ONLY 1229 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 1230 1230 USE indice_sol_mod 1231 USE mod_grid_phy_lmdz, ONLY 1231 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 1232 1232 USE time_phylmdz_mod, ONLY: start_time, itau_phy 1233 1233 USE config_ocean_skin_m, only: activate_ocean_skin … … 1554 1554 ! champ_out champ sur la grille 'gatherd' 1555 1555 1556 USE mod_grid_phy_lmdz, ONLY 1556 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 1557 1557 1558 1558 ! Input … … 1601 1601 ! champ_out champ sur la grille 2D 1602 1602 1603 USE mod_grid_phy_lmdz, ONLY 1603 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 1604 1604 1605 1605 ! Input arguments -
LMDZ6/branches/Amaury_dev/libf/phylmd/create_etat0_limit_unstruct_mod.F90
r4938 r5101 14 14 USE mod_phys_lmdz_para, ONLY: is_omp_master 15 15 USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured 16 USE ioipsl, ONLY 16 USE ioipsl, ONLY: ioget_year_len 17 17 USE ioipsl_getin_p_mod, ONLY: getin_p 18 USE time_phylmdz_mod, ONLY 18 USE time_phylmdz_mod, ONLY: annee_ref 19 19 USE create_etat0_unstruct_mod, ONLY: init_create_etat0_unstruct 20 20 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/create_limit_unstruct_mod.F90
r5099 r5101 11 11 USE dimphy 12 12 USE lmdz_xios 13 USE ioipsl, ONLY 14 USE time_phylmdz_mod, ONLY 13 USE ioipsl, ONLY: ioget_year_len 14 USE time_phylmdz_mod, ONLY: annee_ref 15 15 USE indice_sol_mod 16 16 USE phys_state_var_mod … … 181 181 USE pchfe_95_m, only: pchfe_95 182 182 USE arth_m, only: arth 183 USE dimphy, ONLY 184 USE ioipsl, ONLY 185 USE time_phylmdz_mod, ONLY 183 USE dimphy, ONLY: klon 184 USE ioipsl, ONLY: ioget_year_len 185 USE time_phylmdz_mod, ONLY: annee_ref 186 186 USE mod_phys_lmdz_para 187 187 IMPLICIT NONE … … 231 231 232 232 !------------------------------------------------------------------------------- 233 USE ioipsl, ONLY 233 USE ioipsl, ONLY: ioget_calendar,ioconf_calendar,lock_calendar,ioget_year_len 234 234 IMPLICIT NONE 235 235 !------------------------------------------------------------------------------- … … 264 264 265 265 !------------------------------------------------------------------------------- 266 USE ioipsl, ONLY 266 USE ioipsl, ONLY: ioget_calendar,ioconf_calendar,lock_calendar,ioget_mon_len 267 267 IMPLICIT NONE 268 268 !------------------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/ctstar.F90
r2386 r5101 61 61 62 62 !USE PARKIND1 63 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/parkind1.F90.php#parkind1> ,ONLY 63 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/parkind1.F90.php#parkind1> ,ONLY: JPIM ,JPRB 64 64 !USE YOMHOOK 65 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/yomhook.F90.php#yomhook> ,ONLY 65 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/yomhook.F90.php#yomhook> ,ONLY: LHOOK, DR_HOOK 66 66 67 !USE YOMCST, ONLY : RG, RD68 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomcst.F90.php#yomcst> , ONLY 67 !USE YOMCST, ONLY: RG, RD 68 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomcst.F90.php#yomcst> , ONLY: RG 69 69 70 70 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/support/rg.F.php#rg> ,RD 71 71 !USE YOMSTA 72 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomsta.F90.php#yomsta> , ONLY 72 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomsta.F90.php#yomsta> , ONLY: RDTDZ1 73 73 74 74 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90
r5093 r5101 1129 1129 ! K Emanuel fix 1130 1130 1131 ! callzilch(byp,ncum)1131 ! CALL zilch(byp,ncum) 1132 1132 ! do 530 k=minorig+1,nl-1 1133 1133 ! do 520 i=1,ncum … … 1156 1156 ! J Teixeira fix 1157 1157 1158 ! ori callzilch(byp,ncum)1158 ! ori CALL zilch(byp,ncum) 1159 1159 ! ori do 515 i=1,ncum 1160 1160 ! ori lcape(i)=.true. … … 1624 1624 ! ===================================================================== 1625 1625 1626 ! ym callzilch(asum,ncum*nd)1627 ! ym callzilch(bsum,ncum*nd)1628 ! ym callzilch(csum,ncum*nd)1626 ! ym CALL zilch(asum,ncum*nd) 1627 ! ym CALL zilch(bsum,ncum*nd) 1628 ! ym CALL zilch(csum,ncum*nd) 1629 1629 CALL zilch(asum, nloc*nd) 1630 1630 CALL zilch(csum, nloc*nd) -
LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90
r5099 r5101 7 7 SUBROUTINE cv3_param(nd, k_upper, delt) 8 8 9 USE ioipsl_getin_p_mod, ONLY 9 USE ioipsl_getin_p_mod, ONLY: getin_p 10 10 use mod_phys_lmdz_para 11 11 IMPLICIT NONE … … 309 309 cpnk, hnk, nk, icb, icbmax, iflag, gznk, plcl) 310 310 311 USE mod_phys_lmdz_transfert_para, ONLY 311 USE mod_phys_lmdz_transfert_para, ONLY: bcast 312 312 USE add_phys_tend_mod, ONLY: fl_cor_ebil 313 313 USE print_control_mod, ONLY: prt_level … … 382 382 PRINT *, ' ok_new_feed: ', ok_new_feed 383 383 !$OMP END MASTER 384 callbcast(ok_new_feed)384 CALL bcast(ok_new_feed) 385 385 first = .FALSE. 386 386 END IF … … 579 579 ! print*,'icb dans cv3_feed ' 580 580 ! write(*,'(64i2)') icb(2:len-1) 581 ! calldump2d(64,43,'plcl dans cv3_feed ',plcl(2:len-1))581 ! CALL dump2d(64,43,'plcl dans cv3_feed ',plcl(2:len-1)) 582 582 583 583 DO i = 1, len … … 1912 1912 ! K Emanuel fix 1913 1913 1914 ! callzilch(byp,ncum)1914 ! CALL zilch(byp,ncum) 1915 1915 ! do 530 k=minorig+1,nl-1 1916 1916 ! do 520 i=1,ncum … … 1939 1939 ! J Teixeira fix 1940 1940 1941 ! ori callzilch(byp,ncum)1941 ! ori CALL zilch(byp,ncum) 1942 1942 ! ori do 515 i=1,ncum 1943 1943 ! ori lcape(i)=.true. … … 3450 3450 3451 3451 USE print_control_mod, ONLY: lunout, prt_level 3452 USE add_phys_tend_mod, only: fl_cor_ebil3452 USE add_phys_tend_mod, ONLY: fl_cor_ebil 3453 3453 3454 3454 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/cv_driver.F90
r5099 r5101 565 565 566 566 ! epmax_cape 567 callcv30_epmax_fn_cape(nloc,ncum,nd &567 CALL cv30_epmax_fn_cape(nloc,ncum,nd & 568 568 ,cape,ep,hp,icb,inb,clw,nk,t,h,lv & 569 569 ,epmax_diag) … … 682 682 SUBROUTINE cv_flag(iflag_ice_thermo) 683 683 684 USE ioipsl_getin_p_mod, ONLY 684 USE ioipsl_getin_p_mod, ONLY: getin_p 685 685 686 686 IMPLICIT NONE … … 708 708 ! calculee en une seule iteration. 709 709 icvflag_Tpa=0 710 callgetin_p('icvflag_Tpa', icvflag_Tpa)710 CALL getin_p('icvflag_Tpa', icvflag_Tpa) 711 711 712 712 RETURN -
LMDZ6/branches/Amaury_dev/libf/phylmd/cv_routines.F90
r2311 r5101 686 686 ! K Emanuel fix 687 687 688 ! callzilch(byp,ncum)688 ! CALL zilch(byp,ncum) 689 689 ! do 530 k=minorig+1,nl-1 690 690 ! do 520 i=1,ncum -
LMDZ6/branches/Amaury_dev/libf/phylmd/cva_driver.F90
r5099 r5101 814 814 ! write(*,'(64i1)') iflag1(2:len-1) 815 815 816 ! calldump2d(iim,jjm-1,sig1(2)816 ! CALL dump2d(iim,jjm-1,sig1(2) 817 817 END IF 818 818 … … 942 942 if (prt_level >= 9) & 943 943 PRINT *, 'cva_driver -> cv3_epmax_cape' 944 callcv3_epmax_fn_cape(nloc,ncum,nd &944 CALL cv3_epmax_fn_cape(nloc,ncum,nd & 945 945 , ep,hp,icb,inb,clw,nk,t,h,hnk,lv,lf,frac_s & 946 946 , pbase, p, ph, tv, buoy, sig, w0,iflag & -
LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr.F90
r5099 r5101 12 12 USE IOIPSL 13 13 USE dimphy 14 USE infotrac_phy, ONLY 14 USE infotrac_phy, ONLY: nbtr 15 15 IMPLICIT NONE 16 16 !===================================================================== … … 125 125 126 126 !$OMP MASTER 127 callgetin('ccntrAA_coef',ccntrAA_coef)128 callgetin('ccntrENV_coef',ccntrENV_coef)129 callgetin('coefcoli',coefcoli)127 CALL getin('ccntrAA_coef',ccntrAA_coef) 128 CALL getin('ccntrENV_coef',ccntrENV_coef) 129 CALL getin('coefcoli',coefcoli) 130 130 !$OMP END MASTER 131 131 !$OMP BARRIER -
LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr_noscav.F90
r5099 r5101 4 4 SUBROUTINE cvltr_noscav(it,pdtime,da, phi, mp,wght_cvfd,paprs,pplay,x,upd,dnd,dx) 5 5 USE dimphy 6 USE infotrac_phy, ONLY 6 USE infotrac_phy, ONLY: nbtr 7 7 IMPLICIT NONE 8 8 !===================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr_scav.F90
r5099 r5101 13 13 USE IOIPSL 14 14 USE dimphy 15 USE infotrac_phy, ONLY 15 USE infotrac_phy, ONLY: nbtr 16 16 IMPLICIT NONE 17 17 !===================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr_spl.F90
r5099 r5101 13 13 USE IOIPSL 14 14 USE dimphy 15 USE infotrac_phy, ONLY 15 USE infotrac_phy, ONLY: nbtr 16 16 IMPLICIT NONE 17 17 !===================================================================== … … 152 152 153 153 !!$OMP MASTER 154 ! callgetin('ccntrAA_coef',ccntrAA_coef)155 ! callgetin('ccntrENV_coef',ccntrENV_coef)156 ! callgetin('coefcoli',coefcoli)154 ! CALL getin('ccntrAA_coef',ccntrAA_coef) 155 ! CALL getin('ccntrENV_coef',ccntrENV_coef) 156 ! CALL getin('coefcoli',coefcoli) 157 157 !!$OMP END MASTER 158 158 !!$OMP BARRIER -
LMDZ6/branches/Amaury_dev/libf/phylmd/cvltrorig.F90
r5099 r5101 4 4 SUBROUTINE cvltrorig(it,pdtime,da, phi, mp,paprs,pplay,x,upd,dnd,dx) 5 5 USE dimphy 6 USE infotrac_phy, ONLY 6 USE infotrac_phy, ONLY: nbtr 7 7 IMPLICIT NONE 8 8 !===================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/1DUTILS.h
r5099 r5101 1007 1007 1008 1008 write(*,*) 'in abort_gcm' 1009 callhistclo1010 ! callhistclo(2)1011 ! callhistclo(3)1012 ! callhistclo(4)1013 ! callhistclo(5)1009 CALL histclo 1010 ! CALL histclo(2) 1011 ! CALL histclo(3) 1012 ! CALL histclo(4) 1013 ! CALL histclo(5) 1014 1014 write(*,*) 'out of histclo' 1015 1015 write(*,*) 'Stopping in ', modname 1016 1016 write(*,*) 'Reason = ',message 1017 callgetin_dump1017 CALL getin_dump 1018 1018 1019 1019 if (ierr .eq. 0) then … … 1271 1271 !! This subroutine specifies the surface temperature to be used in 1D simulations 1272 1272 1273 ! USE dimphy, ONLY 1273 ! USE dimphy, ONLY: klon 1274 1274 1275 1275 ! INTEGER, INTENT(IN) :: knon ! nomber of points on compressed grid … … 1485 1485 ! endif 1486 1486 1487 ! callwritefield_phy('d_t_va',d_t_va,llm)1487 ! CALL writefield_phy('d_t_va',d_t_va,llm) 1488 1488 1489 1489 return -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/1D_read_forc_cases.h
r5099 r5101 14 14 if (forcing_SCM) then 15 15 16 write(*,*),'avant callread_SCM'17 callread_SCM_cas16 write(*,*),'avant CALL read_SCM' 17 CALL read_SCM_cas 18 18 write(*,*) 'Forcing read' 19 19 print*,'PS ps_cas',ps_cas -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/1Dconv.h
r5099 r5101 187 187 open(99,FILE=file_fordat,FORM='UNFORMATTED', & 188 188 & ACCESS='DIRECT',RECL=8) 189 callrdgrads(99,icompt,nblvlm,z,ht_mes,hq_mes,hw_mes &189 CALL rdgrads(99,icompt,nblvlm,z,ht_mes,hq_mes,hw_mes & 190 190 & ,hu_mes,hv_mes,hthturb_mes,hqturb_mes & 191 191 & ,ts_fcg,ts_subr,imp_fcg,Turb_fcg) … … 333 333 aaa=spaces(aaa,1) 334 334 print*,'aaa',aaa 335 callgetsch(aaa,' ',' ',5,atemps,nch)335 CALL getsch(aaa,' ',' ',5,atemps,nch) 336 336 print*,'atemps est',atemps 337 337 atemps=atemps(1:nch-2) … … 340 340 dt=imn*60 341 341 print*,'le pas de temps dt',dt 342 callgetsch(aaa,' ',' ',2,apasmax,nch)342 CALL getsch(aaa,' ',' ',2,apasmax,nch) 343 343 apasmax=apasmax(1:nch) 344 344 read(apasmax,*) ipa … … 365 365 IF (Turb_fcg) icomp1 = icomp1 + nblvlm*2 366 366 icompt = icomp1*(in-1) 367 callrdgrads(99,icompt,nblvlm,z,ht_mes,hq_mes,hw_mes &367 CALL rdgrads(99,icompt,nblvlm,z,ht_mes,hq_mes,hw_mes & 368 368 & ,hu_mes,hv_mes,hthturb_mes,hqturb_mes & 369 369 & ,ts_fcg,ts_subr,imp_fcg,Turb_fcg) … … 518 518 !--------------------------------------------------------------------- 519 519 ! pression au milieu des couches du gcm dans la physiq 520 ! (SB: remplace le callconv_lipress_gcm(playgcm) )520 ! (SB: remplace le CALL conv_lipress_gcm(playgcm) ) 521 521 !--------------------------------------------------------------------- 522 522 … … 532 532 !---------------------------------------------------------------------- 533 533 534 callmesolupbis(file_forctl)534 CALL mesolupbis(file_forctl) 535 535 536 536 print*,'la valeur de nblvlm est:',nblvlm … … 542 542 !---------------------------------------------------------------------- 543 543 544 callcorresbis(psolgcm)544 CALL corresbis(psolgcm) 545 545 546 546 !--------------------------------------------------------- … … 603 603 2000 format (a80) 604 604 aaa=spaces(aaa,1) 605 callgetsch(aaa,' ',' ',2,anblvl,nch)605 CALL getsch(aaa,' ',' ',2,anblvl,nch) 606 606 read(anblvl,*) nblvlm 607 607 -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_lmdz1d.F90
r5099 r5101 5 5 6 6 USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin 7 USE phys_state_var_mod, ONLY 7 USE phys_state_var_mod, ONLY: phys_state_var_init, phys_state_var_end, & 8 8 clwcon, detr_therm, & 9 9 qsol, fevap, z0m, z0h, agesno, & … … 23 23 24 24 USE dimphy 25 USE surface_data, only: type_ocean,ok_veget26 USE pbl_surface_mod, only: ftsoil, pbl_surface_init, &25 USE surface_data, ONLY: type_ocean,ok_veget 26 USE pbl_surface_mod, ONLY: ftsoil, pbl_surface_init, & 27 27 pbl_surface_final 28 USE fonte_neige_mod, only : fonte_neige_init, fonte_neige_final28 USE fonte_neige_mod, ONLY: fonte_neige_init, fonte_neige_final 29 29 30 30 USE infotrac ! new … … 43 43 USE temps_mod, ONLY: annee_ref, calend, day_end, day_ini, day_ref, & 44 44 itau_dyn, itau_phy, start_time, year_len 45 USE phys_cal_mod, ONLY 45 USE phys_cal_mod, ONLY: year_len_phys_cal_mod => year_len 46 46 USE mod_1D_cases_read, ONLY: interp_case_time ! used in included old_1D_read_forc_cases.h 47 47 -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/replay1d.F90
r5099 r5101 4 4 USE inigeomphy_mod, ONLY: inigeomphy 5 5 USE comvert_mod, ONLY: presnivs 6 USE comvert_mod, only: preff, pa6 USE comvert_mod, ONLY: preff, pa 7 7 USE ioipsl, only: getin 8 8 -
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/scm.F90
r5099 r5101 2 2 3 3 USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin 4 USE phys_state_var_mod, ONLY 4 USE phys_state_var_mod, ONLY: phys_state_var_init, phys_state_var_end, & 5 5 clwcon, detr_therm, & 6 6 qsol, fevap, z0m, z0h, agesno, & … … 20 20 21 21 USE dimphy 22 USE surface_data, only: type_ocean,ok_veget23 USE pbl_surface_mod, only: ftsoil, pbl_surface_init, &22 USE surface_data, ONLY: type_ocean,ok_veget 23 USE pbl_surface_mod, ONLY: ftsoil, pbl_surface_init, & 24 24 pbl_surface_final 25 USE fonte_neige_mod, only : fonte_neige_init, fonte_neige_final25 USE fonte_neige_mod, ONLY: fonte_neige_init, fonte_neige_final 26 26 27 27 USE infotrac ! new … … 40 40 USE temps_mod, ONLY: annee_ref, calend, day_end, day_ini, day_ref, & 41 41 itau_dyn, itau_phy, start_time, year_len 42 USE phys_cal_mod, ONLY : year_len_phys_cal_mod => year_len42 USE phys_cal_mod, ONLY: year_len_phys_cal_mod => year_len 43 43 44 44 implicit none -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecumev6_flux.F90
r5099 r5101 86 86 USE dimphy 87 87 USE indice_sol_mod 88 USE MODD_CSTS, ONLY 88 USE MODD_CSTS, ONLY: XPI, XDAY, XKARMAN, XG, XP00, XSTEFAN, XRD, XRV, & 89 89 XCPD, XCPV, XCL, XTT, XLVTT 90 90 91 91 92 !USE MODD_SURF_PAR, ONLY 93 !USE MODD_SURF_ATM, ONLY 94 !USE MODD_REPROD_OPER, ONLY 92 !USE MODD_SURF_PAR, ONLY: XUNDEF 93 !USE MODD_SURF_ATM, ONLY: XVCHRNK, XVZ0CM 94 !USE MODD_REPROD_OPER, ONLY: CCHARNOCK 95 95 96 96 !USE MODE_THERMOS … … 98 98 !USE MODI_SURFACE_RI 99 99 100 !USE YOMHOOK, ONLY 101 !USE PARKIND1, ONLY 100 !USE YOMHOOK, ONLY: LHOOK, DR_HOOK 101 !USE PARKIND1, ONLY: JPRB 102 102 103 103 !USE MODI_ABOR1_SFX -
LMDZ6/branches/Amaury_dev/libf/phylmd/ener_conserv.F90
r5099 r5101 20 20 21 21 ! From module 22 USE phys_local_var_mod, ONLY 22 USE phys_local_var_mod, ONLY: d_u_vdf,d_v_vdf,d_t_vdf,d_u_ajs,d_v_ajs,d_t_ajs, & 23 23 d_u_con,d_v_con,d_t_con,d_t_diss 24 USE phys_local_var_mod, ONLY 25 USE phys_local_var_mod, ONLY 26 USE phys_local_var_mod, ONLY 27 USE phys_state_var_mod, ONLY 28 USE phys_output_var_mod, ONLY 29 USE add_phys_tend_mod, ONLY : fl_cor_ebil24 USE phys_local_var_mod, ONLY: d_t_eva,d_t_lsc,d_q_eva,d_q_lsc 25 USE phys_local_var_mod, ONLY: d_u_oro,d_v_oro,d_u_lif,d_v_lif 26 USE phys_local_var_mod, ONLY: du_gwd_hines,dv_gwd_hines,dv_gwd_front,dv_gwd_rando 27 USE phys_state_var_mod, ONLY: du_gwd_front,du_gwd_rando 28 USE phys_output_var_mod, ONLY: bils_ec,bils_ech,bils_tke,bils_kinetic,bils_enthalp,bils_latent,bils_diss 29 USE add_phys_tend_mod, ONLY: fl_cor_ebil 30 30 USE infotrac_phy, ONLY: nqtot 31 31 -
LMDZ6/branches/Amaury_dev/libf/phylmd/flott_gwd_rando_m.F90
r5099 r5101 21 21 use dimphy, only: klon, klev 22 22 use assert_m, only: assert 23 USE ioipsl_getin_p_mod, ONLY 24 USE vertical_layers_mod, ONLY 23 USE ioipsl_getin_p_mod, ONLY: getin_p 24 USE vertical_layers_mod, ONLY: presnivs 25 25 CHARACTER (LEN=20) :: modname='flott_gwd_rando' 26 26 CHARACTER (LEN=80) :: abort_message -
LMDZ6/branches/Amaury_dev/libf/phylmd/fonte_neige_mod.F90
r5099 r5101 8 8 9 9 !**************************************************************************************** 10 USE dimphy, ONLY 10 USE dimphy, ONLY: klon 11 11 USE indice_sol_mod 12 12 -
LMDZ6/branches/Amaury_dev/libf/phylmd/freinage.F90
r5099 r5101 8 8 use dimphy, only: klon, klev 9 9 ! USE control, ONLY: nvm 10 ! USE indice_sol_mod, only: nvm_orch10 ! USE indice_sol_mod, ONLY: nvm_orch 11 11 12 12 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/grid_index.F90
r5074 r5101 5 5 ! Please do not put this function in a m*odule not to complexify the replay script 6 6 !-------------------------------------------------------------------------------- 7 USE dimphy, only: klon7 USE dimphy, ONLY: klon 8 8 USE geometry_mod, ONLY: latitude_deg, longitude_deg 9 9 implicit none -
LMDZ6/branches/Amaury_dev/libf/phylmd/grid_noro_m.F90
r5100 r5101 436 436 !=============================================================================== 437 437 USE netcdf, ONLY: nf90_open, nf90_inq_dimid, nf90_inquire_dimension, & 438 nf90_noerr, nf90_close, NF90_INQ_VARID, nf90_get_var, NF90_STRERROR, &438 nf90_noerr, nf90_close, nf90_inq_varid, nf90_get_var, nf90_strerror, & 439 439 nf90_nowrite 440 440 IMPLICIT NONE … … 498 498 CHARACTER(LEN=*), INTENT(IN) :: var 499 499 REAL, INTENT(INOUT) :: fld(:,:) 500 CALL ncerr( NF90_INQ_VARID(fid,var,vid),var)500 CALL ncerr(nf90_inq_varid(fid,var,vid),var) 501 501 CALL ncerr(nf90_get_var(fid,vid,fld(1:imar,:)),var) 502 502 fld(imar+1,:)=fld(1,:) … … 511 511 CALL ncerr(nf90_inq_dimid(fid,dimd,did)) 512 512 CALL ncerr(nf90_inquire_dimension(fid,did,len=n)); ALLOCATE(tmp(n)) 513 CALL ncerr( NF90_INQ_VARID(fid,dimd,did))513 CALL ncerr(nf90_inq_varid(fid,dimd,did)) 514 514 CALL ncerr(nf90_get_var(fid,did,tmp)) 515 515 IF(MAXVAL(tmp)>xpi) tmp=tmp*d2r … … 529 529 IF(PRESENT(var)) mess=TRIM(mess)//' and variable "'//TRIM(var)//'"' 530 530 WRITE(lunout,*)TRIM(mess)//'.' 531 CALL abort_physic(modname, NF90_STRERROR(ncres),1)531 CALL abort_physic(modname,nf90_strerror(ncres),1) 532 532 END IF 533 533 END SUBROUTINE ncerr -
LMDZ6/branches/Amaury_dev/libf/phylmd/ice_sursat_mod.F90
r5099 r5101 59 59 60 60 USE print_control_mod, ONLY: lunout 61 USE ioipsl_getin_p_mod, ONLY 61 USE ioipsl_getin_p_mod, ONLY: getin_p 62 62 63 63 IMPLICIT NONE … … 91 91 USE mod_grid_phy_lmdz, ONLY: klon_glo 92 92 USE geometry_mod, ONLY: cell_area 93 USE phys_cal_mod, ONLY 93 USE phys_cal_mod, ONLY: mth_cur 94 94 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 95 95 USE mod_phys_lmdz_omp_data, ONLY: is_omp_root -
LMDZ6/branches/Amaury_dev/libf/phylmd/infotrac_phy.F90
r5099 r5101 3 3 MODULE infotrac_phy 4 4 5 USE strings_mod, ONLY 6 USE readTracFiles_mod, ONLY 5 USE strings_mod, ONLY: msg, fmsg, maxlen, cat, dispTable, int2str, bool2str, strStack, strParse, strIdx 6 USE readTracFiles_mod, ONLY: trac_type, readTracersFiles, tracers, setGeneration, itZonIso, nzone, tran0, isoZone, & 7 7 delPhase, niso, getKey, isot_type, processIsotopes, isotope, maxTableWidth, iqIsoPha, nphas, ixIso, isoPhas, & 8 8 addPhase, iH2O, addKey, isoSelect, testTracersFiles, isoKeys, indexUpdate, iqWIsoPha, nbIso, ntiso, isoName, isoCheck … … 128 128 129 129 SUBROUTINE init_infotrac_phy 130 USE ioipsl_getin_p_mod, ONLY 130 USE ioipsl_getin_p_mod, ONLY: getin_p 131 131 #ifdef REPROBUS 132 132 USE CHEM_REP, ONLY: Init_chem_rep_trac 133 133 #endif 134 USE lmdz_cppkeys_wrapper, ONLY 134 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_STRATAER 135 135 IMPLICIT NONE 136 136 !============================================================================================================================== -
LMDZ6/branches/Amaury_dev/libf/phylmd/ini_csts.F90
r5099 r5101 57 57 58 58 59 !USE YOMHOOK ,ONLY 60 !USE PARKIND1 ,ONLY 59 !USE YOMHOOK ,ONLY: LHOOK, DR_HOOK 60 !USE PARKIND1 ,ONLY: JPRB 61 61 62 62 !USE MODI_INI_CTURBS -
LMDZ6/branches/Amaury_dev/libf/phylmd/inifis_mod.F90
r3435 r5101 6 6 SUBROUTINE inifis(punjours, prad, pg, pr, pcpp) 7 7 ! Initialize some physical constants and settings 8 USE init_print_control_mod, ONLY 8 USE init_print_control_mod, ONLY: init_print_control 9 9 USE print_control_mod, ONLY: lunout 10 10 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/init_be.F90
r5099 r5101 6 6 USE dimphy 7 7 USE indice_sol_mod 8 USE geometry_mod, ONLY 8 USE geometry_mod, ONLY: longitude, latitude 9 9 10 10 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/initrrnpb.F90
r5099 r5101 4 4 SUBROUTINE initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr,vdeptr,scavtr) 5 5 USE dimphy 6 USE infotrac_phy, ONLY 7 USE traclmdz_mod, ONLY 6 USE infotrac_phy, ONLY: nbtr 7 USE traclmdz_mod, ONLY: id_rn, id_pb 8 8 USE indice_sol_mod 9 9 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VAR0SV.F90
r3792 r5101 1 1 MODULE VAR0SV 2 2 3 USE VAR_SV, only : nsol, nsno4 USE VARdSV, only: nsot3 USE VAR_SV, ONLY: nsol, nsno 4 USE VARdSV, ONLY: nsot 5 5 6 6 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARtSV.F90
r3900 r5101 1 1 MODULE VARtSV 2 2 3 USE VAR_SV, only: nsol, nsno, klonv3 USE VAR_SV, ONLY: nsol, nsno, klonv 4 4 5 5 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARxSV.F90
r5099 r5101 1 1 MODULE VARxSV 2 2 3 USE VAR_SV, only: klonv, nsol, nsno, nb_wri3 USE VAR_SV, ONLY: klonv, nsol, nsno, nb_wri 4 4 5 5 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARySV.F90
r5099 r5101 1 1 MODULE VARySV 2 2 3 USE VAR_SV, only: klonv, nsol, nsno3 USE VAR_SV, ONLY: klonv, nsol, nsno 4 4 5 5 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90
r5099 r5101 50 50 USE VARtSV 51 51 USE VARphy 52 USE surface_data, only: iflag_tsurf_inlandsis, SnoMod, BloMod, ok_outfor52 USE surface_data, ONLY: iflag_tsurf_inlandsis, SnoMod, BloMod, ok_outfor 53 53 54 54 IMPLICIT NONE … … 751 751 ! of SISVAT, therefore it's needed here. 752 752 753 USE mod_phys_lmdz_mpi_data, ONLY 753 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 754 754 USE mod_phys_lmdz_para 755 755 USE VAR_SV -
LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90
r5099 r5101 166 166 mpi_size, mpi_rank 167 167 USE ioipsl, ONLY: flio_dom_set 168 USE mod_grid_phy_lmdz, ONLY 168 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 169 169 170 170 IMPLICIT NONE … … 211 211 USE mod_phys_lmdz_para, ONLY: is_sequential, is_using_mpi, is_mpi_root, & 212 212 jj_begin, jj_end, jj_nb 213 USE mod_grid_phy_lmdz, ONLY 213 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 214 214 USE ioipsl, ONLY: histbeg 215 215 USE wxios, ONLY: wxios_add_file, using_xios … … 251 251 252 252 USE mod_phys_lmdz_para, ONLY: jj_begin, jj_end, jj_nb, is_sequential 253 USE mod_grid_phy_lmdz, ONLY 253 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 254 254 USE ioipsl, ONLY: histbeg 255 255 … … 461 461 nid_files, nhorim, swaero_diag, dryaod_diag, nfiles, & 462 462 ok_4xCO2atm 463 USE mod_grid_phy_lmdz, ONLY 464 USE aero_mod, ONLY 463 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 464 USE aero_mod, ONLY: naero_tot, name_aero_tau 465 465 USE print_control_mod, ONLY: prt_level,lunout 466 466 … … 538 538 nhorim, zdtime_moy, levmin, levmax, & 539 539 nvertm, nfiles 540 USE mod_grid_phy_lmdz, ONLY 540 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 541 541 USE print_control_mod, ONLY: prt_level,lunout 542 542 IMPLICIT NONE … … 590 590 ok_4xCO2atm 591 591 USE print_control_mod, ONLY: prt_level,lunout 592 USE mod_grid_phy_lmdz, ONLY 593 USE aero_mod, ONLY 592 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 593 USE aero_mod, ONLY: naero_tot, name_aero_tau 594 594 USE wxios, ONLY: wxios_add_field_to_file, using_xios 595 595 USE print_control_mod, ONLY: prt_level,lunout … … 706 706 levmax, nvertm 707 707 USE print_control_mod, ONLY: prt_level,lunout 708 USE mod_grid_phy_lmdz, ONLY 708 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 709 709 USE wxios, ONLY: wxios_add_field_to_file, using_xios 710 710 USE print_control_mod, ONLY: prt_level,lunout … … 811 811 USE ioipsl, ONLY: histwrite 812 812 USE print_control_mod, ONLY: prt_level,lunout 813 USE mod_grid_phy_lmdz, ONLY 813 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 814 814 815 815 IMPLICIT NONE … … 876 876 is_sequential, klon_mpi_begin, klon_mpi_end, & 877 877 jj_nb, klon_mpi, is_master 878 USE mod_grid_phy_lmdz, ONLY 878 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 879 879 USE ioipsl, ONLY: histwrite 880 880 USE print_control_mod, ONLY: prt_level,lunout … … 954 954 nid_files, swaerofree_diag, swaero_diag, dryaod_diag, ok_4xCO2atm 955 955 USE print_control_mod, ONLY: prt_level,lunout 956 USE mod_grid_phy_lmdz, ONLY 956 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, unstructured, regular_lonlat 957 957 USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios 958 958 USE print_control_mod, ONLY: lunout, prt_level … … 1166 1166 nfiles, vars_defined, clef_stations, & 1167 1167 nid_files, swaerofree_diag 1168 USE mod_grid_phy_lmdz, ONLY 1168 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, regular_lonlat, unstructured 1169 1169 USE lmdz_xios, ONLY: xios_send_field, xios_field_is_active, using_xios 1170 1170 USE print_control_mod, ONLY: prt_level,lunout … … 1362 1362 is_sequential, klon_mpi_begin, klon_mpi_end, & 1363 1363 jj_nb, klon_mpi, is_master 1364 USE mod_grid_phy_lmdz, ONLY 1364 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, unstructured 1365 1365 USE lmdz_xios, ONLY: xios_send_field 1366 1366 USE print_control_mod, ONLY: prt_level,lunout … … 1446 1446 jj_nb, klon_mpi, is_master 1447 1447 USE lmdz_xios, ONLY: xios_send_field 1448 USE mod_grid_phy_lmdz, ONLY 1448 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, unstructured 1449 1449 USE print_control_mod, ONLY: prt_level,lunout 1450 1450 -
LMDZ6/branches/Amaury_dev/libf/phylmd/iophys.F90
r5099 r5101 77 77 78 78 USE mod_phys_lmdz_para, ONLY: klon_omp 79 USE dimphy, ONLY 79 USE dimphy, ONLY: klon 80 80 USE mod_grid_phy_lmdz, ONLY: klon_glo 81 81 IMPLICIT NONE … … 115 115 USE dimphy, ONLY: klev 116 116 USE mod_grid_phy_lmdz, ONLY: klon_glo 117 USE time_phylmdz_mod, ONLY 118 USE phys_cal_mod, ONLY 117 USE time_phylmdz_mod, ONLY: annee_ref, day_ref, day_ini 118 USE phys_cal_mod, ONLY: calend 119 119 120 120 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/iostart.F90
r5100 r5101 70 70 71 71 IF (is_mpi_root .AND. is_omp_root) THEN 72 ierr= NF90_INQ_VARID(nid_start,Field_name,varid)72 ierr=nf90_inq_varid(nid_start,Field_name,varid) 73 73 IF (ierr==nf90_noerr) THEN 74 74 Inquire_field=.TRUE. … … 146 146 IF (is_master) THEN 147 147 148 ierr= NF90_INQ_VARID(nid_start,Field_name,varid)148 ierr=nf90_inq_varid(nid_start,Field_name,varid) 149 149 150 150 IF (ierr==nf90_noerr) THEN … … 173 173 IF (.NOT. tmp_found) THEN 174 174 PRINT*, 'phyetat0: Le champ <'//field_name//'> est absent' 175 callabort_physic("", "", 1)175 CALL abort_physic("", "", 1) 176 176 ENDIF 177 177 ENDIF … … 193 193 IF (ierr/=nf90_noerr) THEN 194 194 PRINT*, 'phyetat0: Lecture echouee aussi en 2D pour <'//field_name//'>' 195 callabort_physic("", "", 1)195 CALL abort_physic("", "", 1) 196 196 ELSE 197 197 PRINT*, 'phyetat0: La variable <'//field_name//'> lu sur surface seulement'!, selon ancien format, le reste mis a zero' 198 198 END IF 199 199 ELSE 200 callabort_physic("", "", 1)200 CALL abort_physic("", "", 1) 201 201 ENDIF 202 202 ENDIF … … 267 267 IF (is_mpi_root .AND. is_omp_root) THEN 268 268 269 ierr= NF90_INQ_VARID(nid_start,var_name,varid)269 ierr=nf90_inq_varid(nid_start,var_name,varid) 270 270 271 271 IF (ierr==nf90_noerr) THEN … … 273 273 IF (ierr/=nf90_noerr) THEN 274 274 PRINT*, 'phyetat0: Lecture echouee pour <'//var_name//'>' 275 callabort_physic("", "", 1)275 CALL abort_physic("", "", 1) 276 276 ENDIF 277 277 tmp_found=.TRUE. … … 293 293 IF (.NOT. tmp_found) THEN 294 294 PRINT*, 'phyetat0: La variable champ <'//var_name//'> est absente' 295 callabort_physic("", "", 1)295 CALL abort_physic("", "", 1) 296 296 ENDIF 297 297 ENDIF … … 425 425 ENDIF 426 426 427 ! ierr = NF90_REDEF(nid_restart)427 ! ierr = nf90_redef (nid_restart) 428 428 ierr = nf90_def_var (nid_restart, field_name, nf90_format,(/ idim /),nvarid) 429 429 IF (LEN_TRIM(title) > 0) ierr = nf90_put_att (nid_restart,nvarid,"title", title) … … 453 453 ENDDO 454 454 455 ierr = NF90_INQ_VARID(nid_restart, field_name, nvarid)455 ierr = nf90_inq_varid(nid_restart, field_name, nvarid) 456 456 ierr = nf90_put_var(nid_restart,nvarid,RESHAPE(field_glo,(/klon_glo*field_size/))) 457 457 ENDIF … … 528 528 IF (var_size/=length) THEN 529 529 PRINT *, "erreur phyredem : probleme de dimension" 530 callabort_physic("", "", 1)530 CALL abort_physic("", "", 1) 531 531 ENDIF 532 532 … … 534 534 IF (pass==1) THEN 535 535 536 ! ierr = NF90_REDEF(nid_restart)536 ! ierr = nf90_redef (nid_restart) 537 537 538 538 ierr = nf90_def_var (nid_restart, var_name, nf90_format,(/ idim1 /),nvarid) … … 542 542 ! second pass : write 543 543 ELSE IF (pass==2) THEN 544 ierr = NF90_INQ_VARID(nid_restart, var_name, nvarid)544 ierr = nf90_inq_varid(nid_restart, var_name, nvarid) 545 545 ierr = nf90_put_var(nid_restart,nvarid,var) 546 546 ENDIF -
LMDZ6/branches/Amaury_dev/libf/phylmd/iotd_fin.F90
r5099 r5101 1 1 SUBROUTINE iotd_fin 2 USE netcdf, ONLY 2 USE netcdf, ONLY: nf90_close 3 3 4 4 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/limit_read_mod.F90
r5100 r5101 131 131 ! This subroutine returns the sea surface temperature already read from limit.nc. 132 132 133 USE dimphy, ONLY 133 USE dimphy, ONLY: klon 134 134 135 135 INTEGER, INTENT(IN) :: knon ! nomber of points on compressed grid … … 164 164 USE mod_grid_phy_lmdz 165 165 USE mod_phys_lmdz_para 166 USE surface_data, ONLY 166 USE surface_data, ONLY: type_ocean, ok_veget 167 167 USE netcdf, ONLY:nf90_get_var,nf90_inq_varid,nf90_close,nf90_inquire_dimension,& 168 168 nf90_inquire,nf90_get_att,nf90_inq_dimid,nf90_nowrite,nf90_noerr,nf90_open 169 169 USE indice_sol_mod 170 USE phys_cal_mod, ONLY 170 USE phys_cal_mod, ONLY: calend, year_len 171 171 USE print_control_mod, ONLY: lunout, prt_level 172 172 USE lmdz_xios, ONLY: xios_recv_field, using_xios … … 239 239 240 240 !--- WARNING IF CALENDAR IS KNOWN AND DOES NOT MATCH THE ONE OF LMDZ 241 ierr= NF90_INQ_VARID(nid, 'TEMPS', nvarid)241 ierr=nf90_inq_varid(nid, 'TEMPS', nvarid) 242 242 ierr=nf90_get_att(nid, nvarid, 'calendar', calendar) 243 243 IF(ierr==nf90_noerr.AND.calendar/=calend.AND.prt_level>=1) THEN … … 251 251 ierr=nf90_inq_dimid(nid,"time_year",ndimid) 252 252 ELSE 253 ierr= NF90_INQUIRE(nid, UnlimitedDimID=ndimid)253 ierr=nf90_inquire(nid, UnlimitedDimID=ndimid) 254 254 ENDIF 255 255 ierr=nf90_inquire_dimension(nid, ndimid, len=nn) … … 351 351 352 352 ! Ocean fraction 353 ierr = NF90_INQ_VARID(nid, 'FOCE', nvarid)353 ierr = nf90_inq_varid(nid, 'FOCE', nvarid) 354 354 IF (ierr /= nf90_noerr) CALL abort_physic(modname, 'Le champ <FOCE> est absent',1) 355 355 … … 358 358 359 359 ! Sea-ice fraction 360 ierr = NF90_INQ_VARID(nid, 'FSIC', nvarid)360 ierr = nf90_inq_varid(nid, 'FSIC', nvarid) 361 361 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <FSIC> est absent',1) 362 362 … … 369 369 370 370 ! Land fraction 371 ierr = NF90_INQ_VARID(nid, 'FTER', nvarid)371 ierr = nf90_inq_varid(nid, 'FTER', nvarid) 372 372 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <FTER> est absent',1) 373 373 … … 376 376 377 377 ! Continentale ice fraction 378 ierr = NF90_INQ_VARID(nid, 'FLIC', nvarid)378 ierr = nf90_inq_varid(nid, 'FLIC', nvarid) 379 379 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <FLIC> est absent',1) 380 380 … … 391 391 IF ( type_ocean /= 'couple') THEN 392 392 393 ierr = NF90_INQ_VARID(nid, 'SST', nvarid)393 ierr = nf90_inq_varid(nid, 'SST', nvarid) 394 394 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <SST> est absent',1) 395 395 … … 407 407 408 408 ! Read albedo 409 ierr = NF90_INQ_VARID(nid, 'ALB', nvarid)409 ierr = nf90_inq_varid(nid, 'ALB', nvarid) 410 410 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <ALB> est absent',1) 411 411 … … 414 414 415 415 ! Read rugosity 416 ierr = NF90_INQ_VARID(nid, 'RUG', nvarid)416 ierr = nf90_inq_varid(nid, 'RUG', nvarid) 417 417 IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Le champ <RUG> est absent',1) 418 418 -
LMDZ6/branches/Amaury_dev/libf/phylmd/limit_slab.F90
r5100 r5101 84 84 ! First read first layer 85 85 ! try first "BILS_OCE01" 86 ierr = NF90_INQ_VARID(nid, 'BILS_OCE01', nvarid)86 ierr = nf90_inq_varid(nid, 'BILS_OCE01', nvarid) 87 87 IF (ierr /= nf90_noerr) THEN 88 88 ! Else BILS_OCE 89 ierr = NF90_INQ_VARID(nid, 'BILS_OCE', nvarid)89 ierr = nf90_inq_varid(nid, 'BILS_OCE', nvarid) 90 90 IF (ierr /= nf90_noerr) THEN 91 91 read_bils=.FALSE. … … 102 102 DO i=2,nslay 103 103 WRITE(str2,'(i2.2)') i 104 ierr = NF90_INQ_VARID(nid,'BILS_OCE'//str2, nvarid)104 ierr = nf90_inq_varid(nid,'BILS_OCE'//str2, nvarid) 105 105 IF (ierr==nf90_noerr) THEN 106 106 ierr = nf90_get_var(nid,nvarid,bils_glo(:,i),start,epais) … … 113 113 114 114 ! Read sst_glo for this day 115 ierr = NF90_INQ_VARID(nid, 'SST', nvarid)115 ierr = nf90_inq_varid(nid, 'SST', nvarid) 116 116 IF (ierr /= nf90_noerr) THEN 117 117 read_sst=.FALSE. … … 127 127 128 128 ! Read siv_glo for this day 129 ierr = NF90_INQ_VARID(nid, 'SICV', nvarid)129 ierr = nf90_inq_varid(nid, 'SICV', nvarid) 130 130 IF (ierr /= nf90_noerr) THEN 131 131 read_siv=.FALSE. -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_atke_exchange_coeff.F90
r5099 r5101 23 23 !======================================================================= 24 24 25 USE lmdz_atke_turbulence_ini, ONLY 26 USE lmdz_atke_turbulence_ini, ONLY 27 USE lmdz_atke_turbulence_ini, ONLY 25 USE lmdz_atke_turbulence_ini, ONLY: iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual, ri0, ri1 26 USE lmdz_atke_turbulence_ini, ONLY: cepsilon, pr_slope, pr_asym, pr_neut, ctkes, rg, rd, rv, atke_ok_vdiff 27 USE lmdz_atke_turbulence_ini, ONLY: viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin, cn 28 28 29 29 !!------------------------------------------------------------------------------------------------------------- … … 416 416 ! E Vignon, July 2023 417 417 418 USE lmdz_atke_turbulence_ini, ONLY 418 USE lmdz_atke_turbulence_ini, ONLY: rd, cke, viscom 419 419 420 420 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_atke_turbulence_ini.F90
r5082 r5101 53 53 !!---------------------------------------------------------------------- 54 54 55 USE ioipsl_getin_p_mod, ONLY : getin_p55 USE ioipsl_getin_p_mod, ONLY: getin_p 56 56 57 57 ! input arguments (universal constants for planet) -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_blowing_snow_ini.F90
r4835 r5101 34 34 RVTMP2_in, RTT_in,RD_in,RG_in, RV_in, RPI_in) 35 35 36 USE ioipsl_getin_p_mod, ONLY 36 USE ioipsl_getin_p_mod, ONLY: getin_p 37 37 38 38 real, intent(in) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, RPI_in -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90
r5099 r5101 11 11 12 12 13 use lmdz_blowing_snow_ini, only: iflag_sublim_bs, iflag_sedim_bs, coef_sub_bs,RTT,RD,RG,fallv_bs14 use lmdz_blowing_snow_ini, only: qbmin, RCPD, RLSTT, RLMLT, RLVTT, RVTMP2, RV, RPI, tbsmelt, taumeltbs0, rhobs, r_bs15 USE lmdz_lscp_tools, only: calc_qsat_ecmwf13 use lmdz_blowing_snow_ini, ONLY: iflag_sublim_bs, iflag_sedim_bs, coef_sub_bs,RTT,RD,RG,fallv_bs 14 use lmdz_blowing_snow_ini, ONLY: qbmin, RCPD, RLSTT, RLMLT, RLVTT, RVTMP2, RV, RPI, tbsmelt, taumeltbs0, rhobs, r_bs 15 USE lmdz_lscp_tools, ONLY: calc_qsat_ecmwf 16 16 17 17 implicit none -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_call_atke.F90
r5082 r5101 1 1 module lmdz_call_atke 2 2 3 USE lmdz_atke_exchange_coeff, ONLY 3 USE lmdz_atke_exchange_coeff, ONLY: atke_compute_km_kh 4 4 5 5 implicit none … … 15 15 16 16 17 USE lmdz_atke_turbulence_ini, ONLY 17 USE lmdz_atke_turbulence_ini, ONLY: iflag_num_atke, rg, rd 18 18 USE phys_local_var_mod, ONLY: tke_shear, tke_buoy, tke_trans 19 19 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_call_blowing_snow.F90
r4724 r5101 6 6 dtemp_bs,dq_bs,dqbs_bs,bsfl,precip_bs) 7 7 8 use lmdz_blowing_snow_sublim_sedim, only: blowing_snow_sublim_sedim8 use lmdz_blowing_snow_sublim_sedim, ONLY: blowing_snow_sublim_sedim 9 9 implicit none 10 10 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_call_cloud_optics_prop.F90
r5087 r5101 18 18 !------------------ 19 19 20 USE lmdz_cloud_optics_prop_ini , ONLY 21 USE lmdz_cloud_optics_prop_ini , ONLY 22 USE lmdz_cloud_optics_prop_ini , ONLY 23 USE lmdz_cloud_optics_prop, ONLY 20 USE lmdz_cloud_optics_prop_ini , ONLY: ok_cdnc 21 USE lmdz_cloud_optics_prop_ini , ONLY: bl95_b0, bl95_b1 22 USE lmdz_cloud_optics_prop_ini , ONLY: iflag_ice_thermo, ok_new_lscp, iflag_t_glace 23 USE lmdz_cloud_optics_prop, ONLY: cloud_optics_prop 24 24 25 25 USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14) -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_cloud_optics_prop.F90
r5099 r5101 11 11 icefrac_optics, dNovrN, ptconv,rnebcon, ccwcon) 12 12 13 USE lmdz_cloud_optics_prop_ini , ONLY 14 USE lmdz_cloud_optics_prop_ini , ONLY 15 USE lmdz_cloud_optics_prop_ini , ONLY 16 USE lmdz_cloud_optics_prop_ini , ONLY 17 USE lmdz_cloud_optics_prop_ini , ONLY 18 USE lmdz_cloud_optics_prop_ini , ONLY 19 USE lmdz_cloud_optics_prop_ini , ONLY 20 USE lmdz_cloud_optics_prop_ini , ONLY 21 USE lmdz_cloud_optics_prop_ini , ONLY 22 USE lmdz_cloud_optics_prop_ini , ONLY 23 USE lmdz_cloud_optics_prop_ini , ONLY 24 USE lmdz_cloud_optics_prop_ini , ONLY 25 USE lmdz_cloud_optics_prop_ini , ONLY 26 USE lmdz_cloud_optics_prop_ini , ONLY 27 USE lmdz_cloud_optics_prop_ini , ONLY 28 USE lmdz_cloud_optics_prop_ini , ONLY 29 USE lmdz_cloud_optics_prop_ini , ONLY 13 USE lmdz_cloud_optics_prop_ini , ONLY: flag_aerosol, ok_cdnc 14 USE lmdz_cloud_optics_prop_ini , ONLY: lunout 15 USE lmdz_cloud_optics_prop_ini , ONLY: bl95_b0, bl95_b1 16 USE lmdz_cloud_optics_prop_ini , ONLY: latitude_deg 17 USE lmdz_cloud_optics_prop_ini , ONLY: iflag_t_glace 18 USE lmdz_cloud_optics_prop_ini , ONLY: cdnc_max, cdnc_max_m3 19 USE lmdz_cloud_optics_prop_ini , ONLY: cdnc_min, cdnc_min_m3 20 USE lmdz_cloud_optics_prop_ini , ONLY: thres_tau, thres_neb 21 USE lmdz_cloud_optics_prop_ini , ONLY: prmhc, prlmc 22 USE lmdz_cloud_optics_prop_ini , ONLY: coef_froi, coef_chau 23 USE lmdz_cloud_optics_prop_ini , ONLY: seuil_neb 24 USE lmdz_cloud_optics_prop_ini , ONLY: t_glace_min_old, t_glace_max_old 25 USE lmdz_cloud_optics_prop_ini , ONLY: k_ice0, df 26 USE lmdz_cloud_optics_prop_ini , ONLY: rg, rd, rpi 27 USE lmdz_cloud_optics_prop_ini , ONLY: rad_chau1, rad_chau2, rad_froid, iflag_rei 28 USE lmdz_cloud_optics_prop_ini , ONLY: ok_icefra_lscp, rei_max, rei_min 29 USE lmdz_cloud_optics_prop_ini , ONLY: zepsec, novlp, iflag_ice_thermo, ok_new_lscp 30 30 31 31 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_cloud_optics_prop_ini.F90
r5087 r5101 48 48 iflag_ice_thermo_in, ok_new_lscp_in) 49 49 50 USE ioipsl_getin_p_mod, ONLY 50 USE ioipsl_getin_p_mod, ONLY: getin_p 51 51 52 52 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_cloudth.F90
r5099 r5101 819 819 !=========================================================================== 820 820 821 use lmdz_cloudth_ini, only: iflag_cloudth_vert,iflag_ratqs822 use lmdz_cloudth_ini, only: vert_alpha,vert_alpha_th, sigma1s_factor, sigma1s_power , sigma2s_factor , sigma2s_power , cloudth_ratqsmin , iflag_cloudth_vert_noratqs821 use lmdz_cloudth_ini, ONLY: iflag_cloudth_vert,iflag_ratqs 822 use lmdz_cloudth_ini, ONLY: vert_alpha,vert_alpha_th, sigma1s_factor, sigma1s_power , sigma2s_factor , sigma2s_power , cloudth_ratqsmin , iflag_cloudth_vert_noratqs 823 823 824 824 IMPLICIT NONE … … 2123 2123 !============================================================================= 2124 2124 2125 use phys_state_var_mod, ONLY 2125 use phys_state_var_mod, ONLY: fm_therm, detr_therm, entr_therm 2126 2126 2127 2127 IMPLICIT none -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_cloudth_ini.F90
r4651 r5101 12 12 subroutine cloudth_ini(iflag_cloudth_vert_in,iflag_ratqs_in) 13 13 14 use ioipsl_getin_p_mod, ONLY 14 use ioipsl_getin_p_mod, ONLY: getin_p 15 15 implicit none 16 16 integer, intent(in) :: iflag_cloudth_vert_in,iflag_ratqs_in -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp.F90
r5099 r5101 95 95 96 96 ! USE de modules contenant des fonctions. 97 USE lmdz_cloudth, ONLY : cloudth, cloudth_v3, cloudth_v6, cloudth_mpc98 USE lmdz_lscp_tools, ONLY 99 USE lmdz_lscp_tools, ONLY 100 USE lmdz_lscp_tools, ONLY 101 USE ice_sursat_mod, ONLY 102 USE lmdz_lscp_poprecip, ONLY 97 USE lmdz_cloudth, ONLY: cloudth, cloudth_v3, cloudth_v6, cloudth_mpc 98 USE lmdz_lscp_tools, ONLY: calc_qsat_ecmwf, calc_gammasat 99 USE lmdz_lscp_tools, ONLY: icefrac_lscp, icefrac_lscp_turb 100 USE lmdz_lscp_tools, ONLY: fallice_velocity, distance_to_cloud_top 101 USE ice_sursat_mod, ONLY: ice_sursat 102 USE lmdz_lscp_poprecip, ONLY: poprecip_precld, poprecip_postcld 103 103 104 104 ! Use du module lmdz_lscp_ini contenant les constantes 105 USE lmdz_lscp_ini, ONLY 106 USE lmdz_lscp_ini, ONLY 107 USE lmdz_lscp_ini, ONLY 108 USE lmdz_lscp_ini, ONLY : iflag_cloudth_vert, iflag_rain_incloud_vol, iflag_t_glace, t_glace_min109 USE lmdz_lscp_ini, ONLY 110 USE lmdz_lscp_ini, ONLY 111 USE lmdz_lscp_ini, ONLY 112 USE lmdz_lscp_ini, ONLY 113 USE lmdz_lscp_ini, ONLY 114 USE lmdz_lscp_ini, ONLY 105 USE lmdz_lscp_ini, ONLY: prt_level, lunout 106 USE lmdz_lscp_ini, ONLY: seuil_neb, niter_lscp, iflag_evap_prec, t_coup, DDT0, ztfondue, rain_int_min 107 USE lmdz_lscp_ini, ONLY: ok_radocond_snow, a_tr_sca, cld_expo_con, cld_expo_lsc 108 USE lmdz_lscp_ini, ONLY: iflag_cloudth_vert, iflag_rain_incloud_vol, iflag_t_glace, t_glace_min 109 USE lmdz_lscp_ini, ONLY: coef_eva, coef_sub,cld_tau_lsc, cld_tau_con, cld_lc_lsc, cld_lc_con 110 USE lmdz_lscp_ini, ONLY: iflag_bergeron, iflag_fisrtilp_qsat, iflag_vice, cice_velo, dice_velo 111 USE lmdz_lscp_ini, ONLY: iflag_autoconversion, ffallv_con, ffallv_lsc, min_frac_th_cld 112 USE lmdz_lscp_ini, ONLY: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG 113 USE lmdz_lscp_ini, ONLY: ok_poprecip 114 USE lmdz_lscp_ini, ONLY: iflag_icefrac 115 115 IMPLICIT NONE 116 116 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_ini.F90
r5082 r5101 251 251 252 252 253 USE ioipsl_getin_p_mod, ONLY 253 USE ioipsl_getin_p_mod, ONLY: getin_p 254 254 USE ice_sursat_mod, ONLY: ice_sursat_init 255 USE lmdz_cloudth_ini, ONLY 255 USE lmdz_cloudth_ini, ONLY: cloudth_ini 256 256 257 257 REAL, INTENT(IN) :: dtime -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_old.F90
r5099 r5101 14 14 15 15 USE icefrac_lsc_mod ! compute ice fraction (JBM 3/14) 16 USE lmdz_cloudth, only: cloudth, cloudth_v3, cloudth_v616 USE lmdz_cloudth, ONLY: cloudth, cloudth_v3, cloudth_v6 17 17 18 18 USE lmdz_lscp_ini, ONLY: prt_level, lunout 19 USE lmdz_lscp_ini, ONLY : fl_cor_ebil19 USE lmdz_lscp_ini, ONLY: fl_cor_ebil 20 20 USE lmdz_lscp_ini, ONLY: iflag_t_glace,t_glace_min, t_glace_max, exposant_glace 21 USE lmdz_lscp_ini, ONLY 21 USE lmdz_lscp_ini, ONLY: seuil_neb, rain_int_min, iflag_evap_prec, iflag_oldbug_fisrtilp,a_tr_sca 22 22 USE lmdz_lscp_ini, ONLY: iflag_cloudth_vert, iflag_rain_incloud_vol 23 23 USE lmdz_lscp_ini, ONLY: coef_eva, ffallv_lsc, ffallv_con -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_poprecip.F90
r5099 r5101 22 22 ) 23 23 24 USE lmdz_lscp_ini, ONLY 25 USE lmdz_lscp_ini, ONLY 26 USE lmdz_lscp_ini, ONLY 27 USE lmdz_lscp_ini, ONLY 28 USE lmdz_lscp_tools, ONLY 24 USE lmdz_lscp_ini, ONLY: prt_level, lunout 25 USE lmdz_lscp_ini, ONLY: coef_eva, coef_sub, expo_eva, expo_sub, thresh_precip_frac 26 USE lmdz_lscp_ini, ONLY: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG 27 USE lmdz_lscp_ini, ONLY: ok_corr_vap_evasub 28 USE lmdz_lscp_tools, ONLY: calc_qsat_ecmwf 29 29 30 30 IMPLICIT NONE … … 280 280 dqsauto, dqsagg, dqsrim, dqsmelt, dqsfreez) 281 281 282 USE lmdz_lscp_ini, ONLY 283 USE lmdz_lscp_ini, ONLY 284 USE lmdz_lscp_tools, ONLY 285 286 USE lmdz_lscp_ini, ONLY 282 USE lmdz_lscp_ini, ONLY: prt_level, lunout 283 USE lmdz_lscp_ini, ONLY: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RPI 284 USE lmdz_lscp_tools, ONLY: calc_qsat_ecmwf 285 286 USE lmdz_lscp_ini, ONLY: cld_lc_con, cld_tau_con, cld_expo_con, seuil_neb, & 287 287 cld_lc_lsc, cld_tau_lsc, cld_expo_lsc, rain_int_min, & 288 288 thresh_precip_frac, gamma_col, gamma_agg, gamma_rim, & -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_tools.F90
r5082 r5101 114 114 USE print_control_mod, ONLY: lunout, prt_level 115 115 USE lmdz_lscp_ini, ONLY: t_glace_min, t_glace_max, exposant_glace, iflag_t_glace 116 USE lmdz_lscp_ini, ONLY 116 USE lmdz_lscp_ini, ONLY: RTT, dist_liq, temp_nowater 117 117 118 118 IMPLICIT NONE … … 247 247 248 248 249 USE lmdz_lscp_ini, ONLY 250 USE lmdz_lscp_ini, ONLY 251 USE lmdz_lscp_ini, ONLY 252 USE lmdz_lscp_ini, ONLY 253 USE lmdz_lscp_ini, ONLY 249 USE lmdz_lscp_ini, ONLY: prt_level, lunout 250 USE lmdz_lscp_ini, ONLY: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG, RV, RPI 251 USE lmdz_lscp_ini, ONLY: seuil_neb, temp_nowater 252 USE lmdz_lscp_ini, ONLY: tau_mixenv, lmix_mpc, naero5, gamma_snwretro, gamma_taud, capa_crystal 253 USE lmdz_lscp_ini, ONLY: eps 254 254 255 255 IMPLICIT NONE … … 651 651 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 652 652 653 USE lmdz_lscp_ini, ONLY 653 USE lmdz_lscp_ini, ONLY: rd,rg,tresh_cl 654 654 655 655 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_ratqs_ini.F90
r5081 r5101 25 25 SUBROUTINE ratqs_ini(klon,klev,iflag_thermals,lunout_in,nbsrf_in,is_lic_in,is_ter_in,RG_in,RV_in,RD_in,RCPD_in,RLSTT_in,RLVTT_in,RTT_in) 26 26 27 USE ioipsl_getin_p_mod, ONLY 27 USE ioipsl_getin_p_mod, ONLY: getin_p 28 28 29 29 integer, intent(in) :: klon,klev,iflag_thermals,lunout_in -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_ratqs_multi.F90
r5099 r5101 30 30 ratqs_inter_) 31 31 32 USE lmdz_ratqs_ini, ONLY 33 USE lmdz_ratqs_ini, ONLY 34 USE lmdz_ratqs_ini, ONLY 35 USE lmdz_thermcell_dq, ONLY 32 USE lmdz_ratqs_ini, ONLY: a_ratqs_cv,tau_var,fac_tau,tau_cumul,a_ratqs_wake, dqimpl 33 USE lmdz_ratqs_ini, ONLY: RG 34 USE lmdz_ratqs_ini, ONLY: povariance, var_conv 35 USE lmdz_thermcell_dq, ONLY: thermcell_dq 36 36 37 37 implicit none … … 214 214 ! Etienne Vignon, November 2021: effect of subgrid orography on ratqs 215 215 216 USE lmdz_ratqs_ini, ONLY 216 USE lmdz_ratqs_ini, ONLY: RG,RV,RD,RLSTT,RLVTT,RTT,nbsrf,is_lic,is_ter 217 217 218 218 IMPLICIT NONE … … 296 296 USE lmdz_lscp_tools, ONLY: CALC_QSAT_ECMWF 297 297 298 USE lmdz_ratqs_ini, ONLY 298 USE lmdz_ratqs_ini, ONLY: RG,RD,RTT,nbsrf 299 299 300 300 IMPLICIT NONE … … 377 377 ! model. J. Atmos. Sci. 73, 279???291. https://doi.org/10.1175/JAS-D-15-0021. 378 378 379 USE lmdz_ratqs_ini, ONLY 379 USE lmdz_ratqs_ini, ONLY: RG,RV,RD,RCPD,RLSTT,RLVTT,RTT 380 380 381 381 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_alp.F90
r5099 r5101 20 20 21 21 USE indice_sol_mod 22 USE lmdz_thermcell_main, ONLY 22 USE lmdz_thermcell_main, ONLY: thermcell_tke_transport 23 23 IMPLICIT NONE 24 24 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_down.F90
r5099 r5101 226 226 227 227 228 USE lmdz_thermcell_ini, ONLY 228 USE lmdz_thermcell_ini, ONLY: prt_level,RLvCp,RKAPPA,RETV,fact_thermals_down 229 229 IMPLICIT NONE 230 230 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_env.F90
r5099 r5101 11 11 12 12 13 USE lmdz_thermcell_ini, ONLY 14 USE lmdz_thermcell_qsat, ONLY 13 USE lmdz_thermcell_ini, ONLY: prt_level,RLvCp,RKAPPA,RETV 14 USE lmdz_thermcell_qsat, ONLY: thermcell_qsat 15 15 IMPLICIT NONE 16 16 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_flux2.F90
r5099 r5101 16 16 !--------------------------------------------------------------------------- 17 17 18 USE lmdz_thermcell_ini, ONLY 18 USE lmdz_thermcell_ini, ONLY: prt_level,iflag_thermals_optflux 19 19 IMPLICIT NONE 20 20 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_ini.F90
r5087 r5101 50 50 RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in) 51 51 52 USE ioipsl_getin_p_mod, ONLY 52 USE ioipsl_getin_p_mod, ONLY: getin_p 53 53 54 54 integer, intent(in) :: iflag_thermals,prt_level_in,lunout_in -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_main.F90
r5099 r5101 42 42 43 43 ! USE necessaires pour les lignes importees de thermcell_env 44 USE lmdz_thermcell_ini, ONLY 45 USE lmdz_thermcell_qsat, ONLY 44 USE lmdz_thermcell_ini, ONLY: RLvCp,RKAPPA,RETV 45 USE lmdz_thermcell_qsat, ONLY: thermcell_qsat 46 46 47 47 48 48 #ifdef ISO 49 USE infotrac_phy, ONLY 49 USE infotrac_phy, ONLY: ntiso 50 50 #ifdef ISOVERIF 51 USE isotopes_mod, ONLY 51 USE isotopes_mod, ONLY: iso_eau,iso_HDO 52 52 USE isotopes_verif_mod, ONLY: iso_verif_egalite, & 53 53 iso_verif_aberrant_encadre -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_old.F90
r5093 r5101 8 8 USE dimphy 9 9 USE write_field_phy 10 USE lmdz_thermcell_dv2, ONLY 11 USE lmdz_thermcell_dq, ONLY 10 USE lmdz_thermcell_dv2, ONLY: thermcell_dv2 11 USE lmdz_thermcell_dq, ONLY: thermcell_dq 12 12 IMPLICIT NONE 13 13 … … 1012 1012 ! si sature: ql est surestime, d'ou la sous-relax 1013 1013 dt = 0.5*rlvcp*qlbef 1014 ! write(18,*) ,'DT0=',DT1014 ! write(18,*) 'DT0=',DT 1015 1015 ! on pourra enchainer 2 ou 3 calculs sans Do while 1016 1016 DO WHILE (abs(dt)>ddt0) -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.F90
r5099 r5101 29 29 USE lmdz_thermcell_ini, ONLY: detr_min, entr_min, detr_q_coef, detr_q_power 30 30 USE lmdz_thermcell_ini, ONLY: mix0, thermals_flag_alim 31 USE lmdz_thermcell_alim, ONLY 32 USE lmdz_thermcell_qsat, ONLY 31 USE lmdz_thermcell_alim, ONLY: thermcell_alim 32 USE lmdz_thermcell_qsat, ONLY: thermcell_qsat 33 33 34 34 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume_6A.F90
r5099 r5101 19 19 USE lmdz_thermcell_ini, ONLY: detr_min, entr_min, detr_q_coef, detr_q_power 20 20 USE lmdz_thermcell_ini, ONLY: mix0, thermals_flag_alim 21 USE lmdz_thermcell_alim, ONLY 22 USE lmdz_thermcell_qsat, ONLY 21 USE lmdz_thermcell_alim, ONLY: thermcell_alim 22 USE lmdz_thermcell_qsat, ONLY: thermcell_qsat 23 23 24 24 … … 724 724 725 725 USE lmdz_thermcell_ini, ONLY: prt_level,fact_thermals_ed_dz,iflag_thermals_ed,RLvCP,RETV,RG 726 USE lmdz_thermcell_qsat, ONLY 726 USE lmdz_thermcell_qsat, ONLY: thermcell_qsat 727 727 IMPLICIT NONE 728 728 -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake.F90
r5099 r5101 27 27 28 28 29 USE lmdz_wake_ini , ONLY 30 USE lmdz_wake_ini , ONLY 31 USE lmdz_wake_ini , ONLY 32 USE lmdz_wake_ini , ONLY 33 USE lmdz_wake_ini , ONLY 34 USE lmdz_wake_ini , ONLY 35 USE lmdz_wake_ini , ONLY 36 USE lmdz_wake_ini , ONLY 37 USE lmdz_wake_ini , ONLY 29 USE lmdz_wake_ini , ONLY: wake_ini 30 USE lmdz_wake_ini , ONLY: prt_level,epsim1,RG,RD 31 USE lmdz_wake_ini , ONLY: stark, wdens_ref, coefgw, alpk, wk_pupper 32 USE lmdz_wake_ini , ONLY: crep_upper, crep_sol, tau_cv, rzero, aa0, flag_wk_check_trgl 33 USE lmdz_wake_ini , ONLY: ok_bug_gfl 34 USE lmdz_wake_ini , ONLY: iflag_wk_act, iflag_wk_check_trgl, iflag_wk_pop_dyn, wdensinit, wdensthreshold 35 USE lmdz_wake_ini , ONLY: sigmad, hwmin, wapecut, cstart, sigmaw_max, dens_rate, epsilon_loc 36 USE lmdz_wake_ini , ONLY: iflag_wk_profile 37 USE lmdz_wake_ini , ONLY: smallestreal,wk_nsub 38 38 39 39 … … 2414 2414 ktop, wk_adv, h_zzz, ptop1, ktop1) 2415 2415 2416 USE lmdz_wake_ini , ONLY 2417 USE lmdz_wake_ini , ONLY 2418 USE lmdz_wake_ini , ONLY 2419 USE lmdz_wake_ini , ONLY 2420 USE lmdz_wake_ini , ONLY 2416 USE lmdz_wake_ini , ONLY: wk_pupper 2417 USE lmdz_wake_ini , ONLY: RG 2418 USE lmdz_wake_ini , ONLY: hwmin 2419 USE lmdz_wake_ini , ONLY: iflag_wk_new_ptop, wk_delta_t_min, wk_frac_int_delta_t 2420 USE lmdz_wake_ini , ONLY: wk_int_delta_t_min 2421 2421 2422 2422 IMPLICIT NONE … … 2817 2817 2818 2818 2819 USE lmdz_wake_ini , ONLY 2820 USE lmdz_wake_ini , ONLY 2821 USE lmdz_wake_ini , ONLY 2822 USE lmdz_wake_ini , ONLY 2823 !! USE lmdz_wake_ini , ONLY 2824 USE lmdz_wake_ini , ONLY 2825 USE lmdz_wake_ini , ONLY 2819 USE lmdz_wake_ini , ONLY: wake_ini 2820 USE lmdz_wake_ini , ONLY: prt_level,RG 2821 USE lmdz_wake_ini , ONLY: stark, wdens_ref 2822 USE lmdz_wake_ini , ONLY: tau_cv, rzero, aa0 2823 !! USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn, wdensmin 2824 USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn 2825 USE lmdz_wake_ini , ONLY: sigmad, cstart, sigmaw_max 2826 2826 2827 2827 IMPLICIT NONE … … 2981 2981 2982 2982 2983 USE lmdz_wake_ini , ONLY 2984 USE lmdz_wake_ini , ONLY 2985 USE lmdz_wake_ini , ONLY 2986 USE lmdz_wake_ini , ONLY 2987 !! USE lmdz_wake_ini , ONLY 2988 USE lmdz_wake_ini , ONLY 2989 USE lmdz_wake_ini , ONLY 2983 USE lmdz_wake_ini , ONLY: wake_ini 2984 USE lmdz_wake_ini , ONLY: prt_level,RG 2985 USE lmdz_wake_ini , ONLY: stark, wdens_ref 2986 USE lmdz_wake_ini , ONLY: tau_cv, rzero, aa0 2987 !! USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn, wdensmin 2988 USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn 2989 USE lmdz_wake_ini , ONLY: sigmad, cstart, sigmaw_max 2990 2990 2991 2991 IMPLICIT NONE … … 3132 3132 3133 3133 3134 USE lmdz_wake_ini , ONLY 3135 USE lmdz_wake_ini , ONLY 3136 USE lmdz_wake_ini , ONLY 3137 USE lmdz_wake_ini , ONLY 3138 !! USE lmdz_wake_ini , ONLY 3139 USE lmdz_wake_ini , ONLY 3140 USE lmdz_wake_ini , ONLY 3141 USE lmdz_wake_ini , ONLY 3134 USE lmdz_wake_ini , ONLY: wake_ini 3135 USE lmdz_wake_ini , ONLY: prt_level,RG 3136 USE lmdz_wake_ini , ONLY: stark, wdens_ref 3137 USE lmdz_wake_ini , ONLY: tau_cv, rzero, aa0 3138 !! USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn, wdensmin 3139 USE lmdz_wake_ini , ONLY: iflag_wk_pop_dyn 3140 USE lmdz_wake_ini , ONLY: sigmad, cstart, sigmaw_max 3141 USE lmdz_wake_ini , ONLY: smallestreal 3142 3142 3143 3143 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake_ini.F90
r5099 r5101 102 102 ! ------------------------------------------------------------------------- 103 103 104 USE ioipsl_getin_p_mod, ONLY 104 USE ioipsl_getin_p_mod, ONLY: getin_p 105 105 real eps 106 106 integer, intent(in) :: prt_lev -
LMDZ6/branches/Amaury_dev/libf/phylmd/lsc_scav.F90
r5099 r5101 10 10 USE mod_phys_lmdz_para 11 11 USE traclmdz_mod 12 USE infotrac_phy,ONLY 12 USE infotrac_phy,ONLY: nbtr 13 13 USE iophy 14 14 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/macv2sp.F90
r5099 r5101 14 14 15 15 USE mo_simple_plumes, ONLY: sp_aop_profile 16 USE phys_cal_mod, ONLY 16 USE phys_cal_mod, ONLY: year_cur, day_cur, year_len 17 17 USE dimphy 18 18 USE aero_mod 19 19 USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer, dNovrN 20 !!USE YOMCST, ONLY 20 !!USE YOMCST, ONLY: RD, RG 21 21 22 22 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/methox.F90
r5087 r5101 59 59 ! ------------------------------------------------------------------ 60 60 61 USE YOEMETH , ONLY 61 USE YOEMETH , ONLY: RALPHA1 ,RALPHA2 ,RQLIM ,& 62 62 RPBOTOX, RPBOTPH ,RPTOPOX ,RPTOPPH ,& 63 63 RALPHA3, RLOGPPH -
LMDZ6/branches/Amaury_dev/libf/phylmd/nuage.F90
r5099 r5101 7 7 USE lmdz_lscp_tools, only: icefrac_lscp 8 8 USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14) 9 USE lmdz_lscp_ini, only: iflag_t_glace9 USE lmdz_lscp_ini, ONLY: iflag_t_glace 10 10 USE phys_local_var_mod, ONLY: ptconv 11 11 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/oasis.F90
r5099 r5101 112 112 113 113 USE IOIPSL 114 USE surface_data, ONLY 115 USE carbon_cycle_mod, ONLY 116 USE wxios, ONLY : wxios_context_init117 USE chemistry_cycle_mod, ONLY 114 USE surface_data, ONLY: version_ocean 115 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 116 USE wxios, ONLY: wxios_context_init 117 USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl 118 118 USE lmdz_xios 119 119 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/ocean_albedo.F90
r5099 r5101 44 44 USE ocean_albedo_para 45 45 USE dimphy 46 USE phys_state_var_mod, ONLY 46 USE phys_state_var_mod, ONLY: chl_con 47 47 48 48 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ocean_cpl_mod.F90
r5099 r5101 23 23 ! Allocate fields for this module and initailize the module mod_cpl 24 24 25 USE dimphy, ONLY 25 USE dimphy, ONLY: klon 26 26 USE cpl_mod 27 27 … … 61 61 ! surface is done and finally it sends some fields to the coupler. 62 62 63 USE dimphy, ONLY 63 USE dimphy, ONLY: klon 64 64 USE calcul_fluxs_mod 65 65 USE indice_sol_mod 66 USE phys_output_var_mod, ONLY 67 USE cpl_mod, ONLY 66 USE phys_output_var_mod, ONLY: sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o 67 USE cpl_mod, ONLY: gath2cpl, cpl_receive_ocean_fields, & 68 68 cpl_send_ocean_fields 69 69 use config_ocean_skin_m, only: activate_ocean_skin … … 252 252 ! some fields to the coupler. 253 253 254 USE dimphy, ONLY 254 USE dimphy, ONLY: klon 255 255 USE cpl_mod 256 256 USE calcul_fluxs_mod 257 257 USE indice_sol_mod 258 USE phys_output_var_mod, ONLY 258 USE phys_output_var_mod, ONLY: sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o 259 259 260 260 INCLUDE "YOMCST.h" -
LMDZ6/branches/Amaury_dev/libf/phylmd/ocean_forced_mod.F90
r5099 r5101 39 39 USE mod_grid_phy_lmdz 40 40 USE indice_sol_mod 41 USE phys_output_var_mod, ONLY 41 USE phys_output_var_mod, ONLY: sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o 42 42 use config_ocean_skin_m, only: activate_ocean_skin 43 43 #ifdef ISO … … 262 262 USE geometry_mod, ONLY: longitude,latitude 263 263 USE calcul_fluxs_mod 264 USE surface_data, ONLY 264 USE surface_data, ONLY: calice, calsno 265 265 USE limit_read_mod 266 USE fonte_neige_mod, ONLY 266 USE fonte_neige_mod, ONLY: fonte_neige 267 267 USE indice_sol_mod 268 USE phys_output_var_mod, ONLY 268 USE phys_output_var_mod, ONLY: sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o 269 269 #ifdef ISO 270 270 USE infotrac_phy, ONLY: niso, ntiso … … 511 511 ! This subroutine specifies the surface temperature to be used in 1D simulations 512 512 513 ! USE dimphy, ONLY 513 ! USE dimphy, ONLY: klon 514 514 515 515 ! INTEGER, INTENT(IN) :: knon ! nomber of points on compressed grid -
LMDZ6/branches/Amaury_dev/libf/phylmd/ocean_slab_mod.F90
r5099 r5101 128 128 !, seaice_rst etc 129 129 130 USE ioipsl_getin_p_mod, ONLY 131 USE mod_phys_lmdz_transfert_para, ONLY 132 USE slab_heat_transp_mod, ONLY 130 USE ioipsl_getin_p_mod, ONLY: getin_p 131 USE mod_phys_lmdz_transfert_para, ONLY: gather 132 USE slab_heat_transp_mod, ONLY: ini_slab_transp 133 133 134 134 ! Input variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/open_climoz_m.F90
r5100 r5101 54 54 nlev = SIZE(press_in_cen) 55 55 CALL nf95_inq_varid(ncID, "time", varID) 56 CALL NF95_GW_VAR(ncid, varid, time_in)56 CALL nf95_gw_var(ncid, varid, time_in) 57 57 ntim = SIZE(time_in) 58 58 -
LMDZ6/branches/Amaury_dev/libf/phylmd/output_physiqex_mod.F90
r4547 r5101 7 7 SUBROUTINE output_physiqex(debut,zjulian,pdtphys,presnivs,paprs,u,v,t,qx,cf,zqr,zqs,zqg,ptke,theta) 8 8 9 USE dimphy, only: klon,klev10 USE iophy, only: histbeg_phy,histwrite_phy11 USE ioipsl, only: histvert,histdef,histend,histsync12 USE mod_phys_lmdz_para, only: jj_nb13 USE ioipsl_getin_p_mod, ONLY 9 USE dimphy, ONLY: klon,klev 10 USE iophy, ONLY: histbeg_phy,histwrite_phy 11 USE ioipsl, ONLY: histvert,histdef,histend,histsync 12 USE mod_phys_lmdz_para, ONLY: jj_nb 13 USE ioipsl_getin_p_mod, ONLY: getin_p 14 14 USE mod_grid_phy_lmdz, ONLY: nbp_lon,nbp_lat 15 USE iophy, ONLY 15 USE iophy, ONLY: init_iophy_new 16 16 USE geometry_mod, ONLY: latitude_deg, longitude_deg 17 USE infotrac_phy, only: nqtot17 USE infotrac_phy, ONLY: nqtot 18 18 19 19 -
LMDZ6/branches/Amaury_dev/libf/phylmd/paramlmdz_phy_mod.F90
r5099 r5101 151 151 152 152 #ifdef CPP_RRTM 153 USE YOESW, ONLY 153 USE YOESW, ONLY: RSUN 154 154 #endif 155 155 -
LMDZ6/branches/Amaury_dev/libf/phylmd/pbl_surface_mod.F90
r5099 r5101 11 11 12 12 USE dimphy 13 USE mod_phys_lmdz_para, ONLY 14 USE mod_grid_phy_lmdz, ONLY 13 USE mod_phys_lmdz_para, ONLY: mpi_size 14 USE mod_grid_phy_lmdz, ONLY: klon_glo 15 15 USE ioipsl 16 USE surface_data, ONLY 17 USE surf_land_mod, ONLY 18 USE surf_landice_mod, ONLY 19 USE surf_ocean_mod, ONLY 20 USE surf_seaice_mod, ONLY 21 USE cpl_mod, ONLY 22 USE climb_hq_mod, ONLY 23 USE climb_qbs_mod, ONLY 24 USE climb_wind_mod, ONLY 25 USE coef_diff_turb_mod, ONLY 26 USE lmdz_call_atke, ONLY 27 USE ioipsl_getin_p_mod, ONLY 16 USE surface_data, ONLY: type_ocean, ok_veget, landice_opt 17 USE surf_land_mod, ONLY: surf_land 18 USE surf_landice_mod, ONLY: surf_landice 19 USE surf_ocean_mod, ONLY: surf_ocean 20 USE surf_seaice_mod, ONLY: surf_seaice 21 USE cpl_mod, ONLY: gath2cpl 22 USE climb_hq_mod, ONLY: climb_hq_down, climb_hq_up 23 USE climb_qbs_mod, ONLY: climb_qbs_down, climb_qbs_up 24 USE climb_wind_mod, ONLY: climb_wind_down, climb_wind_up 25 USE coef_diff_turb_mod, ONLY: coef_diff_turb 26 USE lmdz_call_atke, ONLY: call_atke 27 USE ioipsl_getin_p_mod, ONLY: getin_p 28 28 USE cdrag_mod 29 29 USE stdlevvar_mod 30 USE wx_pbl_var_mod, ONLY 30 USE wx_pbl_var_mod, ONLY: wx_pbl_init, wx_pbl_final, & 31 31 wx_pbl_prelim_0, wx_pbl_prelim_beta 32 USE wx_pbl_mod, ONLY 32 USE wx_pbl_mod, ONLY: wx_pbl0_merge, wx_pbl_split, wx_pbl_dts_merge, & 33 33 wx_pbl_check, wx_pbl_dts_check, wx_evappot 34 34 use config_ocean_skin_m, only: activate_ocean_skin … … 84 84 USE indice_sol_mod 85 85 USE print_control_mod, ONLY: lunout 86 USE ioipsl_getin_p_mod, ONLY 86 USE ioipsl_getin_p_mod, ONLY: getin_p 87 87 IMPLICIT NONE 88 88 … … 391 391 ! treedrg--output-R- tree drag (m) 392 392 393 USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, carbon_cycle_tr, level_coupling_esm394 USE carbon_cycle_mod, ONLY 393 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, carbon_cycle_tr, level_coupling_esm 394 USE carbon_cycle_mod, ONLY: co2_send, nbcf_out, fields_out, yfields_out, cfname_out 395 395 use hbtm_mod, only: hbtm 396 396 USE indice_sol_mod 397 USE time_phylmdz_mod, ONLY 398 USE mod_grid_phy_lmdz, ONLY 399 USE print_control_mod, ONLY 397 USE time_phylmdz_mod, ONLY: day_ini,annee_ref,itau_phy 398 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid1dto2d_glo 399 USE print_control_mod, ONLY: prt_level,lunout 400 400 #ifdef ISO 401 401 USE isotopes_mod, ONLY: Rdefault,iso_eau … … 407 407 #endif 408 408 #endif 409 USE ioipsl_getin_p_mod, ONLY 409 USE ioipsl_getin_p_mod, ONLY: getin_p 410 410 use phys_state_var_mod, only: ds_ns, dt_ns, delta_sst, delta_sal, dter, & 411 411 dser, dt_ds, zsig, zmea 412 412 use phys_output_var_mod, only: tkt, tks, taur, sss 413 use lmdz_blowing_snow_ini, only: zeta_bs413 use lmdz_blowing_snow_ini, ONLY: zeta_bs 414 414 USE wxios, ONLY: missing_val_xios => missing_val, using_xios 415 415 USE netcdf, only: missing_val_netcdf => nf90_fill_real -
LMDZ6/branches/Amaury_dev/libf/phylmd/phyaqua_mod.F90
r5099 r5101 23 23 24 24 USE dimphy, ONLY: klon 25 USE geometry_mod, ONLY 25 USE geometry_mod, ONLY: latitude 26 26 USE surface_data, ONLY: type_ocean, ok_veget 27 27 USE pbl_surface_mod, ONLY: pbl_surface_init … … 36 36 USE mod_phys_lmdz_transfert_para, ONLY: bcast 37 37 USE mod_grid_phy_lmdz 38 USE ioipsl_getin_p_mod, ONLY 38 USE ioipsl_getin_p_mod, ONLY: getin_p 39 39 USE phys_cal_mod , ONLY: calend, year_len_phy => year_len 40 40 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90
r5099 r5101 11 11 12 12 USE dimphy, only: klon, zmasq, klev 13 USE iophy, ONLY 14 USE ocean_cpl_mod, ONLY 15 USE fonte_neige_mod, ONLY 16 USE pbl_surface_mod, ONLY 17 USE surface_data, ONLY 18 USE phyetat0_get_mod, ONLY 19 USE phys_state_var_mod, ONLY 13 USE iophy, ONLY: init_iophy_new 14 USE ocean_cpl_mod, ONLY: ocean_cpl_init 15 USE fonte_neige_mod, ONLY: fonte_neige_init 16 USE pbl_surface_mod, ONLY: pbl_surface_init 17 USE surface_data, ONLY: type_ocean, version_ocean 18 USE phyetat0_get_mod, ONLY: phyetat0_get, phyetat0_srf 19 USE phys_state_var_mod, ONLY: ancien_ok, clwcon, detr_therm, phys_tstep, & 20 20 qsol, fevap, z0m, z0h, agesno, & 21 21 du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, & -
LMDZ6/branches/Amaury_dev/libf/phylmd/phyredem.F90
r5099 r5101 10 10 !------------------------------------------------------------------------------- 11 11 USE dimphy, ONLY: klon, klev 12 USE fonte_neige_mod, ONLY 13 USE pbl_surface_mod, ONLY 12 USE fonte_neige_mod, ONLY: fonte_neige_final 13 USE pbl_surface_mod, ONLY: pbl_surface_final 14 14 USE phys_state_var_mod, ONLY: radpas, zmasq, pctsrf, & 15 15 ftsol, beta_aridity, delta_tsurf, falb_dir, & … … 32 32 delta_sst, ratqs_inter_, dter, dser, dt_ds 33 33 34 USE geometry_mod, ONLY 34 USE geometry_mod, ONLY: longitude_deg, latitude_deg 35 35 USE iostart, ONLY: open_restartphy, close_restartphy, enddef_restartphy, put_field, put_var 36 USE traclmdz_mod, ONLY 36 USE traclmdz_mod, ONLY: traclmdz_to_restart 37 37 USE infotrac_phy, ONLY: type_trac, nqtot, tracers, nbtr 38 USE carbon_cycle_mod, ONLY 38 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, co2_send, carbon_cycle_rad, RCO2_glo 39 39 USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic, epsfra 40 40 USE surface_data, ONLY: type_ocean, version_ocean 41 USE ocean_slab_mod, ONLY 41 USE ocean_slab_mod, ONLY: nslay, tslab, seaice, tice, fsic 42 42 USE time_phylmdz_mod, ONLY: annee_ref, day_end, itau_phy, pdtphys 43 43 use config_ocean_skin_m, only: activate_ocean_skin -
LMDZ6/branches/Amaury_dev/libf/phylmd/phys_local_var_mod.F90
r5099 r5101 3 3 4 4 MODULE phys_local_var_mod 5 USE lmdz_cppkeys_wrapper, ONLY 5 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER 6 6 ! Variables locales pour effectuer les appels en serie 7 7 !====================================================================== … … 815 815 SUBROUTINE phys_local_var_init 816 816 USE dimphy 817 USE infotrac_phy, ONLY 818 #ifdef ISO 819 USE infotrac_phy, ONLY 817 USE infotrac_phy, ONLY: nbtr, nqtot 818 #ifdef ISO 819 USE infotrac_phy, ONLY: ntraciso=>ntiso,niso 820 820 #endif 821 821 USE aero_mod … … 823 823 USE phys_output_var_mod 824 824 USE phys_state_var_mod 825 USE infotrac_phy, ONLY 825 USE infotrac_phy, ONLY: nbtr_bin 826 826 827 827 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/phys_output_mod.F90
r5099 r5101 4 4 USE indice_sol_mod 5 5 USE phys_output_var_mod 6 USE phys_output_write_mod, ONLY 6 USE phys_output_write_mod, ONLY: phys_output_write 7 7 REAL, DIMENSION(nfiles),SAVE :: ecrit_files 8 8 … … 37 37 USE strings_mod, ONLY: maxlen 38 38 USE ioipsl 39 USE phys_cal_mod, only: hour, calend39 USE phys_cal_mod, ONLY: hour, calend 40 40 USE mod_phys_lmdz_para 41 41 !Martin 42 USE surface_data, ONLY 42 USE surface_data, ONLY: landice_opt 43 43 USE phys_output_ctrlout_mod 44 44 USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat -
LMDZ6/branches/Amaury_dev/libf/phylmd/phys_output_write_mod.F90
r5099 r5101 4 4 MODULE phys_output_write_mod 5 5 6 USE phytrac_mod, ONLY 6 USE phytrac_mod, ONLY: d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, & 7 7 d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls, & 8 8 d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry … … 28 28 USE strings_mod, ONLY: maxlen 29 29 USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy 30 USE mod_grid_phy_lmdz, ONLY 30 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 31 31 USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy 32 USE vertical_layers_mod, ONLY 32 USE vertical_layers_mod, ONLY: ap, bp, aps, bps 33 33 USE phystokenc_mod, ONLY: offline 34 34 USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, & … … 458 458 USE lmdz_xios 459 459 USE wxios, ONLY: wxios_closedef, missing_val_xios=>missing_val, wxios_set_context 460 USE phys_cal_mod, ONLY 460 USE phys_cal_mod, ONLY: mth_len 461 461 462 462 #ifdef CPP_RRTM 463 USE YOESW, ONLY 463 USE YOESW, ONLY: RSUN 464 464 #endif 465 465 USE tracinca_mod, ONLY: config_inca -
LMDZ6/branches/Amaury_dev/libf/phylmd/phys_state_var_mod.F90
r5099 r5101 535 535 USE dimphy 536 536 USE aero_mod 537 USE infotrac_phy, ONLY 537 USE infotrac_phy, ONLY: nbtr 538 538 #ifdef ISO 539 USE infotrac_phy, ONLY 539 USE infotrac_phy, ONLY: ntraciso=>ntiso,niso 540 540 #endif 541 541 USE indice_sol_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/physiq_mod.F90
r5099 r5101 19 19 ! PLEASE try to follow this rule 20 20 21 USE ACAMA_GWD_rando_m, only: ACAMA_GWD_rando21 USE ACAMA_GWD_rando_m, ONLY: ACAMA_GWD_rando 22 22 USE aero_mod 23 USE add_phys_tend_mod, only: add_pbl_tend, add_phys_tend, diag_phys_tend, prt_enerbil, &23 USE add_phys_tend_mod, ONLY: add_pbl_tend, add_phys_tend, diag_phys_tend, prt_enerbil, & 24 24 fl_ebil, fl_cor_ebil 25 USE assert_m, only: assert25 USE assert_m, ONLY: assert 26 26 USE change_srf_frac_mod 27 USE conf_phys_m, only: conf_phys28 USE carbon_cycle_mod, ONLY 27 USE conf_phys_m, ONLY: conf_phys 28 USE carbon_cycle_mod, ONLY: infocfields_init, RCO2_glo, carbon_cycle_rad 29 29 USE CFMIP_point_locations ! IM stations CFMIP 30 30 USE cmp_seri_mod 31 31 USE dimphy 32 32 USE etat0_limit_unstruct_mod 33 USE FLOTT_GWD_rando_m, only: FLOTT_GWD_rando34 USE fonte_neige_mod, ONLY 35 USE geometry_mod, ONLY 36 USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, &33 USE FLOTT_GWD_rando_m, ONLY: FLOTT_GWD_rando 34 USE fonte_neige_mod, ONLY: fonte_neige_get_vars 35 USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg 36 USE ioipsl, ONLY: histbeg, histvert, histdef, histend, histsync, & 37 37 histwrite, ju2ymds, ymds2ju, getin 38 USE ioipsl_getin_p_mod, ONLY 38 USE ioipsl_getin_p_mod, ONLY: getin_p 39 39 USE indice_sol_mod 40 USE infotrac_phy, ONLY 41 USE readTracFiles_mod, ONLY 42 USE strings_mod, ONLY 40 USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac 41 USE readTracFiles_mod, ONLY: addPhase 42 USE strings_mod, ONLY: strIdx 43 43 USE iophy 44 USE limit_read_mod, ONLY 45 USE mod_grid_phy_lmdz, ONLY 46 USE mod_phys_lmdz_mpi_data, only: is_mpi_root44 USE limit_read_mod, ONLY: init_limit_read 45 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev, klon_glo, grid1dTo2d_glo, grid_type, unstructured 46 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 47 47 USE mod_phys_lmdz_para 48 USE netcdf95, only: nf95_close49 USE netcdf, only: nf90_fill_real ! IM for NMC files50 USE open_climoz_m, only: open_climoz ! ozone climatology from a file51 USE ozonecm_m, only: ozonecm ! ozone of J.-F. Royer52 USE pbl_surface_mod, ONLY 53 USE phyaqua_mod, only: zenang_an54 USE phyetat0_mod, only: phyetat055 USE phystokenc_mod, ONLY 56 USE phys_cal_mod, only: year_len, mth_len, days_elapsed, jh_1jan, &48 USE netcdf95, ONLY: nf95_close 49 USE netcdf, ONLY: nf90_fill_real ! IM for NMC files 50 USE open_climoz_m, ONLY: open_climoz ! ozone climatology from a file 51 USE ozonecm_m, ONLY: ozonecm ! ozone of J.-F. Royer 52 USE pbl_surface_mod, ONLY: pbl_surface 53 USE phyaqua_mod, ONLY: zenang_an 54 USE phyetat0_mod, ONLY: phyetat0 55 USE phystokenc_mod, ONLY: offline, phystokenc 56 USE phys_cal_mod, ONLY: year_len, mth_len, days_elapsed, jh_1jan, & 57 57 year_cur, mth_cur, jD_cur, jH_cur, jD_ref, day_cur, hour, calend 58 !! USE phys_local_var_mod, ONLY 58 !! USE phys_local_var_mod, ONLY: a long list of variables 59 59 !! ==> see below, after "CPP Keys" section 60 60 USE phys_state_var_mod ! Variables sauvegardees de la physique 61 61 USE phys_output_mod 62 62 USE phys_output_ctrlout_mod 63 USE print_control_mod, ONLY 63 USE print_control_mod, ONLY: mydebug => debug, lunout, prt_level, & 64 64 alert_first_call, call_alert, prt_alerte 65 USE readaerosol_mod, ONLY 66 USE readaerosolstrato_m, ONLY 67 USE radlwsw_m, only: radlwsw68 USE regr_horiz_time_climoz_m, ONLY 69 USE regr_pr_time_av_m, only: regr_pr_time_av70 USE surface_data, ONLY 71 USE time_phylmdz_mod, only: current_time, itau_phy, pdtphys, raz_date, update_time72 USE tracinca_mod, ONLY 73 USE tropopause_m, ONLY 74 USE ice_sursat_mod, ONLY 65 USE readaerosol_mod, ONLY: init_aero_fromfile 66 USE readaerosolstrato_m, ONLY: init_readaerosolstrato 67 USE radlwsw_m, ONLY: radlwsw 68 USE regr_horiz_time_climoz_m, ONLY: regr_horiz_time_climoz 69 USE regr_pr_time_av_m, ONLY: regr_pr_time_av 70 USE surface_data, ONLY: type_ocean, ok_veget 71 USE time_phylmdz_mod, ONLY: current_time, itau_phy, pdtphys, raz_date, update_time 72 USE tracinca_mod, ONLY: config_inca 73 USE tropopause_m, ONLY: dyn_tropopause 74 USE ice_sursat_mod, ONLY: flight_init, airplane 75 75 USE vampir 76 76 USE write_field_phy 77 USE wxios, ONLY 78 USE lmdz_lscp, ONLY 79 USE lmdz_call_cloud_optics_prop, ONLY 80 USE lmdz_lscp_old, ONLY 81 USE lmdz_call_blowing_snow, ONLY 82 USE lmdz_wake_ini, ONLY 83 USE yamada_ini_mod, ONLY 84 USE lmdz_atke_turbulence_ini, ONLY 85 USE lmdz_thermcell_ini, ONLY 86 USE lmdz_thermcell_dtke, ONLY 87 USE lmdz_blowing_snow_ini, ONLY 88 USE lmdz_lscp_ini, ONLY 89 USE lmdz_ratqs_main, ONLY 90 USE lmdz_ratqs_ini, ONLY 91 USE lmdz_cloud_optics_prop_ini, ONLY 92 USE phys_output_var_mod, ONLY 93 USE phys_output_var_mod, ONLY 94 95 USE geometry_mod, ONLY 96 USE time_phylmdz_mod, ONLY 97 USE infotrac_phy, ONLY 98 99 USE time_phylmdz_mod, ONLY 100 USE vertical_layers_mod, ONLY 101 102 USE lmdz_xios, ONLY 77 USE wxios, ONLY: g_ctx, wxios_set_context 78 USE lmdz_lscp, ONLY: lscp 79 USE lmdz_call_cloud_optics_prop, ONLY: call_cloud_optics_prop 80 USE lmdz_lscp_old, ONLY: fisrtilp 81 USE lmdz_call_blowing_snow, ONLY: call_blowing_snow_sublim_sedim 82 USE lmdz_wake_ini, ONLY: wake_ini 83 USE yamada_ini_mod, ONLY: yamada_ini 84 USE lmdz_atke_turbulence_ini, ONLY: atke_ini 85 USE lmdz_thermcell_ini, ONLY: thermcell_ini, iflag_thermals_tenv 86 USE lmdz_thermcell_dtke, ONLY: thermcell_dtke 87 USE lmdz_blowing_snow_ini, ONLY: blowing_snow_ini, qbst_bs 88 USE lmdz_lscp_ini, ONLY: lscp_ini 89 USE lmdz_ratqs_main, ONLY: ratqs_main 90 USE lmdz_ratqs_ini, ONLY: ratqs_ini 91 USE lmdz_cloud_optics_prop_ini, ONLY: cloud_optics_prop_ini 92 USE phys_output_var_mod, ONLY: cloudth_sth, cloudth_senv, cloudth_sigmath, cloudth_sigmaenv 93 USE phys_output_var_mod, ONLY: cloud_cover_sw, cloud_cover_sw_s2 94 95 USE geometry_mod, ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo 96 USE time_phylmdz_mod, ONLY: ndays 97 USE infotrac_phy, ONLY: nqCO2 98 99 USE time_phylmdz_mod, ONLY: annee_ref, day_ini, day_ref, start_time 100 USE vertical_layers_mod, ONLY: aps, bps, ap, bp 101 102 USE lmdz_xios, ONLY: xios_update_calendar, xios_context_finalize, xios_get_field_attr, & 103 103 xios_field_is_active, xios_context, xios_set_current_context 104 USE wxios, ONLY 104 USE wxios, ONLY: missing_val, using_xios 105 105 USE paramLMDZ_phy_mod 106 106 107 USE phytracr_spl_mod, ONLY 107 USE phytracr_spl_mod, ONLY: phytracr_spl, phytracr_spl_out_init 108 108 USE phys_output_write_spl_mod 109 USE phytrac_mod, ONLY 109 USE phytrac_mod, ONLY: phytrac_init, phytrac 110 110 USE phys_output_write_mod 111 111 112 USE lmdz_cppkeys_wrapper, ONLY 112 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DUST 113 113 114 114 !!!!!!!!!!!!!!!!!! "USE" section for CPP keys !!!!!!!!!!!!!!!!!!!!!!!! … … 123 123 124 124 #ifdef CPP_RRTM 125 USE YOERAD, ONLY 125 USE YOERAD, ONLY: NRADLP 126 126 #endif 127 127 128 USE phys_local_var_mod, ONLY 128 USE phys_local_var_mod, ONLY: d_q_emiss 129 129 USE strataer_local_var_mod 130 USE strataer_nuc_mod, ONLY 131 USE strataer_emiss_mod, ONLY 132 USE lmdz_cppkeys_wrapper, ONLY 130 USE strataer_nuc_mod, ONLY: strataer_nuc_init 131 USE strataer_emiss_mod, ONLY: strataer_emiss_init 132 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER 133 133 134 134 135 135 !!!!!!!!!!!!!!!!!! END "USE" for CPP keys !!!!!!!!!!!!!!!!!!!!!! 136 136 137 USE physiqex_mod, ONLY 138 USE phys_local_var_mod, ONLY 137 USE physiqex_mod, ONLY: physiqex 138 USE phys_local_var_mod, ONLY: phys_local_var_init, phys_local_var_end, & 139 139 ! [Variables internes non sauvegardees de la physique] 140 140 ! Variables locales pour effectuer les appels en serie … … 345 345 zxfluxt, zxfluxq 346 346 347 USE phys_local_var_mod, ONLY 348 USE phys_output_var_mod, ONLY 347 USE phys_local_var_mod, ONLY: zfice, dNovrN, ptconv 348 USE phys_output_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, & 349 349 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra 350 USE output_physiqex_mod, ONLY 351 USE lmdz_cppkeys_wrapper, ONLY 350 USE output_physiqex_mod, ONLY: output_physiqex 351 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 352 352 353 353 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/physiqex_mod.F90
r5099 r5101 13 13 d_u, d_v, d_t, d_qx, d_ps) 14 14 15 USE dimphy, only: klon,klev16 USE infotrac_phy, only: nqtot17 USE geometry_mod, only: latitude18 ! USE comcstphy, only: rg19 USE ioipsl, only: ymds2ju20 USE phys_state_var_mod, only: phys_state_var_init15 USE dimphy, ONLY: klon,klev 16 USE infotrac_phy, ONLY: nqtot 17 USE geometry_mod, ONLY: latitude 18 ! USE comcstphy, ONLY: rg 19 USE ioipsl, ONLY: ymds2ju 20 USE phys_state_var_mod, ONLY: phys_state_var_init 21 21 USE phyetat0_mod, only: phyetat0 22 22 USE output_physiqex_mod, ONLY: output_physiqex -
LMDZ6/branches/Amaury_dev/libf/phylmd/phystokenc_mod.F90
r5099 r5101 34 34 USE ioipsl 35 35 USE dimphy 36 USE infotrac_phy, ONLY 36 USE infotrac_phy, ONLY: nqtot 37 37 USE iophy 38 38 USE indice_sol_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/phytrac_mod.F90
r5099 r5101 55 55 56 56 USE dimphy 57 USE infotrac_phy, ONLY 58 USE tracco2i_mod, ONLY 57 USE infotrac_phy, ONLY: nbtr, type_trac 58 USE tracco2i_mod, ONLY: tracco2i_init 59 59 60 60 IMPLICIT NONE … … 120 120 121 121 USE ioipsl 122 USE phys_cal_mod, only: hour122 USE phys_cal_mod, ONLY: hour 123 123 USE dimphy 124 USE infotrac_phy, ONLY 125 USE strings_mod, ONLY 124 USE infotrac_phy, ONLY: nbtr, nqCO2, type_trac, conv_flg, pbl_flg 125 USE strings_mod, ONLY: int2str 126 126 USE mod_grid_phy_lmdz 127 127 USE mod_phys_lmdz_para … … 131 131 USE tracreprobus_mod 132 132 USE indice_sol_mod 133 USE mod_phys_lmdz_mpi_data, ONLY 134 USE print_control_mod, ONLY 135 USE aero_mod, ONLY 136 USE lmdz_thermcell_dq, ONLY 133 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 134 USE print_control_mod, ONLY: lunout 135 USE aero_mod, ONLY: naero_grp 136 USE lmdz_thermcell_dq, ONLY: thermcell_dq 137 137 138 138 USE tracco2i_mod 139 139 140 140 USE traccoag_mod 141 USE phys_local_var_mod, ONLY 142 USE phys_local_var_mod, ONLY 143 USE phys_local_var_mod, ONLY 144 USE phys_local_var_mod, ONLY 145 USE phys_local_var_mod, ONLY 146 USE infotrac_phy, ONLY 147 USE strataer_nuc_mod, ONLY 141 USE phys_local_var_mod, ONLY: mdw 142 USE phys_local_var_mod, ONLY: budg_dep_dry_ocs, budg_dep_wet_ocs 143 USE phys_local_var_mod, ONLY: budg_dep_dry_so2, budg_dep_wet_so2 144 USE phys_local_var_mod, ONLY: budg_dep_dry_h2so4, budg_dep_wet_h2so4 145 USE phys_local_var_mod, ONLY: budg_dep_dry_part, budg_dep_wet_part 146 USE infotrac_phy, ONLY: nbtr_sulgas, id_OCS_strat, id_SO2_strat, id_H2SO4_strat 147 USE strataer_nuc_mod, ONLY: tracstrataer_init 148 148 USE aerophys 149 USE lmdz_cppkeys_wrapper, ONLY 149 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER 150 150 USE lmdz_yomcst 151 151 -
LMDZ6/branches/Amaury_dev/libf/phylmd/pppmer.F90
r2386 r5101 53 53 54 54 ! USE PARKIND1 55 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/parkind1.F90.php#parkind1> ,ONLY 55 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/parkind1.F90.php#parkind1> ,ONLY: JPIM ,JPRB 56 56 ! USE YOMHOOK 57 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/yomhook.F90.php#yomhook> ,ONLY 57 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/module/yomhook.F90.php#yomhook> ,ONLY: LHOOK, DR_HOOK 58 58 59 !USE YOMCST, ONLY 60 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomcst.F90.php#yomcst> , ONLY 59 !USE YOMCST, ONLY: RG, RD 60 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomcst.F90.php#yomcst> , ONLY: RG 61 61 62 62 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/xrd/support/rg.F.php#rg> ,RD 63 63 ! USE YOMSTA 64 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomsta.F90.php#yomsta> , ONLY 64 !<http://intra.cnrm.meteo.fr/eac/ARPCLI5.2/doci/code/arpcli5.2/arp/module/yomsta.F90.php#yomsta> , ONLY: RDTDZ1 65 65 66 66 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/qsat_seawater2.F90
r5099 r5101 49 49 ! ------------ 50 50 51 USE MODD_CSTS, ONLY 51 USE MODD_CSTS, ONLY: XRD, XRV 52 52 USE dimphy 53 53 USE indice_sol_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/radio_decay.F90
r5099 r5101 7 7 8 8 USE dimphy 9 USE infotrac_phy, ONLY 10 USE traclmdz_mod, ONLY 9 USE infotrac_phy, ONLY: nbtr 10 USE traclmdz_mod, ONLY: id_rn, id_pb 11 11 IMPLICIT NONE 12 12 !----------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/radlwsw_m.F90
r5099 r5101 49 49 ! Modules necessaires 50 50 USE DIMPHY 51 USE assert_m, ONLY 52 USE infotrac_phy, ONLY 51 USE assert_m, ONLY: assert 52 USE infotrac_phy, ONLY: type_trac 53 53 USE write_field_phy 54 54 55 55 #ifdef REPROBUS 56 USE CHEM_REP, ONLY 56 USE CHEM_REP, ONLY: solaireTIME, ok_SUNTIME, ndimozon 57 57 #endif 58 58 … … 60 60 ! modules necessaires au rayonnement 61 61 ! ----------------------------------------- 62 USE YOERAD , ONLY 62 USE YOERAD , ONLY: NLW, LRRTM ,LCCNL ,LCCNO ,& 63 63 NRADIP , NRADLP , NICEOPT, NLIQOPT ,RCCNLND , RCCNSEA 64 USE YOELW , ONLY 65 USE YOESW , ONLY : RYFWCA ,RYFWCB ,RYFWCC ,RYFWCD,&64 USE YOELW , ONLY: NSIL ,NTRA ,NUA ,TSTAND ,XP 65 USE YOESW , ONLY: RYFWCA ,RYFWCB ,RYFWCC ,RYFWCD,& 66 66 RYFWCE ,RYFWCF ,REBCUA ,REBCUB ,REBCUC,& 67 67 REBCUD ,REBCUE ,REBCUF ,REBCUI ,REBCUJ,& … … 71 71 RFLDD1 ,RFLDD2 ,RFLDD3 ,RFUETA ,RASWCA,& 72 72 RASWCB ,RASWCC ,RASWCD ,RASWCE ,RASWCF 73 USE YOERDU , ONLY 74 USE YOERRTWN , ONLY : DELWAVE ,TOTPLNK75 USE YOMPHY3 , ONLY 73 USE YOERDU , ONLY: NUAER ,NTRAER ,REPLOG ,REPSC ,REPSCW ,DIFF 74 USE YOERRTWN , ONLY: DELWAVE ,TOTPLNK 75 USE YOMPHY3 , ONLY: RII0 76 76 #endif 77 77 USE aero_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/read_map2D.F90
r5100 r5101 33 33 IF (ierr /= nf90_noerr) CALL write_err_mess('Problem in opening file') 34 34 35 ierr = NF90_INQ_VARID(nid, trim(varname), nvarid)35 ierr = nf90_inq_varid(nid, trim(varname), nvarid) 36 36 IF (ierr /= nf90_noerr) CALL write_err_mess('The variable is absent in file') 37 37 -
LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_interp.F90
r5099 r5101 14 14 15 15 USE ioipsl 16 USE dimphy, ONLY 17 USE mod_phys_lmdz_para, ONLY : mpi_rank16 USE dimphy, ONLY: klev,klon 17 USE mod_phys_lmdz_para, ONLY: mpi_rank 18 18 USE readaerosol_mod 19 USE aero_mod, ONLY 19 USE aero_mod, ONLY: naero_spc, name_aero 20 20 USE write_field_phy 21 21 USE phys_cal_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_mod.F90
r5100 r5101 271 271 grid2Dto1D_glo, grid_type, unstructured 272 272 USE mod_phys_lmdz_para 273 USE iophy, ONLY 273 USE iophy, ONLY: io_lon, io_lat 274 274 USE print_control_mod, ONLY: lunout 275 275 USE lmdz_xios -
LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosolstrato.F90
r5099 r5101 5 5 use netcdf, only: nf90_get_var, nf90_noerr, nf90_nowrite 6 6 7 USE phys_cal_mod, ONLY 7 USE phys_cal_mod, ONLY: mth_cur 8 8 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo, & 9 9 grid2dto1d_glo, grid_type, unstructured 10 USE mod_phys_lmdz_mpi_data, ONLY 11 USE mod_phys_lmdz_omp_data, ONLY 10 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root 11 USE mod_phys_lmdz_omp_data, ONLY: is_omp_root 12 12 USE mod_phys_lmdz_para 13 13 USE phys_state_var_mod -
LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosolstrato_m.F90
r5088 r5101 29 29 USE mod_phys_lmdz_para 30 30 USE lmdz_xios 31 ! USE YOERAD, ONLY 31 ! USE YOERAD, ONLY: NLW 32 32 IMPLICIT NONE 33 33 REAL, allocatable:: latitude(:) … … 72 72 USE mod_phys_lmdz_para 73 73 USE lmdz_xios 74 ! USE YOERAD, ONLY 74 ! USE YOERAD, ONLY: NLW 75 75 IMPLICIT NONE 76 76 REAL, allocatable:: latitude(:) -
LMDZ6/branches/Amaury_dev/libf/phylmd/reevap.F90
r5099 r5101 3 3 4 4 ! flag to include modifications to ensure energy conservation (if flag >0) 5 USE add_phys_tend_mod, only : fl_cor_ebil5 USE add_phys_tend_mod, ONLY: fl_cor_ebil 6 6 7 7 IMPLICIT none -
LMDZ6/branches/Amaury_dev/libf/phylmd/regr_horiz_time_climoz_m.F90
r5100 r5101 1 1 MODULE regr_horiz_time_climoz_m 2 2 3 USE interpolation, ONLY 4 USE mod_grid_phy_lmdz, ONLY 5 USE nrtype, ONLY 6 USE netcdf, ONLY 3 USE interpolation, ONLY: locate 4 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, unstructured 5 USE nrtype, ONLY: pi 6 USE netcdf, ONLY: nf90_clobber, nf90_float, nf90_open, & 7 7 nf90_nowrite, nf90_noerr, nf90_get_att, nf90_global 8 USE netcdf95, ONLY 9 nf95_def_var, nf95_inq_varid, NF95_INQUIRE_VARIABLE, &10 NF95_OPEN, NF95_CREATE, NF95_GET_ATT, NF95_GW_VAR, nf95_get_var, &11 NF95_CLOSE, NF95_ENDDEF, NF95_PUT_ATT, NF95_PUT_VAR, NF95_COPY_ATT12 USE print_control_mod, ONLY 8 USE netcdf95, ONLY: nf95_def_dim, nf95_inq_dimid, nf95_inquire_dimension, & 9 nf95_def_var, nf95_inq_varid, nf95_inquire_variable, & 10 nf95_open, nf95_create, nf95_get_att, nf95_gw_var, nf95_get_var, & 11 nf95_close, nf95_enddef, nf95_put_att, nf95_put_var, nf95_copy_att 12 USE print_control_mod, ONLY: lunout 13 13 USE dimphy 14 14 IMPLICIT NONE … … 57 57 ! Attribute "missing_value" or "_FillValue" must be present in input file. 58 58 !------------------------------------------------------------------------------- 59 USE assert_m, ONLY 60 USE cal_tools_m, ONLY 59 USE assert_m, ONLY: assert 60 USE cal_tools_m, ONLY: year_len, mid_month 61 61 !! USE control_mod, ONLY: anneeref 62 USE time_phylmdz_mod, ONLY 63 USE ioipsl, ONLY 64 USE regr_conserv_m, ONLY 65 USE regr_lint_m, ONLY 66 USE regular_lonlat_mod, ONLY 67 USE slopes_m, ONLY 62 USE time_phylmdz_mod, ONLY: annee_ref 63 USE ioipsl, ONLY: ioget_year_len, ioget_calendar 64 USE regr_conserv_m, ONLY: regr_conserv 65 USE regr_lint_m, ONLY: regr_lint 66 USE regular_lonlat_mod, ONLY: boundslon_reg, boundslat_reg, south, west, east 67 USE slopes_m, ONLY: slopes 68 68 USE lmdz_xios 69 USE mod_phys_lmdz_para, ONLY 70 USE geometry_mod, ONLY 71 USE mod_grid_phy_lmdz, ONLY 69 USE mod_phys_lmdz_para, ONLY: is_mpi_root, is_master, is_omp_master, gather, gather_mpi, bcast_mpi, klon_mpi 70 USE geometry_mod, ONLY: latitude_deg, ind_cell_glo 71 USE mod_grid_phy_lmdz, ONLY: klon_glo 72 72 73 73 !------------------------------------------------------------------------------- … … 164 164 CALL assert(read_climoz == 1 .OR. read_climoz == 2, "regr_lat_time_climoz") 165 165 166 CALL NF95_OPEN("climoz.nc", nf90_nowrite, fID_in)166 CALL nf95_open("climoz.nc", nf90_nowrite, fID_in) 167 167 lprev = nf90_open("climoz_m.nc", nf90_nowrite, fID_in_m)==nf90_noerr 168 168 lnext = nf90_open("climoz_p.nc", nf90_nowrite, fID_in_p)==nf90_noerr … … 171 171 ! Few inversions because "regr_conserv" and gcm need ascending vectors. 172 172 CALL nf95_inq_varid(fID_in, vars_in(1), varid) 173 CALL NF95_INQUIRE_VARIABLE(fID_in, varid, dimids = dIDs, ndims = ndims)173 CALL nf95_inquire_variable(fID_in, varid, dimids = dIDs, ndims = ndims) 174 174 l3D = ndims==4; l2D = ndims==3 175 175 IF(l3D) WRITE(lunout, *)"Input files contain full 3D ozone fields." … … 181 181 SELECT CASE(ii) 182 182 CASE(1) !--- LONGITUDE 183 CALL NF95_GW_VAR(fID_in, varid, lon_in)183 CALL nf95_gw_var(fID_in, varid, lon_in) 184 184 ldec_lon = lon_in(1)>lon_in(dln); IF(ldec_lon) lon_in = lon_in(dln:1:-1) 185 185 nlon_in = dln; lon_in = lon_in * deg2rad 186 186 CASE(2) !--- LATITUDE 187 CALL NF95_GW_VAR(fID_in, varid, lat_in)187 CALL nf95_gw_var(fID_in, varid, lat_in) 188 188 ldec_lat = lat_in(1)>lat_in(dln); IF(ldec_lat) lat_in = lat_in(dln:1:-1) 189 189 nlat_in = dln; lat_in = lat_in * deg2rad 190 190 CASE(3) !--- PRESSURE LEVELS 191 CALL NF95_GW_VAR(fID_in, varid, lev_in)191 CALL nf95_gw_var(fID_in, varid, lev_in) 192 192 ldec_lev = lev_in(1)>lev_in(dln); IF(ldec_lev) lev_in = lev_in(dln:1:-1) 193 193 nlev_in = dln 194 CALL NF95_GET_ATT(fID_in, varid, "units", press_unit)194 CALL nf95_get_att(fID_in, varid, "units", press_unit) 195 195 k = LEN_TRIM(press_unit) 196 196 DO WHILE(ICHAR(press_unit(k:k))==0) … … 329 329 DO iv = 1, read_climoz 330 330 CALL nf95_inq_varid(fID_in, vars_in(1), vID_in(iv)) 331 IF(l3D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in3(1:nlon_in, :, :, 1:12, iv))332 IF(l2D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in2(:, :, 1:12, iv))331 IF(l3D) call nf95_get_var(fID_in, vID_in(iv), o3_in3(1:nlon_in, :, :, 1:12, iv)) 332 IF(l2D) call nf95_get_var(fID_in, vID_in(iv), o3_in2(:, :, 1:12, iv)) 333 333 IF(lprev) THEN; sta(ndims) = 12 334 334 CALL nf95_inq_varid(fID_in_m, vars_in(1), vID_in(iv)) 335 IF(l3D) call NF95_GET_VAR(fID_in_m, vID_in(iv), o3_in3(1:nlon_in, :, :, 0, iv), sta, cnt)336 IF(l2d) call NF95_GET_VAR(fID_in_m, vID_in(iv), o3_in2(:, :, 0, iv), sta, cnt)335 IF(l3D) call nf95_get_var(fID_in_m, vID_in(iv), o3_in3(1:nlon_in, :, :, 0, iv), sta, cnt) 336 IF(l2d) call nf95_get_var(fID_in_m, vID_in(iv), o3_in2(:, :, 0, iv), sta, cnt) 337 337 END IF 338 338 IF(lnext) THEN; sta(ndims) = 1 339 339 CALL nf95_inq_varid(fID_in_p, vars_in(1), vID_in(iv)) 340 IF(l3D) call NF95_GET_VAR(fID_in_p, vID_in(iv), o3_in3(1:nlon_in, :, :, 13, iv), sta, cnt)341 IF(l2D) call NF95_GET_VAR(fID_in_p, vID_in(iv), o3_in2(:, :, 13, iv), sta, cnt)340 IF(l3D) call nf95_get_var(fID_in_p, vID_in(iv), o3_in3(1:nlon_in, :, :, 13, iv), sta, cnt) 341 IF(l2D) call nf95_get_var(fID_in_p, vID_in(iv), o3_in2(:, :, 13, iv), sta, cnt) 342 342 END IF 343 343 END DO 344 344 IF(lprev.OR.lnext) DEALLOCATE(sta, cnt) 345 IF(lprev) CALL NF95_CLOSE(fID_in_m)346 IF(lnext) CALL NF95_CLOSE(fID_in_p)345 IF(lprev) CALL nf95_close(fID_in_m) 346 IF(lnext) CALL nf95_close(fID_in_p) 347 347 348 348 !--- Revert decreasing coordinates vector … … 469 469 IF (is_mpi_root) THEN 470 470 !--- Write remaining coordinate variables: 471 CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)472 IF(interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)473 IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)471 CALL nf95_put_var(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in) 472 IF(interpt) CALL nf95_put_var(fID_ou, timID_ou, tmidday) 473 IF(.NOT.interpt) CALL nf95_put_var(fID_ou, timID_ou, tmidmonth) 474 474 475 475 !--- Write to file (the order of "rlatu" is inverted in the output file): … … 482 482 483 483 DO m = 1, read_climoz 484 CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m))484 CALL nf95_put_var(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m)) 485 485 END DO 486 486 487 487 ELSE 488 488 DO m = 1, read_climoz 489 CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out3(:, nlat_ou:1:-1, :, :, m))489 CALL nf95_put_var(fID_ou, vID_ou(m), o3_out3(:, nlat_ou:1:-1, :, :, m)) 490 490 END DO 491 491 ENDIF 492 CALL NF95_CLOSE(fID_ou)492 CALL nf95_close(fID_ou) 493 493 494 494 ENDIF … … 552 552 553 553 !--- Write remaining coordinate variables: 554 CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)555 IF(interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)556 IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)554 CALL nf95_put_var(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in) 555 IF(interpt) CALL nf95_put_var(fID_ou, timID_ou, tmidday) 556 IF(.NOT.interpt) CALL nf95_put_var(fID_ou, timID_ou, tmidmonth) 557 557 558 558 IF (grid_type==unstructured) THEN … … 564 564 565 565 DO m = 1, read_climoz 566 CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m))566 CALL nf95_put_var(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m)) 567 567 END DO 568 568 ELSE 569 569 !--- Write to file (the order of "rlatu" is inverted in the output file): 570 570 DO m = 1, read_climoz 571 CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out2(nlat_ou:1:-1, :, :, m))571 CALL nf95_put_var(fID_ou, vID_ou(m), o3_out2(nlat_ou:1:-1, :, :, m)) 572 572 END DO 573 573 ENDIF 574 574 575 CALL NF95_CLOSE(fID_ou)575 CALL nf95_close(fID_ou) 576 576 577 577 ENDIF … … 581 581 !============================================================================= 582 582 583 IF (is_mpi_root) CALL NF95_CLOSE(fID_in)583 IF (is_mpi_root) CALL nf95_close(fID_in) 584 584 585 585 ENDIF ! is_omp_master … … 599 599 ! dimensions and variables, and writes some of the coordinate variables. 600 600 !------------------------------------------------------------------------------- 601 USE regular_lonlat_mod, ONLY 602 USE regular_lonlat_mod, ONLY 603 USE mod_phys_lmdz_para, ONLY 604 USE mod_grid_phy_lmdz, ONLY 601 USE regular_lonlat_mod, ONLY: lon_reg, lat_reg 602 USE regular_lonlat_mod, ONLY: lon_reg, lat_reg 603 USE mod_phys_lmdz_para, ONLY: is_mpi_root 604 USE mod_grid_phy_lmdz, ONLY: klon_glo 605 605 606 606 !------------------------------------------------------------------------------- … … 625 625 sub = "prepare_out" 626 626 WRITE(lunout, *)"CALL sequence information: " // TRIM(sub) 627 CALL NF95_CREATE("climoz_LMDZ.nc", nf90_clobber, fID_ou)627 CALL nf95_create("climoz_LMDZ.nc", nf90_clobber, fID_ou) 628 628 629 629 !--- Dimensions: … … 641 641 CALL nf95_def_var(fID_ou, "time", nf90_float, dtimID, vtimID) 642 642 IF(ndims==4) & 643 CALL NF95_PUT_ATT(fID_ou, vlonID, "units", "degrees_east")644 CALL NF95_PUT_ATT(fID_ou, vlatID, "units", "degrees_north")645 CALL NF95_PUT_ATT(fID_ou, vlevID, "units", "millibar")646 CALL NF95_PUT_ATT(fID_ou, vtimID, "units", "days since 2000-1-1")643 CALL nf95_put_att(fID_ou, vlonID, "units", "degrees_east") 644 CALL nf95_put_att(fID_ou, vlatID, "units", "degrees_north") 645 CALL nf95_put_att(fID_ou, vlevID, "units", "millibar") 646 CALL nf95_put_att(fID_ou, vtimID, "units", "days since 2000-1-1") 647 647 IF(ndims==4) & 648 CALL NF95_PUT_ATT(fID_ou, vlonID, "standard_name", "longitude")649 CALL NF95_PUT_ATT(fID_ou, vlatID, "standard_name", "latitude")650 CALL NF95_PUT_ATT(fID_ou, vlevID, "standard_name", "air_pressure")651 CALL NF95_PUT_ATT(fID_ou, vtimID, "standard_name", "time")652 CALL NF95_PUT_ATT(fID_ou, vlevID, "long_name", "air pressure")653 CALL NF95_PUT_ATT(fID_ou, vtimID, "calendar", cal_ou)648 CALL nf95_put_att(fID_ou, vlonID, "standard_name", "longitude") 649 CALL nf95_put_att(fID_ou, vlatID, "standard_name", "latitude") 650 CALL nf95_put_att(fID_ou, vlevID, "standard_name", "air_pressure") 651 CALL nf95_put_att(fID_ou, vtimID, "standard_name", "time") 652 CALL nf95_put_att(fID_ou, vlevID, "long_name", "air pressure") 653 CALL nf95_put_att(fID_ou, vtimID, "calendar", cal_ou) 654 654 655 655 !--- Define the main variables: … … 657 657 IF(ndims==4) dIDs = [dlonID, dlatID, dlevID, dtimID] 658 658 CALL nf95_def_var(fID_ou, vars_in(1), nf90_float, dIDs(1:ndims), vID_ou(1)) 659 CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "long_name", "ozone mole fraction")660 CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "standard_name", "mole_fraction_of_ozone&659 CALL nf95_put_att(fID_ou, vID_ou(1), "long_name", "ozone mole fraction") 660 CALL nf95_put_att(fID_ou, vID_ou(1), "standard_name", "mole_fraction_of_ozone& 661 661 _in_air") 662 662 IF(SIZE(vID_ou) == 2) THEN 663 663 CALL nf95_def_var(fID_ou, vars_in(2), nf90_float, dIDs(1:ndims), vID_ou(2)) 664 CALL NF95_PUT_ATT(fID_ou, vID_ou(2), "long_name", "ozone mole fraction in da&664 CALL nf95_put_att(fID_ou, vID_ou(2), "long_name", "ozone mole fraction in da& 665 665 ylight") 666 666 END IF … … 669 669 ! The following commands, copying attributes, may fail. That is OK. 670 670 ! It should just mean that the attribute is not defined in the input file. 671 CALL NF95_COPY_ATT(fID_in, nf90_global, "Conventions", fID_ou, nf90_global, ncerr)671 CALL nf95_copy_att(fID_in, nf90_global, "Conventions", fID_ou, nf90_global, ncerr) 672 672 CALL handle_err_copy_att("Conventions") 673 CALL NF95_COPY_ATT(fID_in, nf90_global, "title", fID_ou, nf90_global, ncerr)673 CALL nf95_copy_att(fID_in, nf90_global, "title", fID_ou, nf90_global, ncerr) 674 674 CALL handle_err_copy_att("title") 675 CALL NF95_COPY_ATT(fID_in, nf90_global, "institution", fID_ou, nf90_global, ncerr)675 CALL nf95_copy_att(fID_in, nf90_global, "institution", fID_ou, nf90_global, ncerr) 676 676 CALL handle_err_copy_att("institution") 677 CALL NF95_COPY_ATT(fID_in, nf90_global, "source", fID_ou, nf90_global, ncerr)677 CALL nf95_copy_att(fID_in, nf90_global, "source", fID_ou, nf90_global, ncerr) 678 678 CALL handle_err_copy_att("source") 679 CALL NF95_PUT_ATT(fID_ou, nf90_global, "comment", "Regridded for LMDZ")680 CALL NF95_ENDDEF(fID_ou)679 CALL nf95_put_att (fID_ou, nf90_global, "comment", "Regridded for LMDZ") 680 CALL nf95_enddef(fID_ou) 681 681 682 682 IF (grid_type==unstructured) THEN … … 685 685 latitude_glo_(ind_cell_glo_glo(i)) = latitude_glo(i) 686 686 ENDDO 687 CALL NF95_PUT_VAR(fID_ou, vlatID, latitude_glo_)687 CALL nf95_put_var(fID_ou, vlatID, latitude_glo_) 688 688 ELSE 689 689 !--- Write one of the coordinate variables: 690 IF(ndims==4) CALL NF95_PUT_VAR(fID_ou, vlonID, lon_reg / deg2rad)691 CALL NF95_PUT_VAR(fID_ou, vlatID, lat_reg(nlat_ou:1:-1) / deg2rad)690 IF(ndims==4) CALL nf95_put_var(fID_ou, vlonID, lon_reg / deg2rad) 691 CALL nf95_put_var(fID_ou, vlatID, lat_reg(nlat_ou:1:-1) / deg2rad) 692 692 ! (convert from rad to degrees and sort in ascending order) 693 693 ENDIF … … 701 701 702 702 !------------------------------------------------------------------------------- 703 USE netcdf, ONLY : nf90_noerr, NF90_strerror703 USE netcdf, ONLY: nf90_noerr, nf90_strerror 704 704 !------------------------------------------------------------------------------- 705 705 ! Arguments: … … 707 707 !------------------------------------------------------------------------------- 708 708 IF(ncerr /= nf90_noerr) & 709 WRITE(lunout, *)TRIM(sub) // " prepare_out NF95_COPY_ATT" // TRIM(att_name) // &710 " -- " // TRIM( NF90_strerror(ncerr))709 WRITE(lunout, *)TRIM(sub) // " prepare_out nf95_copy_att " // TRIM(att_name) // & 710 " -- " // TRIM(nf90_strerror(ncerr)) 711 711 712 712 END SUBROUTINE handle_err_copy_att -
LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90
r5088 r5101 40 40 ! when we regrid in pressure). 41 41 42 use mod_grid_phy_lmdz, ONLY 42 use mod_grid_phy_lmdz, ONLY: nbp_lat 43 43 use regr_conserv_m, only: regr_conserv 44 44 use regr_lint_m, only: regr_lint … … 240 240 ! dimensions and variables, and writes one of the coordinate variables. 241 241 242 use mod_grid_phy_lmdz, ONLY 242 use mod_grid_phy_lmdz, ONLY: nbp_lat 243 243 use assert_eq_m, only: assert_eq 244 244 … … 247 247 use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global 248 248 use nrtype, only: pi 249 use regular_lonlat_mod, only: lat_reg249 use regular_lonlat_mod, ONLY: lat_reg 250 250 251 251 integer, intent(in):: ncid_in, varid_in(:), n_plev -
LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_time_av_m.F90
r5100 r5101 113 113 !------------------------------------------------------------------------------- 114 114 USE dimphy, ONLY: klon 115 USE netcdf95, ONLY: nf95_inq_varid, NF95_INQUIRE_VARIABLE, &115 USE netcdf95, ONLY: nf95_inq_varid, nf95_inquire_variable, & 116 116 nf95_inq_dimid, nf95_inquire_dimension, nf95_get_var 117 USE netcdf, ONLY: NF90_INQ_VARID, nf90_noerr117 USE netcdf, ONLY: nf90_inq_varid, nf90_noerr 118 118 USE assert_m, ONLY: assert 119 119 USE assert_eq_m, ONLY: assert_eq … … 210 210 !=== CHECK WHICH FIELDS ARE AVAILABLE IN THE INPUT FILE 211 211 IF(lfirst) THEN 212 lPrSfile=lAdjTro.AND. NF90_INQ_VARID(fID,"ps" ,vID)==nf90_noerr213 lPrTfile=lAdjTro.AND. NF90_INQ_VARID(fID,"tropopause_air_pressure",vID)==nf90_noerr214 lO3Tfile=lAdjTro.AND. NF90_INQ_VARID(fID,"tro3_at_tropopause" ,vID)==nf90_noerr212 lPrSfile=lAdjTro.AND.nf90_inq_varid(fID,"ps" ,vID)==nf90_noerr 213 lPrTfile=lAdjTro.AND.nf90_inq_varid(fID,"tropopause_air_pressure",vID)==nf90_noerr 214 lO3Tfile=lAdjTro.AND.nf90_inq_varid(fID,"tro3_at_tropopause" ,vID)==nf90_noerr 215 215 CALL nf95_inq_dimid(fID,"time",vID) 216 216 CALL nf95_inquire_dimension(fID,vID,nclen=ntim_in) … … 481 481 !------------------------------------------------------------------------------- 482 482 CALL nf95_inq_varid(fID, TRIM(var), vID) 483 CALL NF95_INQUIRE_VARIABLE(fID, vID, ndims=n_dim)484 IF(n_dim==2) call NF95_GET_VAR(fID,vID,v(1,:), start=[ 1,irec])485 IF(n_dim==3) call NF95_GET_VAR(fID,vID,v(:,:), start=[1,1,irec])483 CALL nf95_inquire_variable(fID, vID, ndims=n_dim) 484 IF(n_dim==2) call nf95_get_var(fID,vID,v(1,:), start=[ 1,irec]) 485 IF(n_dim==3) call nf95_get_var(fID,vID,v(:,:), start=[1,1,irec]) 486 486 487 487 !--- Flip latitudes: ascending in input file, descending in "rlatu". … … 512 512 DO i=1,SIZE(nam) 513 513 CALL nf95_inq_varid(fID, TRIM(nam(i)), vID) 514 CALL NF95_INQUIRE_VARIABLE(fID, vID, ndims=n_dim)515 IF(n_dim==3) call NF95_GET_VAR(fID,vID,v(1,:,:,i), start=[ 1,1,irec])516 IF(n_dim==4) call NF95_GET_VAR(fID,vID,v(:,:,:,i), start=[1,1,1,irec])514 CALL nf95_inquire_variable(fID, vID, ndims=n_dim) 515 IF(n_dim==3) call nf95_get_var(fID,vID,v(1,:,:,i), start=[ 1,1,irec]) 516 IF(n_dim==4) call nf95_get_var(fID,vID,v(:,:,:,i), start=[1,1,1,irec]) 517 517 END DO 518 518 -
LMDZ6/branches/Amaury_dev/libf/phylmd/slab_heat_transp_mod.F90
r5099 r5101 146 146 147 147 ! USE ioipsl_getin_p_mod, only: getin_p 148 USE IOIPSL, ONLY : getin148 USE IOIPSL, ONLY: getin 149 149 IMPLICIT NONE 150 150 -
LMDZ6/branches/Amaury_dev/libf/phylmd/stdlevvar_mod.F90
r5099 r5101 319 319 n2mout) 320 320 321 USE ioipsl_getin_p_mod, ONLY 321 USE ioipsl_getin_p_mod, ONLY: getin_p 322 322 IMPLICIT NONE 323 323 !------------------------------------------------------------------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/sumethox.F90
r5087 r5101 45 45 ! ------------------------------------------------------------------ 46 46 47 USE YOEMETH , ONLY 47 USE YOEMETH , ONLY: RALPHA1 ,RALPHA2 ,RQLIM ,& 48 48 RPBOTOX, RPBOTPH ,RPTOPOX ,RPTOPPH ,& 49 49 RALPHA3, RLOGPPH -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_mod.F90
r5099 r5101 23 23 USE intersurf ! module in ORCHIDEE 24 24 #endif 25 USE cpl_mod, ONLY 26 USE surface_data, ONLY 27 USE geometry_mod, ONLY 25 USE cpl_mod, ONLY: cpl_send_land_fields, cpl_send_landice_fields 26 USE surface_data, ONLY: type_ocean, landice_opt 27 USE geometry_mod, ONLY: dx, dy, boundslon, boundslat,longitude, latitude, cell_area, ind_cell_glo 28 28 USE mod_grid_phy_lmdz 29 29 USE mod_phys_lmdz_para, mpi_root_rank=>mpi_master 30 USE carbon_cycle_mod, ONLY 31 USE nrtype, ONLY 30 USE carbon_cycle_mod, ONLY: nbcf_in_orc, nbcf_out, fields_in, yfields_in, yfields_out, cfname_in, cfname_out 31 USE nrtype, ONLY: PI 32 32 33 33 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_nofrein_mod.F90
r5099 r5101 23 23 USE intersurf ! module d'ORCHIDEE 24 24 #endif 25 USE cpl_mod, ONLY 26 USE surface_data, ONLY 27 USE geometry_mod, ONLY 25 USE cpl_mod, ONLY: cpl_send_land_fields 26 USE surface_data, ONLY: type_ocean 27 USE geometry_mod, ONLY: dx, dy 28 28 USE mod_grid_phy_lmdz 29 29 USE mod_phys_lmdz_para, mpi_root_rank=>mpi_master … … 53 53 USE mod_surf_para 54 54 USE mod_synchro_omp 55 USE carbon_cycle_mod, ONLY 55 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 56 56 USE indice_sol_mod 57 57 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_nolic_mod.F90
r5099 r5101 19 19 USE intersurf ! module d'ORCHIDEE 20 20 #endif 21 USE cpl_mod, ONLY 22 USE surface_data, ONLY 23 USE geometry_mod, ONLY 21 USE cpl_mod, ONLY: cpl_send_land_fields 22 USE surface_data, ONLY: type_ocean 23 USE geometry_mod, ONLY: dx, dy, boundslon, boundslat,longitude, latitude, cell_area, ind_cell_glo 24 24 USE mod_grid_phy_lmdz 25 25 USE mod_phys_lmdz_para, mpi_root_rank=>mpi_master 26 USE carbon_cycle_mod, ONLY 27 USE nrtype, ONLY 26 USE carbon_cycle_mod, ONLY: nbcf_in_orc, nbcf_out, fields_in, yfields_in, yfields_out, cfname_in, cfname_out 27 USE nrtype, ONLY: PI 28 28 29 29 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_noopenmp_mod.F90
r5099 r5101 20 20 USE intersurf ! module d'ORCHIDEE 21 21 #endif 22 USE cpl_mod, ONLY 23 USE surface_data, ONLY 24 USE geometry_mod, ONLY 22 USE cpl_mod, ONLY: cpl_send_land_fields 23 USE surface_data, ONLY: type_ocean 24 USE geometry_mod, ONLY: dx, dy 25 25 USE mod_grid_phy_lmdz 26 26 USE mod_phys_lmdz_para … … 99 99 ! qsurf air moisture at surface 100 100 101 USE carbon_cycle_mod, ONLY 101 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl, fco2_land_inst, fco2_lu_inst 102 102 USE indice_sol_mod 103 103 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_nounstruct_mod.F90
r5099 r5101 19 19 USE intersurf ! module d'ORCHIDEE 20 20 #endif 21 USE cpl_mod, ONLY 22 USE surface_data, ONLY 23 USE geometry_mod, ONLY 21 USE cpl_mod, ONLY: cpl_send_land_fields 22 USE surface_data, ONLY: type_ocean 23 USE geometry_mod, ONLY: dx, dy 24 24 USE mod_grid_phy_lmdz 25 25 USE mod_phys_lmdz_para, mpi_root_rank=>mpi_master … … 49 49 USE mod_surf_para 50 50 USE mod_synchro_omp 51 USE carbon_cycle_mod, ONLY 51 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 52 52 USE indice_sol_mod 53 53 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_orchidee_noz0h_mod.F90
r5099 r5101 22 22 USE intersurf ! module d'ORCHIDEE 23 23 #endif 24 USE cpl_mod, ONLY 25 USE surface_data, ONLY 26 USE geometry_mod, ONLY 24 USE cpl_mod, ONLY: cpl_send_land_fields 25 USE surface_data, ONLY: type_ocean 26 USE geometry_mod, ONLY: dx, dy 27 27 USE mod_grid_phy_lmdz 28 28 USE mod_phys_lmdz_para, mpi_root_rank=>mpi_master … … 51 51 USE mod_surf_para 52 52 USE mod_synchro_omp 53 USE carbon_cycle_mod, ONLY 53 USE carbon_cycle_mod, ONLY: carbon_cycle_cpl 54 54 USE indice_sol_mod 55 55 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_landice_mod.F90
r5099 r5101 31 31 32 32 USE dimphy 33 USE geometry_mod, ONLY 34 USE surface_data, ONLY 35 USE fonte_neige_mod, ONLY 36 USE cpl_mod, ONLY 33 USE geometry_mod, ONLY: longitude,latitude 34 USE surface_data, ONLY: type_ocean, calice, calsno, landice_opt, iflag_albcalc 35 USE fonte_neige_mod, ONLY: fonte_neige,run_off_lic,fqcalving_global,ffonte_global,fqfonte_global,runofflic_global 36 USE cpl_mod, ONLY: cpl_send_landice_fields 37 37 USE calcul_fluxs_mod 38 USE phys_local_var_mod, ONLY 39 USE phys_output_var_mod, ONLY 38 USE phys_local_var_mod, ONLY: zxrhoslic, zxustartlic, zxqsaltlic 39 USE phys_output_var_mod, ONLY: snow_o,zfra_o 40 40 #ifdef ISO 41 USE fonte_neige_mod, ONLY 42 USE infotrac_phy, ONLY 41 USE fonte_neige_mod, ONLY: xtrun_off_lic 42 USE infotrac_phy, ONLY: ntiso,niso 43 43 USE isotopes_routines_mod, ONLY: calcul_iso_surf_lic_vectall 44 44 #ifdef ISOVERIF … … 49 49 50 50 !FC 51 USE ioipsl_getin_p_mod, ONLY 52 USE lmdz_blowing_snow_ini, ONLY 53 USE lmdz_blowing_snow_ini, ONLY 51 USE ioipsl_getin_p_mod, ONLY: getin_p 52 USE lmdz_blowing_snow_ini, ONLY: c_esalt_bs, zeta_bs, pbst_bs, prt_bs, rhoice_bs, rhohard_bs 53 USE lmdz_blowing_snow_ini, ONLY: rhofresh_bs, tau_eqsalt_bs, tau_dens0_bs, tau_densmin_bs 54 54 #ifdef CPP_INLANDSIS 55 USE surf_inlandsis_mod, ONLY 55 USE surf_inlandsis_mod, ONLY: surf_inlandsis 56 56 #endif 57 57 -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90
r5099 r5101 32 32 USE dimphy, ONLY: klon, zmasq 33 33 USE surface_data, ONLY : type_ocean 34 USE ocean_forced_mod, ONLY 34 USE ocean_forced_mod, ONLY: ocean_forced_noice 35 35 USE ocean_slab_mod, ONLY : ocean_slab_noice 36 36 USE ocean_cpl_mod, ONLY : ocean_cpl_noice 37 USE indice_sol_mod, ONLY 38 #ifdef ISO 39 USE infotrac_phy, ONLY 37 USE indice_sol_mod, ONLY: nbsrf, is_oce 38 #ifdef ISO 39 USE infotrac_phy, ONLY: ntraciso=>ntiso,niso 40 40 #ifdef ISOVERIF 41 41 USE isotopes_mod, ONLY: iso_eau,ridicule -
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_seaice_mod.F90
r5099 r5101 30 30 USE dimphy 31 31 USE surface_data 32 USE ocean_forced_mod, ONLY 32 USE ocean_forced_mod, ONLY: ocean_forced_ice 33 33 USE ocean_cpl_mod, ONLY : ocean_cpl_ice 34 34 USE ocean_slab_mod, ONLY : ocean_slab_ice 35 35 USE indice_sol_mod 36 36 #ifdef ISO 37 USE infotrac_phy, ONLY 37 USE infotrac_phy, ONLY: ntiso,niso 38 38 #endif 39 39 -
LMDZ6/branches/Amaury_dev/libf/phylmd/sw_aeroAR4.F90
r5099 r5101 21 21 22 22 USE dimphy 23 USE phys_output_mod, ONLY 23 USE phys_output_mod, ONLY: swaero_diag 24 24 USE print_control_mod, ONLY: lunout 25 USE aero_mod, ONLY 25 USE aero_mod, ONLY: naero_grp 26 26 IMPLICIT NONE 27 27 -
LMDZ6/branches/Amaury_dev/libf/phylmd/time_phylmdz_mod.F90
r5099 r5101 36 36 SUBROUTINE init_time(annee_ref_, day_ref_, day_ini_, start_time_, & 37 37 ndays_, pdtphys_) 38 USE ioipsl_getin_p_mod, ONLY 38 USE ioipsl_getin_p_mod, ONLY: getin_p 39 39 USE phys_cal_mod, ONLY: phys_cal_init 40 40 IMPLICIT NONE … … 79 79 SUBROUTINE update_time(pdtphys_) 80 80 ! This subroutine updates the module saved variables. 81 USE ioipsl, ONLY 81 USE ioipsl, ONLY: ymds2ju 82 82 USE phys_cal_mod, ONLY: phys_cal_update 83 83 USE print_control_mod, ONLY: lunout -
LMDZ6/branches/Amaury_dev/libf/phylmd/tracco2i_mod.F90
r5099 r5101 344 344 USE dimphy 345 345 ! USE infotrac_phy 346 USE geometry_mod, ONLY 346 USE geometry_mod, ONLY: cell_area 347 347 USE mod_grid_phy_lmdz 348 348 USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root … … 353 353 USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite 354 354 355 USE carbon_cycle_mod, ONLY 355 USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean 356 356 357 357 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/tracinca_mod.F90
r5099 r5101 61 61 USE geometry_mod, ONLY: cell_area 62 62 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 63 USE aero_mod, ONLY 63 USE aero_mod, ONLY: naero_grp 64 64 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 65 65 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/tracreprobus_mod.F90
r5099 r5101 13 13 USE infotrac_phy, ONLY: nbtr 14 14 #ifdef REPROBUS 15 USE CHEM_REP, ONLY 15 USE CHEM_REP, ONLY: pdt_rep, & ! pas de temps reprobus 16 16 daynum, iter, & ! jourjulien, iteration chimie 17 17 pdel,& -
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada4.F90
r5099 r5101 4 4 cd, tke, eps, km, kn, kq, ustar, iflag_pbl, drgpro) 5 5 6 USE dimphy, only: klev,klon6 USE dimphy, ONLY: klev,klon 7 7 USE phys_local_var_mod, only: wprime 8 USE yamada_ini_mod, only: new_yamada4,yamada4_num,hboville9 USE yamada_ini_mod, only: prt_level, lunout,pbl_lmixmin_alpha,b1,kap,viscom,viscoh10 USE yamada_ini_mod, only: ric, yun,ydeux,lmixmin,iflag_vdif_q28 USE yamada_ini_mod, ONLY: new_yamada4,yamada4_num,hboville 9 USE yamada_ini_mod, ONLY: prt_level, lunout,pbl_lmixmin_alpha,b1,kap,viscom,viscoh 10 USE yamada_ini_mod, ONLY: ric, yun,ydeux,lmixmin,iflag_vdif_q2 11 11 12 12 IMPLICIT NONE … … 712 712 SUBROUTINE vdif_q2(timestep, gravity, rconst, ngrid, plev, temp, kmy, q2) 713 713 714 USE dimphy, only: klev,klon714 USE dimphy, ONLY: klev,klon 715 715 IMPLICIT NONE 716 716 … … 801 801 SUBROUTINE vdif_q2e(timestep, gravity, rconst, ngrid, plev, temp, kmy, q2) 802 802 803 USE dimphy, only: klev,klon803 USE dimphy, ONLY: klev,klon 804 804 IMPLICIT NONE 805 805 … … 872 872 873 873 874 USE dimphy, only: klev,klon875 USE yamada_ini_mod, only: l0874 USE dimphy, ONLY: klev,klon 875 USE yamada_ini_mod, ONLY: l0 876 876 USE phys_state_var_mod, only: zstd, zsig, zmea 877 877 USE phys_local_var_mod, only: l_mixmin, l_mix 878 USE yamada_ini_mod, only: kap, kapb878 USE yamada_ini_mod, ONLY: kap, kapb 879 879 880 880 ! zstd: ecart type de la'altitud e sous-maille -
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.F90
r5099 r5101 7 7 USE dimphy, ONLY: klon, klev 8 8 USE print_control_mod, ONLY: prt_level 9 USE ioipsl_getin_p_mod, ONLY 9 USE ioipsl_getin_p_mod, ONLY: getin_p 10 10 11 11 IMPLICIT NONE -
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_ini_mod.F90
r5081 r5101 27 27 SUBROUTINE yamada_ini(klon,lunout_in,prt_level_in) 28 28 29 USE ioipsl_getin_p_mod, ONLY 29 USE ioipsl_getin_p_mod, ONLY: getin_p 30 30 31 31 ! local
Note: See TracChangeset
for help on using the changeset viewer.