- Timestamp:
- Jul 24, 2024, 4:23:34 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/phytracr_spl_mod.F90
r5116 r5117 828 828 ! ------- 829 829 830 real, intent(in) :: pdtphys ! pas d'integration pour la physique (seconde)831 REAL, intent(in) :: jD_cur, jH_cur832 real, intent(in) :: ftsol(klon, nbsrf) ! temperature du sol par type833 real, intent(in) :: t_seri(klon, klev) ! temperature834 real, intent(in) :: u_seri(klon, klev) ! vent835 real, intent(in) :: v_seri(klon, klev) ! vent836 real, intent(in) :: q_seri(klon, klev) ! vapeur d eau kg/kg830 REAL, INTENT(IN) :: pdtphys ! pas d'integration pour la physique (seconde) 831 REAL, INTENT(IN) :: jD_cur, jH_cur 832 REAL, INTENT(IN) :: ftsol(klon, nbsrf) ! temperature du sol par type 833 REAL, INTENT(IN) :: t_seri(klon, klev) ! temperature 834 REAL, INTENT(IN) :: u_seri(klon, klev) ! vent 835 REAL, INTENT(IN) :: v_seri(klon, klev) ! vent 836 REAL, INTENT(IN) :: q_seri(klon, klev) ! vapeur d eau kg/kg 837 837 838 838 LOGICAL, INTENT(IN) :: lafin 839 839 840 realtr_seri(klon, klev, nbtr) ! traceur841 realtmp_var(klon, klev) ! auxiliary variable to replace traceur842 realtmp_var2(klon, nbtr) ! auxiliary variable to replace source843 realtmp_var3(klon, klev, nbtr) ! auxiliary variable 3D844 realdummy1d ! JE auxiliary variable845 realaux_var2(klon) ! auxiliary variable to replace traceur846 realaux_var3(klon, klev) ! auxiliary variable to replace traceur847 reald_tr(klon, klev, nbtr) ! traceur tendance848 realsconc_seri(klon, nbtr) ! surface concentration of traceur849 850 integernbjour840 REAL tr_seri(klon, klev, nbtr) ! traceur 841 REAL tmp_var(klon, klev) ! auxiliary variable to replace traceur 842 REAL tmp_var2(klon, nbtr) ! auxiliary variable to replace source 843 REAL tmp_var3(klon, klev, nbtr) ! auxiliary variable 3D 844 REAL dummy1d ! JE auxiliary variable 845 REAL aux_var2(klon) ! auxiliary variable to replace traceur 846 REAL aux_var3(klon, klev) ! auxiliary variable to replace traceur 847 REAL d_tr(klon, klev, nbtr) ! traceur tendance 848 REAL sconc_seri(klon, nbtr) ! surface concentration of traceur 849 850 INTEGER nbjour 851 851 save nbjour 852 852 !$OMP THREADPRIVATE(nbjour) … … 898 898 !JE20150518>> 899 899 900 real, intent(in) :: paprs(klon, klev + 1) ! pression pour chaque inter-couche (en Pa)901 real, intent(in) :: pplay(klon, klev) ! pression pour le mileu de chaque couche (en Pa)902 real, intent(in) :: RHcl(klon, klev) ! humidite relativen ciel clair903 realznivsig(klev) ! indice des couches904 realpaire(klon)905 real, intent(in) :: pphis(klon)906 real, intent(in) :: pctsrf(klon, nbsrf)907 logical, intent(in) :: debutphy ! le flag de l'initialisation de la physique900 REAL, INTENT(IN) :: paprs(klon, klev + 1) ! pression pour chaque inter-couche (en Pa) 901 REAL, INTENT(IN) :: pplay(klon, klev) ! pression pour le mileu de chaque couche (en Pa) 902 REAL, INTENT(IN) :: RHcl(klon, klev) ! humidite relativen ciel clair 903 REAL znivsig(klev) ! indice des couches 904 REAL paire(klon) 905 REAL, INTENT(IN) :: pphis(klon) 906 REAL, INTENT(IN) :: pctsrf(klon, nbsrf) 907 logical, INTENT(IN) :: debutphy ! le flag de l'initialisation de la physique 908 908 909 909 ! Scaling Parameters: … … 950 950 ! ----------- 951 951 952 REAL, intent(in) :: pmfu(klon, klev) ! flux de masse dans le panache montant953 REAL, intent(in) :: pmfd(klon, klev) ! flux de masse dans le panache descendant954 REAL, intent(in) :: pen_u(klon, klev) ! flux entraine dans le panache montant955 REAL, intent(in) :: pde_u(klon, klev) ! flux detraine dans le panache montant956 REAL, intent(in) :: pen_d(klon, klev) ! flux entraine dans le panache descendant957 REAL, intent(in) :: pde_d(klon, klev) ! flux detraine dans le panache descendant952 REAL, INTENT(IN) :: pmfu(klon, klev) ! flux de masse dans le panache montant 953 REAL, INTENT(IN) :: pmfd(klon, klev) ! flux de masse dans le panache descendant 954 REAL, INTENT(IN) :: pen_u(klon, klev) ! flux entraine dans le panache montant 955 REAL, INTENT(IN) :: pde_u(klon, klev) ! flux detraine dans le panache montant 956 REAL, INTENT(IN) :: pen_d(klon, klev) ! flux entraine dans le panache descendant 957 REAL, INTENT(IN) :: pde_d(klon, klev) ! flux detraine dans le panache descendant 958 958 959 959 ! Convection KE scheme: … … 1023 1023 ! --------- 1024 1024 1025 REAL, intent(in) :: pmflxr(klon, klev + 1), pmflxs(klon, klev + 1) !--convection1026 REAL, intent(in) :: prfl(klon, klev + 1), psfl(klon, klev + 1) !--large-scale1025 REAL, INTENT(IN) :: pmflxr(klon, klev + 1), pmflxs(klon, klev + 1) !--convection 1026 REAL, INTENT(IN) :: prfl(klon, klev + 1), psfl(klon, klev + 1) !--large-scale 1027 1027 REAL :: ql_incl ! contenu en eau liquide nuageuse dans le nuage ! ql_incl=oliq/rneb 1028 1028 REAL :: ql_incloud_ref ! ref value of in-cloud condensed water content … … 1047 1047 ! -------------- 1048 1048 1049 REAL, intent(in) :: coefh(klon, klev) ! coeff melange CL1050 REAL, intent(in) :: cdragh(klon), cdragm(klon)1051 REAL, intent(in) :: yu1(klon) ! vent dans la 1iere couche1052 REAL, intent(in) :: yv1(klon) ! vent dans la 1iere couche1049 REAL, INTENT(IN) :: coefh(klon, klev) ! coeff melange CL 1050 REAL, INTENT(IN) :: cdragh(klon), cdragm(klon) 1051 REAL, INTENT(IN) :: yu1(klon) ! vent dans la 1iere couche 1052 REAL, INTENT(IN) :: yv1(klon) ! vent dans la 1iere couche 1053 1053 1054 1054 … … 1071 1071 ! ----------- 1072 1072 1073 REAL, intent(in) :: rlat(klon) ! latitudes pour chaque point1074 REAL, intent(in) :: rlon(klon) ! longitudes pour chaque point1073 REAL, INTENT(IN) :: rlat(klon) ! latitudes pour chaque point 1074 REAL, INTENT(IN) :: rlon(klon) ! longitudes pour chaque point 1075 1075 1076 1076 INTEGER i, k, iq, itr, j, ig … … 1164 1164 ! REAL his_g2paer(klon) ! gastoparticle in aerosol units (check!) 1165 1165 1166 INTEGER, intent(in) :: iflag_conv1166 INTEGER, INTENT(IN) :: iflag_conv 1167 1167 LOGICAL iscm3 ! debug variable. for checkmass ! JE 1168 1168 … … 1213 1213 source_tr = 0. 1214 1214 1215 if(debutphy) THEN1215 IF (debutphy) THEN 1216 1216 #ifdef IOPHYS_DUST 1217 1217 CALL iophys_ini(pdtphys) … … 1394 1394 enddo 1395 1395 ! check consistency with dust emission scheme: 1396 if(ok_chimeredust) THEN1397 if (.not.(id_scdu>0 .and. id_codu>0 .and. id_fine>0)) THEN1396 IF (ok_chimeredust) THEN 1397 IF (.NOT.(id_scdu>0 .AND. id_codu>0 .AND. id_fine>0)) THEN 1398 1398 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 0', 1) 1399 1399 endif 1400 1400 else 1401 if(id_scdu>0) THEN1401 IF (id_scdu>0) THEN 1402 1402 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU', 1) 1403 1403 endif 1404 if ((id_codu <= 0) .or. (id_fine<=0)) THEN1404 IF ((id_codu <= 0) .OR. (id_fine<=0)) THEN 1405 1405 CALL abort_gcm('phytracr_mod', 'pb in ok_chimdust 1', 1) 1406 1406 endif … … 1521 1521 1522 1522 iscm3 = .FALSE. 1523 if(debutphy) THEN1523 IF (debutphy) THEN 1524 1524 !$OMP MASTER 1525 1525 CALL suphel … … 4442 4442 REAL scale_param(nbreg) 4443 4443 !local vars 4444 integernid, ierr, nvarid4445 realrcode, auxreal4446 integerstart(4), count(4), status4444 INTEGER nid, ierr, nvarid 4445 REAL rcode, auxreal 4446 INTEGER start(4), count(4), status 4447 4447 ! local 4448 4448 CHARACTER*104 varname 4449 4449 CHARACTER*2 aux_2s 4450 integeri, j, ig4450 INTEGER i, j, ig 4451 4451 !$OMP MASTER 4452 4452 IF (is_mpi_root .AND. is_omp_root) THEN 4453 4453 ierr = nf90_open(trim(adjustl(filescaleparams)), nf90_nowrite, nid) 4454 if(ierr == nf90_noerr) THEN4454 IF (ierr == nf90_noerr) THEN 4455 4455 do i = 1, nbreg 4456 4456 WRITE(aux_2s, '(i2.2)') i
Note: See TracChangeset
for help on using the changeset viewer.
