/[MITgcm]/MITgcm/pkg/fizhi/update_earth_exports.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/update_earth_exports.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by molod, Thu Jun 10 20:53:19 2004 UTC revision 1.9 by molod, Thu Jun 10 21:50:33 2004 UTC
# Line 155  C*************************************** Line 155  C***************************************
155        _RL minval        _RL minval
156        external     minval        external     minval
157    
158        PARAMETER (  ALVDRS  = 0.100 )  ! Albedo of soil         for visible   direct solar radiation.  C Albedo of soil         for visible   direct solar radiation.
159        PARAMETER (  ALIDRS  = 0.200 )  ! Albedo of soil         for infra-red direct solar radiation.        PARAMETER (  ALVDRS  = 0.100 )  
160        PARAMETER (  ALVDRDL = 0.300 )  ! Albedo of light desert for visible   direct solar radiation.  C Albedo of soil         for infra-red direct solar radiation.
161        PARAMETER (  ALIDRDL = 0.350 )  ! Albedo of light desert for infra-red direct solar radiation.        PARAMETER (  ALIDRS  = 0.200 )  
162        PARAMETER (  ALVDRDD = 0.250 )  ! Albedo of dark  desert for visible   direct solar radiation.  C Albedo of light desert for visible   direct solar radiation.
163        PARAMETER (  ALIDRDD = 0.300 )  ! Albedo of dark  desert for infra-red direct solar radiation.        PARAMETER (  ALVDRDL = 0.300 )  
164        PARAMETER (  ALVDRI  = 0.800 )  ! Albedo of ice          for visible   direct solar radiation.  C Albedo of light desert for infra-red direct solar radiation.
165        PARAMETER (  ALIDRI  = 0.800 )  ! Albedo of ice          for infra-red direct solar radiation.        PARAMETER (  ALIDRDL = 0.350 )  
166    C Albedo of dark  desert for visible   direct solar radiation.
167          PARAMETER (  ALVDRDD = 0.250 )  
168    C Albedo of dark  desert for infra-red direct solar radiation.
169          PARAMETER (  ALIDRDD = 0.300 )  
170    C Albedo of ice          for visible   direct solar radiation.
171          PARAMETER (  ALVDRI  = 0.800 )  
172    C Albedo of ice          for infra-red direct solar radiation.
173          PARAMETER (  ALIDRI  = 0.800 )  
174    
175  * --------------------------------------------------------------------------------------------  * --------------------------------------------------------------------------------------------
176    
# Line 621  C**** ---------------------------------- Line 629  C**** ----------------------------------
629  #if f77  #if f77
630  cfpp$ expand (coeff)  cfpp$ expand (coeff)
631  #endif  #endif
 #if f90  
 !DIR$ inline always coeff  
 #endif  
