Changeset 769 for palm/trunk/SCRIPTS/NCL/profiles.ncl
- Timestamp:
- Oct 26, 2011 2:20:45 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/profiles.ncl
r731 r769 1048 1048 st=i 1049 1049 break 1050 else 1051 st=0 1050 1052 end if 1051 1053 end do 1052 if (.not. isvar("st"))then1053 print(" ")1054 print("'start_time_step' = "+ start_time_step +"h is invalid")1055 print(" ")1056 print("Select another 'start_time_step'")1057 print(" ")1058 exit1059 end if1060 1054 1061 1055 ; **************************************************** … … 1092 1086 et=i 1093 1087 break 1088 else 1089 et=0 1094 1090 end if 1095 1091 end do 1096 if (.not. isvar("et"))then 1097 print(" ") 1098 print("'end_time_step' = "+ end_time_step +"h is invalid") 1099 print(" ") 1100 print("Select another 'end_time_step'") 1101 print(" ") 1102 exit 1103 end if 1104 1092 1105 1093 delete(start_time_step) 1106 1094 start_time_step=round(st,3) … … 1121 1109 end if 1122 1110 1123 if (time_stride .G E. no_time) then1124 print(" ") 1125 print("'time_stride' is greater equalthan number of available "+\1111 if (time_stride .GT. no_time) then 1112 print(" ") 1113 print("'time_stride' is greater than number of available "+\ 1126 1114 "time steps,") 1127 1115 print("=> 'time_stride' is set to 1") … … 2008 1996 end if 2009 1997 2010 if (abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 10)then 2011 min_value = abs(0.001*min(data(varn,:,min_z_int:max_z_int))) 2012 max_value = abs(0.001*max(data(varn,:,min_z_int:max_z_int))) 2013 else 2014 if (abs(min(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND. \ 2015 abs(max(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then 2016 min_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2017 max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2018 else 2019 if (abs(max(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND.\ 2020 abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then 2021 min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2022 max_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2023 else 2024 min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2025 max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2026 end if 2027 end if 2028 end if 1998 1999 ;data can contain missing values in case of output of t=0h (inital profiles) 2000 ;where no output is possible 2001 n_not_ismissing = num(.not.ismissing(data(varn,:,:))) 2002 2003 if (n_not_ismissing .GT. 0 ) then 2004 2005 if (abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 10)then 2006 min_value = abs(0.001*min(data(varn,:,min_z_int:max_z_int))) 2007 max_value = abs(0.001*max(data(varn,:,min_z_int:max_z_int))) 2008 else 2009 if (abs(min(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND. \ 2010 abs(max(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then 2011 min_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2012 max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2013 else 2014 if (abs(max(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND.\ 2015 abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then 2016 min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2017 max_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2018 else 2019 min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int))) 2020 max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int))) 2021 end if 2022 end if 2023 end if 2029 2024 2030 if (min(data(varn,:,min_z_int:max_z_int)) .EQ. 0 .AND. \ 2031 max(data(varn,:,min_z_int:max_z_int)) .EQ. 0)then 2032 min_value = 0.1 2033 max_value = 0.1 2034 end if 2025 if (min(data(varn,:,min_z_int:max_z_int)) .EQ. 0 .AND. \ 2026 max(data(varn,:,min_z_int:max_z_int)) .EQ. 0)then 2027 min_value = 0.1 2028 max_value = 0.1 2029 end if 2030 2031 else 2032 print(" ") 2033 print(vNam(varn) + " contains only missing values") 2034 print(" ") 2035 end if 2035 2036 2036 2037 if (over .EQ. 0) then
Note: See TracChangeset
for help on using the changeset viewer.