Changeset 1684 for trunk/LMDZ.MARS/libf
- Timestamp:
- Mar 31, 2017, 5:03:31 PM (8 years ago)
- Location:
- trunk/LMDZ.MARS/libf
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/chemthermos.F90 ¶
r1266 r1684 466 466 467 467 !Photoabsorption coefficients 468 call flujo(solarcondate) 469 if(solvarmod.eq.0) then 470 call jthermcalc(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit) 471 else if(solvarmod.eq.1) then 472 call jthermcalc_e107(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit,zday) 473 endif 468 call jthermcalc_e107(ig,nlayer,chemthermod,rm,nesptherm,ztemp,zlocal,zenit,zday) 474 469 475 470 -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90 ¶
r1266 r1684 359 359 enddo 360 360 361 !Solar flux calculation362 call flujo(solarcondate)363 364 ! Not recommended for long runs365 ! (e.g. to build the MCD), as the366 ! solar conditions at the end will367 ! be different to the ones initially368 ! set369 361 370 362 do ig=1,ngrid -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/hrtherm.F ¶
r1266 r1684 98 98 99 99 !Calculation of photoabsortion coefficient 100 if(solvarmod.eq.0) then 101 call jthermcalc(ig,nlayer,euvmod,rm,nespeuv,tx,iz,zenit) 102 else if (solvarmod.eq.1) then 103 call jthermcalc_e107(ig,nlayer,euvmod, 100 call jthermcalc_e107(ig,nlayer,euvmod, 104 101 . rm,nespeuv,tx,iz,zenit,zday) 105 endif106 102 107 103 !Total photoabsorption coefficient … … 118 114 end do 119 115 120 121 return122 123 116 end 124 117 -
TabularUnified trunk/LMDZ.MARS/libf/aeronomars/jthermcalc_e107.F ¶
r1266 r1684 20 20 21 21 implicit none 22 23 include "callkeys.h" 22 24 23 25 c input and output variables … … 136 138 end do 137 139 138 !E10.7 for the day: linear interpolation to tabulated values 139 realday=mod(zday,669.) 140 if(realday.lt.date_e107(1)) then 140 if (solvarmod==0) then 141 e107=fixed_euv_value 142 else 143 !E10.7 for the day: linear interpolation to tabulated values 144 realday=mod(zday,669.) 145 if(realday.lt.date_e107(1)) then 141 146 e107=e107_tab(1) 142 else if(realday.ge.date_e107(669)) then147 else if(realday.ge.date_e107(669)) then 143 148 e107=e107_tab(669) 144 else if(realday.ge.date_e107(1).and.149 else if(realday.ge.date_e107(1).and. 145 150 $ realday.lt.date_e107(669)) then 146 151 do i=1,668 … … 153 158 endif 154 159 enddo 155 endif 160 endif 161 endif ! of if (solvarmod==0) 156 162 157 163 !Photoabsorption coefficients at TOA as a function of E10.7 -
TabularUnified trunk/LMDZ.MARS/libf/phymars/callkeys.h ¶
r1617 r1684 19 19 & ,dustbin,nltemodel,nircorr,solvarmod,solvaryear 20 20 21 COMMON/callkeys_r/topdustref,s olarcondate,semi,alphan,euveff,&22 & tke_heat_flux,dustrefir 21 COMMON/callkeys_r/topdustref,semi,alphan,euveff, & 22 & tke_heat_flux,dustrefir,fixed_euv_value 23 23 24 24 LOGICAL callrad,calldifv,calladj,callcond,callsoil, & … … 38 38 real semi 39 39 real alphan 40 real solarcondate40 real fixed_euv_value 41 41 real euveff 42 42 real tke_heat_flux -
TabularUnified trunk/LMDZ.MARS/libf/phymars/conf_phys.F ¶
r1655 r1684 704 704 705 705 write(*,*) "Method to include solar variability" 706 write(*,*) "0-> old method (using solarcondate); ",707 & "1-> variability wit E10.7"706 write(*,*) "0-> fixed value of E10.7 (fixed_euv_value); ", 707 & "1-> daily evolution of E10.7 (for given solvaryear)" 708 708 solvarmod=1 709 709 call getin("solvarmod",solvarmod) 710 710 write(*,*) " solvarmod = ",solvarmod 711 711 712 write(*,*) "date for solar flux calculation:", 713 & " (1985 < date < 2002)", 714 $ " (Only used if solvarmod=0)" 715 write(*,*) "(Solar min=1996.4 ave=1993.4 max=1990.6)" 716 solarcondate=1993.4 ! default value 717 call getin("solarcondate",solarcondate) 718 write(*,*) " solarcondate = ",solarcondate 712 write(*,*) "Fixed euv (for solvarmod==0) 10.7 value?" 713 write(*,*) " (min=80 , ave=140, max=320)" 714 fixed_euv_value=140 ! default value 715 call getin("fixed_euv_value",fixed_euv_value) 716 write(*,*) " fixed_euv_value = ",fixed_euv_value 719 717 720 718 write(*,*) "Solar variability as observed for MY: " -
TabularUnified trunk/LMDZ.MARS/libf/phymars/physiq_mod.F ¶
r1660 r1684 484 484 call allocate_param_thermos(nlayer) 485 485 call allocate_param_iono(nlayer,nreact) 486 if(solvarmod.eq.0) call param_read 487 if(solvarmod.eq.1) call param_read_e107 486 call param_read_e107 488 487 endif 489 488 #endif
Note: See TracChangeset
for help on using the changeset viewer.