Changeset 556 for trunk/LMDZ.MARS/libf
- Timestamp:
- Mar 2, 2012, 4:07:12 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/physiq.F
r550 r556 1059 1059 END IF ! of IF (water) 1060 1060 1061 1062 c 7b. Chemical species 1063 c ------------------ 1064 1065 #ifndef MESOSCALE 1066 c -------------- 1067 c photochemistry : 1068 c -------------- 1069 IF (photochem .or. thermochem) then 1070 !NB: Photochemistry includes condensation of H2O2 1071 PRINT*, 'SURFDUST,SURFICE TO BE IMPLEMENTED. YAAAAAARG.' 1072 call calchim(ptimestep,pplay,pplev,pt,pdt,dist_sol,mu0, 1073 $ zzlev,zzlay,zday,pq,pdq,zdqchim,zdqschim, 1074 $ zdqcloud,zdqscloud,tauref,co2ice, 1075 $ pu,pdu,pv,pdv,surfdust,surfice) 1076 1077 ! increment values of tracers: 1078 DO iq=1,nq ! loop on all tracers; tendencies for non-chemistry 1079 ! tracers is zero anyways 1080 DO l=1,nlayer 1081 DO ig=1,ngrid 1082 pdq(ig,l,iq)=pdq(ig,l,iq)+zdqchim(ig,l,iq) 1083 ENDDO 1084 ENDDO 1085 ENDDO ! of DO iq=1,nq 1086 ! add condensation tendency for H2O2 1087 if (igcm_h2o2.ne.0) then 1088 DO l=1,nlayer 1089 DO ig=1,ngrid 1090 pdq(ig,l,igcm_h2o2)=pdq(ig,l,igcm_h2o2) 1091 & +zdqcloud(ig,l,igcm_h2o2) 1092 ENDDO 1093 ENDDO 1094 endif 1095 1096 ! increment surface values of tracers: 1097 DO iq=1,nq ! loop on all tracers; tendencies for non-chemistry 1098 ! tracers is zero anyways 1099 DO ig=1,ngrid 1100 dqsurf(ig,iq)=dqsurf(ig,iq)+zdqschim(ig,iq) 1101 ENDDO 1102 ENDDO ! of DO iq=1,nq 1103 ! add condensation tendency for H2O2 1104 if (igcm_h2o2.ne.0) then 1105 DO ig=1,ngrid 1106 dqsurf(ig,igcm_h2o2)=dqsurf(ig,igcm_h2o2) 1107 & +zdqscloud(ig,igcm_h2o2) 1108 ENDDO 1109 endif 1110 1111 END IF ! of IF (photochem.or.thermochem) 1112 #endif 1113 1114 c 7c. Aerosol particles 1061 c 7b. Aerosol particles 1115 1062 c ------------------- 1116 1063 … … 1171 1118 END IF ! of IF (sedimentation) 1172 1119 1173 1120 c 1121 c 7c. Chemical species 1122 c ------------------ 1123 1124 #ifndef MESOSCALE 1125 c -------------- 1126 c photochemistry : 1127 c -------------- 1128 IF (photochem .or. thermochem) then 1129 1130 ! dust and ice surface area 1131 call surfacearea(ngrid, nlayer, ptimestep, pplay, zzlay, 1132 $ pt, pq, pdq, nq, 1133 $ rdust, rice, tau, tauscaling, 1134 $ surfdust, surfice) 1135 ! call photochemistry 1136 call calchim(ptimestep,pplay,pplev,pt,pdt,dist_sol,mu0, 1137 $ zzlev,zzlay,zday,pq,pdq,zdqchim,zdqschim, 1138 $ zdqcloud,zdqscloud,tauref,co2ice, 1139 $ pu,pdu,pv,pdv,surfdust,surfice) 1140 1141 ! increment values of tracers: 1142 DO iq=1,nq ! loop on all tracers; tendencies for non-chemistry 1143 ! tracers is zero anyways 1144 DO l=1,nlayer 1145 DO ig=1,ngrid 1146 pdq(ig,l,iq)=pdq(ig,l,iq)+zdqchim(ig,l,iq) 1147 ENDDO 1148 ENDDO 1149 ENDDO ! of DO iq=1,nq 1150 1151 ! add condensation tendency for H2O2 1152 if (igcm_h2o2.ne.0) then 1153 DO l=1,nlayer 1154 DO ig=1,ngrid 1155 pdq(ig,l,igcm_h2o2)=pdq(ig,l,igcm_h2o2) 1156 & +zdqcloud(ig,l,igcm_h2o2) 1157 ENDDO 1158 ENDDO 1159 endif 1160 1161 ! increment surface values of tracers: 1162 DO iq=1,nq ! loop on all tracers; tendencies for non-chemistry 1163 ! tracers is zero anyways 1164 DO ig=1,ngrid 1165 dqsurf(ig,iq)=dqsurf(ig,iq)+zdqschim(ig,iq) 1166 ENDDO 1167 ENDDO ! of DO iq=1,nq 1168 1169 ! add condensation tendency for H2O2 1170 if (igcm_h2o2.ne.0) then 1171 DO ig=1,ngrid 1172 dqsurf(ig,igcm_h2o2)=dqsurf(ig,igcm_h2o2) 1173 & +zdqscloud(ig,igcm_h2o2) 1174 ENDDO 1175 endif 1176 1177 END IF ! of IF (photochem.or.thermochem) 1178 #endif 1174 1179 1175 1180 c 7d. Updates
Note: See TracChangeset
for help on using the changeset viewer.