Changeset 1390
- Timestamp:
- Mar 5, 2015, 12:44:03 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/newstart.F
r1383 r1390 35 35 use comgeomfi_h, only: ini_fillgeom 36 36 37 37 38 implicit none 38 39 … … 168 169 real :: Mair_old,Mair_new,vmr_old,vmr_new 169 170 real,allocatable :: coefvmr(:) ! Correction coefficient when changing composition 171 integer :: iloc(1), iqmax 170 172 171 173 c sortie visu pour les champs dynamiques … … 1021 1023 1022 1024 if (iq.eq.igcm_n2) then 1023 write(*,*) "New vmr(n2)? (MSL: 1.89E-02 at Ls~186)"1025 write(*,*) "New vmr(n2)? (MSL: 2.03e-02 at Ls~184)" 1024 1026 endif 1025 1027 if (iq.eq.igcm_ar) then 1026 write(*,*) "New vmr(ar)? (MSL: 1.93E-02 at Ls~186)"1028 write(*,*) "New vmr(ar)? (MSL: 2.07e-02 at Ls~184)" 1027 1029 endif 1028 1030 if (iq.eq.igcm_o2) then 1029 write(*,*) "New vmr(o2)? (MSL: 1. 46E-03) at Ls~186"1031 write(*,*) "New vmr(o2)? (MSL: 1.73e-03 at Ls~184)" 1030 1032 endif 1031 1033 if (iq.eq.igcm_co) then 1032 write(*,*) "New vmr(co)? ( ~8.E-04)"1034 write(*,*) "New vmr(co)? (MSL: 7.49e-04 at Ls~184)" 1033 1035 endif 1034 1036 302 read(*,*,iostat=ierr) vmr_new … … 1055 1057 end do 1056 1058 1057 ! Recompute mass mixing ratios everywhere, and adjust mmr(CO2) to keep sum of mmr constant. 1059 ! Recompute mass mixing ratios everywhere, and adjust mmr of most abundant species 1060 ! to keep sum of mmr constant. 1058 1061 do l=1,llm 1059 1062 do j=1,jjp1 … … 1069 1072 end if 1070 1073 enddo 1071 q(i,j,l,igcm_co2) = q(i,j,l,igcm_co2) + Smmr_old-Smmr_new 1074 iloc = maxloc(q(i,j,l,:)) 1075 iqmax = iloc(1) 1076 q(i,j,l,iqmax) = q(i,j,l,iqmax) + Smmr_old - Smmr_new 1072 1077 enddo 1073 1078 enddo … … 1075 1080 1076 1081 write(*,*) 1077 & "mmr(CO2) is modified everywhere to keep sum of mmr constant" 1082 & "The most abundant species is modified everywhere to keep "// 1083 & "sum of mmr constant" 1078 1084 write(*,*) 'At reference site vmr(CO2)=', 1079 1085 & q(iref,jref,1,igcm_co2)*Mair_new/mmol(igcm_co2) 1080 write(*,*) "Compared to MSL observation: vmr(CO2)= 0.95 97"1081 1086 write(*,*) "Compared to MSL observation: vmr(CO2)= 0.957 "// 1087 & "at Ls=184" 1082 1088 1083 1089 c wetstart : wet atmosphere with a north to south gradient
Note: See TracChangeset
for help on using the changeset viewer.