632  #endif  #endif
633    
634        DO 100 I=1,IRUN        DO 100 I=1,IRUN
# Line 1026  c      band 12: 35.7 -  oo  um Line 1031  c      band 12: 35.7 -  oo  um
1031  c  c
1032  c-------------------------------------------------------------------------  c-------------------------------------------------------------------------
1033        data ((emis(i,j),i=1,12),j=1,11) /        data ((emis(i,j),i=1,12),j=1,11) /
1034       &   0.9891, 0.9892, 0.9900, 0.9914, 0.9908, 0.9903, ! evergreen needleleaf  C evergreen needleleaf
1035         &   0.9891, 0.9892, 0.9900, 0.9914, 0.9908, 0.9903,
1036       &   0.9898, 0.9948, 1.0000, 1.0000, 1.0000, 1.0000,       &   0.9898, 0.9948, 1.0000, 1.0000, 1.0000, 1.0000,
1037       &   0.9849, 0.9856, 0.9841, 0.9831, 0.9789, 0.9805, ! deciduous needleleaf  C deciduous needleleaf
1038         &   0.9849, 0.9856, 0.9841, 0.9831, 0.9789, 0.9805,
1039       &   0.9733, 0.9869, 1.0000, 1.0000, 1.0000, 1.0000,       &   0.9733, 0.9869, 1.0000, 1.0000, 1.0000, 1.0000,
1040       &   0.9891, 0.9892, 0.9900, 0.9914, 0.9908, 0.9903, ! evergreen needleleaf  C evergreen needleleaf
1041         &   0.9891, 0.9892, 0.9900, 0.9914, 0.9908, 0.9903,
1042       &   0.9898, 0.9948, 1.0000, 1.0000, 1.0000, 1.0000,       &   0.9898, 0.9948, 1.0000, 1.0000, 1.0000, 1.0000,
1043       &   0.9867, 0.9897, 0.9920, 0.9933, 0.9830, 0.9752, ! grasslands  C grasslands
1044         &   0.9867, 0.9897, 0.9920, 0.9933, 0.9830, 0.9752,
1045       &   0.9853, 0.9928, 1.0000, 1.0000, 1.0000, 1.0000,       &   0.9853, 0.9928, 1.0000, 1.0000, 1.0000, 1.0000,
1046       &   0.9490, 0.9697, 0.9738, 0.9712, 0.9474, 0.9582, ! closed shrublands  C closed shrublands
1047         &   0.9490, 0.9697, 0.9738, 0.9712, 0.9474, 0.9582,
1048       &   0.9663, 0.9747, 0.9836, 0.9836, 0.9836, 0.9836,       &   0.9663, 0.9747, 0.9836, 0.9836, 0.9836, 0.9836,
1049       &   0.9469, 0.9670, 0.9883, 0.9795, 0.9751, 0.9767, ! tundra  C tundra
1050         &   0.9469, 0.9670, 0.9883, 0.9795, 0.9751, 0.9767,
1051       &   0.9920, 0.9888, 0.9888, 0.9888, 0.9888, 0.9888,       &   0.9920, 0.9888, 0.9888, 0.9888, 0.9888, 0.9888,
1052       &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766, ! barren  C barren
1053         &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766,
1054       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,
1055       &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766, ! barren  C barren
1056         &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766,
1057       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,
1058       &   0.9998, 0.9998, 0.9998, 0.9998, 0.9998, 0.9999, ! snow/ice  C snow/ice
1059         &   0.9998, 0.9998, 0.9998, 0.9998, 0.9998, 0.9999,
1060       &   0.9997, 0.9994, 0.9995, 0.9995, 0.9995, 0.9995,       &   0.9997, 0.9994, 0.9995, 0.9995, 0.9995, 0.9995,
1061       &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766, ! barren  C barren
1062         &   0.8353, 0.9163, 0.9342, 0.9229, 0.8354, 0.8766,
1063       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,       &   0.9210, 0.9262, 0.9345, 0.9345, 0.9345, 0.9345,
1064       &   0.9788, 0.9833, 0.9819, 0.9820, 0.9835, 0.9865, ! water  C water
1065         &   0.9788, 0.9833, 0.9819, 0.9820, 0.9835, 0.9865,
1066       &   0.9886, 0.9719, 0.9719, 0.9719, 0.9719, 0.9719/       &   0.9886, 0.9719, 0.9719, 0.9719, 0.9719, 0.9719/
1067    
1068        include 'snwmid.h'        include 'snwmid.h'
# Line 1073  c modify emissivity for snow based on sn Line 1089  c modify emissivity for snow based on sn
1089  c-------------------------------------------------------------  c-------------------------------------------------------------
1090          if(snowdepth (i).gt.0.) then          if(snowdepth (i).gt.0.) then
1091           fac = snowdepth(i) / (snowdepth(i) + snwmid(ityp(i)))           fac = snowdepth(i) / (snowdepth(i) + snwmid(ityp(i)))
1092           newemis(i, 1) = newemis(i, 1) + (((emis( 1,9)+emis( 2,9))/2.) - newemis(i, 1)) * fac           newemis(i, 1) = newemis(i, 1) + (((emis( 1,9)+emis( 2,9))/2.)
1093           newemis(i, 2) = newemis(i, 2) + (((emis( 2,9)+emis( 3,9))/2.) - newemis(i, 2)) * fac       .                                           - newemis(i, 1)) * fac
1094           newemis(i, 3) = newemis(i, 3) + (((emis( 4,9)+emis( 5,9))/2.) - newemis(i, 3)) * fac           newemis(i, 2) = newemis(i, 2) + (((emis( 2,9)+emis( 3,9))/2.)
1095           newemis(i, 4) = newemis(i, 4) +              (emis( 6,9)      - newemis(i, 4)) * fac       .                                           - newemis(i, 2)) * fac
1096           newemis(i, 5) = newemis(i, 5) +              (emis( 7,9)      - newemis(i, 5)) * fac           newemis(i, 3) = newemis(i, 3) + (((emis( 4,9)+emis( 5,9))/2.)
1097           newemis(i, 6) = newemis(i, 6) +              (emis( 8,9)      - newemis(i, 6)) * fac       .                                           - newemis(i, 3)) * fac
1098           newemis(i, 7) = newemis(i, 7) +              (emis( 9,9)      - newemis(i, 7)) * fac           newemis(i, 4) = newemis(i, 4) +              (emis( 6,9)      
1099           newemis(i, 8) = newemis(i, 8) + (((emis(10,9)+emis(11,9))/2.) - newemis(i, 8)) * fac       .                                           - newemis(i, 4)) * fac
1100           newemis(i, 9) = newemis(i, 9) +              (emis(12,9)      - newemis(i, 9)) * fac           newemis(i, 5) = newemis(i, 5) +              (emis( 7,9)      
1101           newemis(i,10) = newemis(i,10) +              (emis( 4,9)      - newemis(i,10)) * fac       .                                           - newemis(i, 5)) * fac
1102             newemis(i, 6) = newemis(i, 6) +              (emis( 8,9)      
1103         .                                           - newemis(i, 6)) * fac
1104             newemis(i, 7) = newemis(i, 7) +              (emis( 9,9)      
1105         .                                           - newemis(i, 7)) * fac
1106             newemis(i, 8) = newemis(i, 8) + (((emis(10,9)+emis(11,9))/2.)
1107         .                                           - newemis(i, 8)) * fac
1108             newemis(i, 9) = newemis(i, 9) +              (emis(12,9)      
1109         .                                           - newemis(i, 9)) * fac
1110             newemis(i,10) = newemis(i,10) +              (emis( 4,9)      
1111         .                                           - newemis(i,10)) * fac
1112          endif          endif
1113    
1114  c open water  c open water
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22