[Dart-dev] DART/branches Revision: 13165

[dart at ucar.edu]

thoar at ucar.edu
2019-05-07 11:11:04 -0600 (Tue, 07 May 2019)
130
Fixing some issues that the gfortran compiler found.
namelists getting declared after executable statements, that sort of thing.




Modified: DART/branches/mizzi/models/wrf_chem/run_scripts/RUN_PERT_CHEM/EMISS_PERT/perturb_chem_emiss_CORR_RT_MA.f90
===================================================================
--- DART/branches/mizzi/models/wrf_chem/run_scripts/RUN_PERT_CHEM/EMISS_PERT/perturb_chem_emiss_CORR_RT_MA.f90	2019-05-05 03:42:48 UTC (rev 13164)
+++ DART/branches/mizzi/models/wrf_chem/run_scripts/RUN_PERT_CHEM/EMISS_PERT/perturb_chem_emiss_CORR_RT_MA.f90	2019-05-07 17:11:04 UTC (rev 13165)
@@ -292,9 +292,9 @@
 ! Perturb the members
                    print *, 'perturb the chemi EMISSs'
                    do imem=1,num_mem
-                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)"),imem
-                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)"),imem
-                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)"),imem
+                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)") imem
+                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)") imem
+                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)") imem
                       wrfchem_file=trim(wrfchemi)//trim(cmem)
                       call get_WRFCHEM_emiss_data(wrfchem_file,ch_chem_spc(isp),chem_data3d,nx,ny,nz_chem)
                       do i=1,nx
@@ -412,9 +412,9 @@
 ! Perturb the members
                    print *, 'perturb the fire EMISSs'
                    do imem=1,num_mem
-                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)"),imem
-                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)"),imem
-                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)"),imem
+                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)") imem
+                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)") imem
+                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)") imem
                       wrffire_file=trim(wrffirechemi)//trim(cmem)
                       call get_WRFCHEM_emiss_data(wrffire_file,ch_fire_spc(isp),chem_data3d,nx,ny,nz_fire)
                       do i=1,nx
@@ -532,9 +532,9 @@
 ! Perturb the members
                    print *, 'perturb the biog EMISSs'
                    do imem=1,num_mem
-                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)"),imem
-                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)"),imem
-                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)"),imem
+                      if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)") imem
+                      if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)") imem
+                      if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)") imem
                       wrfbiog_file=trim(wrfbiochemi)//trim(cmem)
                       call get_WRFCHEM_emiss_data(wrfbiog_file,ch_biog_spc(isp),chem_data3d,nx,ny,nz_biog)
                       do i=1,nx
@@ -766,9 +766,9 @@
 ! Loop through members to find ensemble mean geo_ht
              geo_ht(:,:,:)=0.
              do imem=1,nmem
-                if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)"),imem
-                if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)"),imem
-                if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)"),imem
+                if(imem.ge.0.and.imem.lt.10) write(cmem,"('.e00',i1)") imem
+                if(imem.ge.10.and.imem.lt.100) write(cmem,"('.e0',i2)") imem
+                if(imem.ge.100.and.imem.lt.1000) write(cmem,"('.e',i3)") imem
 !
 ! open netcdf file
 
@@ -1036,7 +1036,7 @@
              end if
              s = ieor(t(1), t(2))
 !             pid = getpid() + 1099279 ! Add a prime
-             call pxfgetpid(pid,ierr)
+             call pxfgetpid(pid,ierr)  ! pxfgetpid is an Intel routine
              s = ieor(s, pid)
              if (n >= 3) then
                 aseed(1) = t(1) + 36269

Modified: DART/branches/mizzi/observations/AIRNOW/airnow_co_ascii_to_obs.f90
===================================================================
--- DART/branches/mizzi/observations/AIRNOW/airnow_co_ascii_to_obs.f90	2019-05-05 03:42:48 UTC (rev 13164)
+++ DART/branches/mizzi/observations/AIRNOW/airnow_co_ascii_to_obs.f90	2019-05-07 17:11:04 UTC (rev 13165)
@@ -128,20 +128,20 @@
       logical                      :: use_log_co, use_log_o3
 
       real                         :: pi,err_frac,ran1,ran2,zfac
-!
+
+      namelist /create_airnow_obs_nml/year0,month0,day0,hour0,beg_year,beg_mon,beg_day, &
+      beg_hour,beg_min,beg_sec,end_year,end_mon,end_day,end_hour,end_min,end_sec, &
+      file_in,lat_mn,lat_mx,lon_mn,lon_mx,use_log_co,use_log_o3 
+
 !============================================================
 !obs sequence extra variables
 !============================================================
-!
-!============================================================
+
       pi=4.*atan(1.0)
       fac=1.0
       err_frac=0.4
       obs_qc(1)=0.
-      namelist /create_airnow_obs_nml/year0,month0,day0,hour0,beg_year,beg_mon,beg_day, &
-      beg_hour,beg_min,beg_sec,end_year,end_mon,end_day,end_hour,end_min,end_sec, &
-      file_in,lat_mn,lat_mx,lon_mn,lon_mx,use_log_co,use_log_o3
-!
+
       save_greg_sec=-9999                                                 
 ! Record the current time, date, etc. to the logfile                                       
       call initialize_utilities('create_obs_sequence')

Modified: DART/branches/mizzi/observations/AIRNOW/airnow_o3_ascii_to_obs.f90
===================================================================